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FINAL  REPORT:  SELECTED  TOPICS  IN  COASTAL  AND  OCEAN  TECHNOLOGY 

GRANT  #  N00014-89J-117 

by 

SCOTT  A.  JENKINS 


This  report  covers  work  in  three  distinct  areas  of  research:  1)  sedimentation  control,  2) 
wave-structure  interaction  and  3)  tactical  oceanography.  These  will  be  considered  separately 
below. 

I.  SEDIMENTATION  CONTROL 

Dredging  of  Navy  harbors  represents  an  escalating  cost  of  ownership,  0  ($25  million/year), 
due  to  environmental  restrictions  imposed  by  the  NEPA  Act  and  the  dwindling  dredge  spoils 
pond  capabilities  of  remaining  naval  facilities.  This  work  sought  to  develop  alternative  measures 
to  dredging  that  would  prevent  or  control  sedimentation.  Both  passive  and  active  measures  were 
studied.  The  mechanics  of  either  approach  involve  elevating  fluid  stresses  on  the  bed  above  the 
ambient  levels  of  the  natural  circulation  systems.  Elevated  shear  stresses  are  shown  in 
Appendices  A  and  B  to  cause  resuspension  of  the  unconsolidated  bed  and  erosion  of  the 

consolidated  bed. 

A)  Passive  Methods 

Shear  stresses  are  found  to  become  elevated  passively  in  tidal  flow  by  the 
vortex  wake  of  a  wing  in  ground  effect  with  the  bed.  In  Appendix  A,  a  theory 
and  companion  experiment  are  given  which  describe  the  resuspension  of  fine 
sediments  by  the  action  of  wings  moored  near  the  bottom  of  an  estuary. 
Calculations  by  matched  inner  and  outer  image  expansions  indicate  that  the  vortex 
system  of  a  wing  can  exert  two  distinct  effects  on  estuarial  sedimentation:  (1)  it 
can  elevate  the  shear  stresses  in  the  neighborhood  of  the  wing  by  a  sufficient 
amount  to  induce  erosion  of  partially  consolidated  layers  of  the  bottom,  and,  (2) 
it  can  promote  vertical  transport  to  either  enhance  or  inhibit  downstream 
deposition.  Deposition  in  the  wake  trail  is  found  to  be  dependent  upon  the 
direction  of  the  lift  force  generated  by  the  wing.  Wings  which  produce  a 
downward  lift  force  deplete  the  region  adjacent  to  the  bottom  of  suspended 
sediment  and  thereby  diminish  downstream  deposition.  Conversely,  wings 
generating  an  upward  lift  force  increase  the  density  of  suspended  sediment  near 
the  bottom,  thus  increasing  the  rate  of  downstream  deposition.  Density 
measurements  behind  prototype  wings  are  consistent  with  these  theoretical 
expectations.  Short  term  changes  in  bottom  contours  are  in  general  agreement 
with  predicted  erosion  in  tlie  inner  domain,  and  with  diminished  accretion  in  the 
outer  domain. 

Short  term  wing  induced  perturbations  to  the  bottom  appear  to  be 
qualitatively  explained  by  entrainment  and  diminished  deposition  in  the  far  wake 


and  by  erosion  in  the  near  wake.  Wings  which  generate  a  downward  lift  force 
induce  erosion  in  both  the  near  and  far  wakes.  Wings  which  generate  an  upward 
lift  force  give  rise  to  erosion  in  the  near  wake  followed  by  accretion  in  the  far 
wake. 

Long  term  wing  induced  perturbations  to  the  bottom  are  qualitatively 
explicable  only  in  terms  of  erosion. 

Both  erosion  and  entrainment  are  maximized  by  maximizing  wing  section 
lift  coefficients  while  in  ground  effect  at  low  wing  Reynolds  numbers. 

There  exists  no  unique  selection  of  mooring  elevation  and  wing  aspect 
ratio  which  simultaneously  maximizes  erosion  and  minimizes  downstream 
deposition.  Erosion  is  maximized  by  maximizing  wing  aspect  ratio  for  a  mooring 
elevation,  b,  equal  to  one  quarter  tire  wings  semi  span,  S.  Deposition  is 
minimized  by  minimizing  wing  aspect  ratio  with  b  =  0.8S. 

Wings  which  generate  a  downward  lift  force  are  more  effective  in  causing 
local  erosion  and  retarding  downstream  deposition  than  wings  generating  an 
upward  lift  force.  Because  either  effect  is  limited  to  a  wake  trail  whose  cross¬ 
stream  dimension  is  of  the  order  of  a  wing  span,  large  areas  of  channel  bottom 
will  require  mass  arrays  of  wings  with  overlapping  wakes  (like  a  flock  of  birds), 
in  order  to  control  shoaling. 

Active  Methods 

The  research  also  considered  the  intervention  of  Man-made  horsepower  to 
elevate  bottom  shear  stresses  to  induce  resuspension  and  erosion.  The  results  in 
Appendix  B  indicate  that  deposition  of  fine  grain  sediments  may  be  controlled 
through  the  semi-diurnal  application  of  hydrostatic  jets  producing  applied  shear 
stresses  in  excess  of  the  cohesive  yield  stress  of  the  fluid  mud  layer.  This  may 
be  accomplished  by  either  of  two  means:  (1)  application  of  high  velocity  jets  or 
(2)  application  of  large  diameter  jets.  By  either  approach  the  effect  of  controlling 
deposition  through  jet  entrainment  in  the  far  field  is  maximized  by  maximizing  the 
total  jet  flow  rate.  This  suggests  that  the  total  eddy  momentum  content  of  the  jet 
discharge  is  progressively  diluted  with  increasing  distance  from  the  jet  discharge 
point.  At  some  maximum  range,  the  eddy  momentum  of  the  jet  discharge  is 
diluted  to  such  an  extent  that  it  is  incapable  of  entraining  sediment  from  the  fluid 
mud  layer  across  the  lutocline  interface.  The  unexpected  persistence  of  this 
entrainment  in  the  far  field  of  low  salinity  jets  also  suggests  that  the  jet  salinity 
may  exert  a  modifying  effect  on  the  cohesive  yield  stress  of  the  fluid  mud, 

rendering  it  more  susceptible  to  eddy  transport  across  the  lutocline.  These  results _ 

are  directly  applicable  to  the  problem  of  bottom  resuspensions  by  ship  thrust _ 

wakes.  ^ 

One  of  the  new  technologies  to  result  from  this  work  was  a  novel  method  □ 

for  generating  a  near  bottom  jet  flow.  This  new  method  is  the  turbine-chimney  □ 

and  is  described  in  the  U.S.  Patent  contained  in  Appendix  C.  The  turbine  - - 

chimney  consists  of  a  propeller  mounted  inside  a  vertically  oriented  90°  elbow 

with  a  horizontal  discharge  at  the  bottom.  This  configuration  protects  the - 

propeller  from  bottom  debris  or  slumping  mud  bands,  as  commonly  found  beneath _ 

piers  and  docks.  The  elbow  may  be  hinged  to  a  frame  and  articulated  with  a  -  Codes 
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hydraulic  ram  to  allow  the  jet  to  scan  through  180°  arcs,  thereby  increasing  the 
area  of  influence  of  the  jet  discharge.  Whereas  the  pump  based  systems  generally 
produce  a  high  velocity  low  flow  rate  jet  through  a  small  convergent  nozzle,  the 
propeller  based  system  produces  a  relatively  low  velocity  jet  at  a  very  high  flow 

rate.  • 

Another  concept  to  emerge  from  this  work  is  the  crater-sink  fluidization 

dam  bypassing  concept  described  in  Appendix  D.  This  system  is  intended^  to 
recover  sediments  trapped  behind  inland  dams,  and  to  bypass  them  to  eroding 
sections  of  coastline. 

II  WAVE  STRUCTURE  INTERACTION:  THE  TETRAHEDRON 

SANDBAG  BREAKWATER 

This  study  produced  a  new  approach  to  development  of  mobile  breakwater  systems  and 
is  detailed  in  Appendix  E.  Rather  than  floating  structures,  a  gravity  inertia  structure  is  developed 
which  rests  on  a  broad  footprint  on  the  sea  floor.  Hence  no  connectors  or  moorings  are  required. 
The  stracture  is  modular,  consisting  of  a  pair  of  staggered  piles  of  tetrahedron  sandbags  arranged 
obliquely  to  the  incident  waves  as  shown  in  Appendix  E.  Each  tetrahedron  in  the  array  consists 
of  a  triangular  shaped  sandbag  constructed  of  geotextile  fabric  and  filled  with  sand. 

During  transport  to  the  LOTS  site  the  sandbags  are  readily  stored  aboard  LOTS  ships  in 
a  flat,  folded  and  unfilled  condition.  Hence  many  such  sandbags  can  be  stored  on  a  single  LOTS 
ship  without  the  need  for  occupying  significant  shipboard  space.  Upon  arrival  at  the  LOTS  site 
and  prior  to  unloading  of  supplies  and  material,  each  LOTS  ship  can  deploy  its  own  array  of 
protective  tetrahedron  sandbag  arrays  using  available  shipboard  firepumps  and  crane  support. 
This  is  accomplished  by  attaching  a  small  jet  pump  to  pump  sand  slurries  from  under  the  ship 
into  each  sandbag.  The  tetrahedron  sandbags  are  positioned  around  the  anchorage  using  crane 
support  aboard  the  LOTS  ships.  The  geotextile  fabric  constmction  of  each  tetrahedron  sandbag 
is  strong,  compliant  and  allows  the  water  to  ventilate  out  of  the  bag  during  filling  with  the  sand 
slurry.  If  necessary,  the  bag  might  later  be  emptied  by  means  of  a  Velcro  sealed  outlet  door  near 
the  base  of  the  tetrahedron  sandbag,  or  by  crane  and  clamshell. 

The  theory  of  attenuating  incident  waves  by  means  of  an  array  of  tetrahedron  sandbags 
is  based  on  the  principle  of  partial  Bragg  reflection,  or  Bragg  scattering,  see  Appendix  E. 
Whereas  a  contiguous  saw-toothed  geometry  can  achieve  complete  quarter  wavelength 
cancellation,  the  tetrahedron  array  has  a  certain  amount  of  leakage  between  adjacent  corner 
reflectors.  However,  it  was  found  that  this  leakage  can  be  minimized  for  particular  dimensioning 
and  spacing  of  the  tetrahedrons  based  on  the  incident  frequency  and  direction.  This  can  be 
readily  shown  for  monochromatic  waves  using  the  classical  expansion  in  solid  harmonics  for 
multiple  scattering  centers  (in  this  case  tetrahedrons).  Such  monochromatic  multiple  scattering 
solutions  are  not  new  for  plane,  shallow  water  waves,  except  for  the  triangular  geometry  which 
apparently  does  not  yet  have  a  published  calculation. 

The  new  aspect  of  this  theoretical  approach  for  the  multiple  scattering  problem  is  to 
include  the  effects  of  wave  groups  incident  on  finite  arrays.  This  approach  was  based  on  a 
multiple  scales  perturbation  method  in  which  the  predominant  carrier  swell  waves  are  scaled  in 


terms  of  fast  variables  while  their  groups  are  scaled  in  terms  of  slow  variables  which  are  order 
l/€  times  those  of  the  fast  variables,  where  e  is  the  wave  steepness.  This  scheme  yields 
analytical  solutions  which  differ  from  those  of  the  monochromatic  case  by  the  presence  of  a 
second  order  low  frequency  oscillation  among  the  elements  of  the  tetrahedron  array.  Because 
the  modulations  induced  by  this  second  order  low  frequency  oscillation  are  locked  with  the  wave 
groups  and  propagate  at  the  speed  of  the  group  velocity,  they  appear  seaward  of  the  array  as 
locked  long  waves  in  the  reflected  wave  group,  but  not  in  the  radiated  waves  shoreward  of  the 
array  Inside  the  array  these  locked  long  waves  introduce  a  discontinuity  which  forces  a  second 
order  free  long  wave  propagating  at  a  speed  of  / gh.  The  free  long  waves  appear  shoreward  of 
the  array  and  may  be  resonant  in  a  low  frequency  range  which  can  be  exploited  for  wave  energy 
collection  as  described  in  Appendix  E. 

The  same  discontinuities  arising  from  the  locked  long  waves  among  the  tetrahedrons 
within  the  array  will  also  effect  the  diffraction  around  the  flanks  of  the  array.  This  diffraction 
is  controlled  by  a  second  order  evanescent  mode  decaying  rapidly  with  distance  away  from  the 
array.  The  evanescent  mode  leads  to  various  patterns  of  shadows  and  bright  spots  on  the  flanks 
of  the  array  depending  upon  the  direction  of  the  incident  wave  groups. 

Ill  TACTICAL  OCEANOGRAPHY:  DEVELOPMENT  OF  THE  GLIDING 

SUBMARINE  TACTIC 

This  effort  was  an  offshoot  of  participation  in  the  Tactical  Oceanography  Symposium, 
1990, 1991  sponsored  by  OP-96.  The  results  were  outlined  in  a  report  to  COMSUBDEVRON-12 
contained  in  Appendix  F,  and  were  published  in  unclassified  form  in  the  reviewed  literature  as 
contained  in  Appendix  G.  The  basic  concept  is  one  of  a  means  for  submarines  to  make  low 
speed  translation  between  any  two  points  without  turning  the  screw.  It  is  believed  that  this  would 
be  possible  by  ballasting  in  a  slightly  positive  or  negative  state  and  then  gliding  on  sail-plane  lift 
in  the  manner  of  an  airplane  or  bird.  The  notion  here  is  that  passive  acoustic  detectability  can 
be  reduced  by  trading  the  continual  emittance  of  hydrodynamic  noise  from  propeller  swirl  for  an 
occasional  "burp"  of  ballast  pump  noise  as  the  buoyancy  is  alternately  adjusted  between  positive 
and  negative  states  during  a  dolphin  style  glide  path.  Tactical  advantages  would  likely  arise 
following  high  speed  transits  to  some  patrol  area,  whence  low  speed  dolphin  glides  could  be 
invoked  to  quietly  slip  away  from  suspected  contact. 
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A  theory  and  tnmpanhui  experiment  are  given  winch  describe  the  resuspension  of  line  sediments 
by  the  action  of  wings  moored  near  the  bottom  of  an  estuary.  Calculations  liy  matched  inner 
and  outer  image  expan.sions  indicate  that  the  vortex  .system  of  a  wing  can  exert  two  distinct 
effects  on  estuarial  sedimentation:  ( 1  Mt  can  elevate  the  shear  stresses  in  the  neighborhood  of 
the  wing  liy  a  sufficient  amount  to  induce  erosion  of  partially  consolidated  layers  of  the  bottom: 
and.  (2>  it  can  promote  vertical  transport  to  either  enhance  or  inhibit  downstream  deposition 
Deposition  in  the  wake  trail  is  found  to  be  dependent  upon  the  direction  of  the  lift  force  gen¬ 
erated  bv  the  wing  Wings  which  produce  a  downward  lift  force  deplete  the  region  adjacent  to 
the  bottom  of  .suspended  sediment  and  thereby  diminish  downstream  deposition  Conversely, 
wings  generating  an  upward  lift  force  increase  the  density  of  suspended  sediment  near  the  hot- 
loin,  thus  increasing  the  rate  of  downstream  deposition  Optima!  wing  configurations  are  for- 
niulaleil  winch  maximi/o  erosion  aim  miiiiliu/.i'  Ui-poMtiuii  DeiisUv  measurement?,  beiiiml  pro- 
lot  vfh*  wings  are  consistent  with  these  theoretical  expi*ctations.  Short  term  changes  in  bottom 
contours  are  in  general  agreement  with  predicted  erosion  in  the  inner  domain,  and  with  dimin¬ 
ished  accretion  in  the  outer  domain 
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consolidation,  accretion,  fluid  mud.  entrainment,  liitocline.  Hinpham  fluid,  diffusivity.  scour, 
liftinp  line,  irnape  method,  circulation,  horseshoe-vortices. 


INTRODUCTION 

The  phenomenon  of  current  scour  due  to  vor¬ 
tical  wakes  shed  from  the  support  piles  of 
bridges  and  other  waterfront  structures  is  well 
known,  (COLLINS.  1980,  and  CHIEW  and 
MELVILLE,  1987).  Transport  of  suspended  par¬ 
ticles  by  near-wake  vortices  behind  bluff  bodies 
has  been  found  to  be  especially  vigorous  in  air, 
(MacLENNAN  and  VINCENT,  1982).  Vortices 
from  boat  traffic  have  been  observed  to  scour 
and  resuspend  fine  sediments,  (GARRAD  and 
HEY.  1987).  The  work  herein  begins  wdth  the 
hypothesis  that  the  vortex  system  of  a  moored 
wing  might  perturb  sediment  transport  and 
cause  similar  scour  features  in  the  muddy  bot¬ 
tom  under  steady  or  slowly  varying  estuarine 
flow’.  If  w'ing  induced  scour  is  found  to  he  suffi¬ 
ciently  vigorous,  it  could  emerge  as  a  new'  con¬ 
cept  among  the  paucity  of  passive  technologies 
to  maintain  dredged  harbors  and  w^aterw^ays. 


Presently,  active  technologies  amass  a  $2  bil¬ 
lion  w'orldw'ide  annual  maintenance  dredging 
bill,  (MARINE  BOARD,  1983).  Those  costs  are 
largely  driven  by  dredge  spoils  disposal  prob¬ 
lems  which  would  otherwise  not  exist  with  pas¬ 
sive  methods. 

The  w'ing  vortex  system  appears  to  be  espe¬ 
cially  well  suited  to  influence  sediment  trans¬ 
port  because  it  imparts  a  net  vertical  advection 
to  the  flow  in  reaction  to  the  lift  force,  (LAN- 
CHESTER.  1908.  PRANDTL  and  BETZ,  1927, 
and  PRANDTL,  1931).  This  vertical  advection 
(induced  velocity)  might  be  employed  to  either 
reduce  or  increase  the  net  settling  rates  of  silts 
and  clays,  depending  upon  the  direction  of  the 
lift  force,  (JENKINS  and  SPARKS,  1985,  and 
JENKINS,  1987).  Because  the  preponderance  of 
transport  of  these  fine-grained  sediments 
occurs  near  the  bottom.  (EINSTEIN  and 
KRONE.  1961,  and  MEHTA,  et  a/.,  1989),  we 
are  particularly  interested  in  the  case  of  a  wing 
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in  proximity  to  the  iminolnle  bed.  /  in  lm'uuikI 
efTcct.  (LAGALLY,  1929,  PRESMAK,  1977.  and 
McC'OKMICK.  1979!.  (iround  effect  is  generally 
considered  to  be  the  condition  in  which  the  wing 
i.<  within  one  semi-span  length  of  a  solid  plan(‘ 
boundary. 

We  seek  analytic  solutions  to  the  wing 
induced  perturbation.--  on  sediment  tran-p(;ri 
and  th(*  bed  response  so  that  cause  and  effect 
relations  might  be  readily  established.  There¬ 
fore  it  shall  be  necessary  to  make  a  number  of 
idealizations  with  respect  to  both  the  sediment 
transport  physics  as  well  as  the  flow  field  pro¬ 
duced  by  the  wing.  We  employ  a  PYckian-hased 
advection-diffusion  transport  model  to  provide 
a  deterministic  approach  with  sufficient  kine¬ 
matic  detail  to  allow  analytic  progress  on  the 
problem  at  hand.  These  are  the  most  widely 
adopted  fine-sediment  transport  models  in  the 
published  literature  and  have  been  validated 
bv  both  laboratory  data,  i  APMAN  and  RUMEK, 

1 970),  and  field  data,  i  KP]RSSENS  ct  ai.,  1 979  i. 
I'hey  have  been  employed  to  account  for  trans¬ 
port  over  eroding  or  accreting  bed.--, 

MONT  and  LEBON.  1978,  and  CHENG.  1985i, 
and  for  the  vertical  and  streamwise  variabiiit\ 
m  transport,  (ME I,  1969).  The  latter  was  con¬ 
strained  by  the  assumption  that  the  settling 
velocity  of  the  sediment,  (o„,  is  of  the  same  order 
as  the  mean  flow  velocity,  u,..  This  rather 
unrealistic  assumption  for  the  case  of  fine-sed- 
irnent  transport  has  been  circumvented  herein 
by  coordinate  stretching  based  on  the  strength 
of  the  wing  perturbation. 

With  regard  to  specifying  the  flow  field 
around  a  wing  moored  near  an  estuarine  bot¬ 
tom,  there  are  two  major  difficulties:  (1  '  viscos¬ 
ity:  and  (2)  density  stratification  due  to  the 
abundance  of  suspended  sediment  near  the  bot¬ 
tom.  Unfortunately,  the  recent  wung  theories 
which  explicitly  account  for  the  effects  of  either 
viscosity  or  density  stratification  are  numeri¬ 
cal,  MOORE  (1974),  MATTEl  and  SAN¬ 
TORO  (1974),  STAUFENBIEL  and  KLEINER- 
UAM  (1980),  PULLIN  and  PHILLIPS  *  1981 )  or 
SAP^FMAN  (1972),  HILL  (1975)  and  GREEN 
1985),  respectively.  Analytic  solutions  lor  the 
tiow  field  induced  by  a  wing  of  finite  span  pre¬ 
sent  one  of  the  most  intractable  problems  in 
aerodynamic  theory,  (VAN  DYKPL  1975'. 
Therefore  w'e  shall  neglect  stratification  and 
approximate  the  local  flow  field  as  an  inviscid 
singular  perturbation  in  the  wing  aspect  ratio. 


A^,  alter  V  AN  DYKli  il9()4i.  This  is  based  on 
the  premise  that  the  wung  circulation  is  large 
compared  to  ambient  levels  and  that  the  rota¬ 
tional  and  irrotational  parts  of  the  flow  do  not 
interact.  Difficulties  associated  with  the  no-slip 
condition  at  the  bottom  do  not  explicitly  arise 
herein  because  the  sediment-transport  equa- 
iion.s  dtn  eioped  in  the  following  sections  rely  on 
inviscid  wing-theory  only  to  specify  the  vertical 
component  of  the  flow  field.  Details  of  the  trail¬ 
ing  vortex  roll-up  in  the  near  wake  of  the  wing 
provided  by  the  analytic  results  of  BATCHPL 
LOR  (1964),  SQUIRE  (1965),  GOVINDARAJU 
and  SAFEMAN  (1972!.  SAFFMAN  (1979*  and 
PH1L17IPS  (1981  )  do  not  appear  at  the  order  of 
approximation  taken  in  the  formulation  of  the 
sediment-transport  equations  below. 

Formulation 

Considm'  a  fluid  of  density,  .  transporting  a 
suspension  of  uniform  particles,  each  with  a 
settling  velocity  relative'  to  the  fluid,  co,,.  and  a 
solid  (iensitx .  (»..  2.65  g  em  .  Tlu'  density 

of'  the  fluid-sediment  mixture,  may  be 

written; 

p,„  (),,N  ^  (l-N)p^  (1) 

where  N  is  the  volume  concentration  equal  to 
the  volume  of  sediment  per  volume  of  the  fluid- 
sediment  mixture.  The  density  of  the  sediment 
component  of  this  mixture,  p,  is  sometimes 
referred  to  as  the  excess  or  bulk  density  and  is 
defined: 

|)  (>,|N  (2) 

Assume  a  uniform  channel  flow  with  a  char¬ 
acteristic  velocity  u,,  which  varies  slowly  at 
angular  frequency  il„  over  a  flat  cohesive  sur¬ 
face  on  the  plane  /  0,  see  Figure  1.  At  or 

below  z  0.  the  fluid-sediment  mixture 
exhibits  a  non-zero  cohesive  yield  stress,  t,,  and 
behaves  as  a  viscoplastic.  Above  z  -  0,  the  co¬ 
hesive  yield  stress  vanishes  and  the  fluid-sedi¬ 
ment  mixture  behaves  as  a  fluid,  i.e.,  unable  to 
support  shear  stress  at  equilibrium.  The  sedi¬ 
ment  component  den.^ity  at  z  0  w^here  the 
mixture  first  begins  to  exhibit  a  cohesive  yield 
stress  IS,  ()  . 

The  upper  regions  of  the  fluid,  z  *  h,  arc  rel¬ 
atively  deficient  of  suspended  sediment  and 
ma\  f)e  characterized  by  some  small  ambient 
value,  (>  p,.  Because  of  settling  under  the 
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Figure  1.  Conceptual  diagram  of  the  fluid  bed  interface  with  the  wing  and  its  mooring  details. 


influence  of  gravity,  the  abundance  of  sus¬ 
pended  sediment  begins  to  increase  below  z  = 
8,  forming  a  non-cohesive  floe  layer  as  shown  in 
Figure  1,  (MEHTA  and  PARTHENIADES, 
1975,  and  KRONE,  1978).  The  surface  z  =  0 
may  thus  be  considered  a  “lutocline"’,  (DYER, 
1985,  and  MEHTA,  1989), 

Immediately  below  the  lutocline,  there  is  a 


rather  thin  layer  of  cohesive  fluid  mud  between, 
0  ^  z  <  -m,  which  has  a  non-zero  velocity, 
(KIRBY  and  PARKER,  1974,  PARKER  and 
KIRBY,  1982  and  NICHOLS,  1985).  Thus  the 
lutocline  at  z  =  0  between  the  fluid  and  the  vis¬ 
coplastic  material  is  not  a  no-slip  boundary. 
The  dynamics  of  the  fluid  mud  layer  are  still 
poorly  understood  when  subjected  to  a  moving 
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fluid.  For  analytic  purposes,  we  shall  assume 
that  the  lutocline  at  z  =  0  is  a  stress  free  bound¬ 
ary  much  like  a  free  surface,  with  the  shear 
stresses  resulting  from  the  fluid  motion  appear¬ 
ing  in  the  stress  field  of  the  fluid  mud  layer. 
This  assumption  relieves  the  wing  induced 
velocity  field  on  the  fluid  side  of  the  lutocline 
inlerfaL-e  from  having  to  satisfy  the  no-slip  con¬ 
dition  at  z  =  0.  The  no-slip  boundary  thus 
appears  at  z  =  -m  where  the  immobile,  par¬ 
tially  consolidated  bottom  begins. 

Below  the  fluid  mud  layer,  z  <  -m,  the  sed¬ 
iment  is  at  rest  and  there  is  a  progressive  build 
up  in  sediment  component  density  to  some  ulti¬ 
mate  saturation  value,  p  =  p„.  This  is  due  to 
consolidation  and  compaction  under  the  weight 
of  the  overburden  which  proceeds  at  a  rate  char¬ 
acterized  by  pcK.g,  where  g  is  the  acceleration 
of  gravity  and  K,  is  the  sedimentation  coeffi¬ 
cient  in  seconds,  (FUJITA,  1962).  The  portion 
of  the  cohesive  bottom  below  the  fluid  mud  layer 
for  which  p^  ^  p  <  p«  shall  be  referred  to  as  par¬ 
tially  consolidated  mud,  while  the  fully  consol¬ 
idated  mud  in  Figure  1  represents  the  oldest 
deposition  which  has  reached  final  saturation 
density  p  =  p,. 

Now  superimpose  a  wing  on  the  river  flow  at 
a  distance  z  =  b  above  the  lutocline.  Take  the 
wing  semi-span  to  be  S.  Assume  the  wing  to  be 
in  ground  effect  so  that  b  <  S.  Let  the  mean 
semi-cord  be  c  so  that  the  wing  aspect  ratio  is 
A,^  =  S/c.  We  shall  prescribe  the  sediment 
transport  and  the  wing  induced  velocity  field  in 
an  Eulerian  frame  with  respect  to  Cartesian 
coordinates  (x,y,z).  The  origin  of  this  system  is 
located  at  the  intersection  between  the  lutoc¬ 
line  plane  and  a  vertical  projection  from  the 
mean  aerodynamic  center  of  the  wing,  see  Fig¬ 
ure  1.  The  axis  o-z  shall  be  positive  vertically 
upward;  the  axis  o-x  is  positive  in  the  direction 
of  the  river  channel  flow;  while  the  axis  o-y  is 
transverse  to  the  mean  channel  flow  in  a  right 
hand  system. 

With  the  problem  posed  in  this  w^ay,  the  only 
horizontal  variability  will  be  that  imposed  by 
the  wing  over  length  scales  the  size  of  a  wing 
semi-span.  On  the  other  hand,  the  density  field 
is  structured  vertically,  with  strong  variability 
near  the  lutocline  over  length  scales  of  the 
order  of  the  floe  layer  thickness.  Therefore,  we 
prescribe  the  following  scalings  for  the  dimen¬ 
sional  variables: 


X  =  xS 
y  =  yS 

z  =  z6  =  zS/Ajr 
t  =  t/a, 

(u,v,w)  =  (u,v,w)Uo 
n  =  p.p 


(3) 


where  carets  denote  dimensionless  variables 
and  Af  =  S/8  is  the  density  field  aspect  ratio. 

Sediment  transport  by  the  fluid  velocity  field, 
(u,  V,  w),  can  then  be  expressed  by  the  non- 
dimensional  Fickian-based  advection-diffusion 
equation  as  follows: 


6p 

€ 


^  R 


.  0Jo\dp  ,5w 

^  +  r  ^  +  p « 


Rr 

Ap 


u, 


-H  U“T  + 

dx  dx 


dz 


dz 
dv 

^  dyl 


+ 


Ap  \cix^ 


-f- 


dyV 


(4) 


where  R,.  =  u„5/€  is  the  mass  flux  Reynolds 
number  equal  to  a  ratio  of  advective  to  diffusive 
mass  fluxes  and  e  is  the  mass  diffusivity.  Under 
the  prototype  conditions  typical  of  the  lower 
depths  of  estuaries  along  the  U.S.  west  coast  as 
encountered  in  Section  7,  we  take: 


Uo  =  0(10  cm/sec) 

(u„  =  0(10"^  cm/sec) 
a  =  0(10  "  sec^') 
h  =  0(10"  cm  -  10"  cm)  (6) 

5  =  0(10  cm) 

€  =  0(1  cm"/sec) 
c  »  0(8) 


The  final  order  of  magnitude  argument  in  (5) 
sets  a  minimum  size  restriction  on  the  wing. 
This  is  necessary  to  avert  Hele-Shaw  flow, 
fWERLE,  1973  and  McMASTERS,  1974)  and  to 
assure  that  the  wing  circulation  is  large  com¬ 
pared  to  that  of  the  ambient  near-bottom  flow. 


OUTER  SOLUTION 


In  the  outer  domain,  x  >  0(1),  the  lowest 
order  solution  follows  from  (4)  when  the  chan¬ 
nel  flow  is  subjected  to  an  0(1/A^)  perturbation 
from  a  wing  at  a  small  angle  of  attack  p.  The 
perturbation  from  a  wing  in  ground  effect  at  a 
distance  z  =  b  above  a  plane  surface  was  first 
represented  by  an  image  pair  of  lifting  lines  by 
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BETZ  (1912).  The  velocity  potential  of  this  fa- 
milar  vortex  system  is  written  in  terms  of  outer 
variables  (3)  as: 


4)  =  X  +  3 


3f(x)  f’  _ 

2A„  J-i  /  ; 


Is  -  > 


£  -  I)  +  (y  -  0" 


-  Hi 

li  b 

lAr  s 


+  (y  -  0“ 


d? 


T(p 

_  f* 


>/*'  ^  (i  ^  I 


+  (y  -  If 


=  dc 


(6) 


where  TC^)  is  a  unit  function  defining  the  span- 
wise  distribution  of  circulation;  fix)  is  a  unit 
function  defining  the  rate  of  downstream  decay 
of  circulation  in  the  trailing  vortices  and  J  is  a 
dummy  variable.  Terms  containing  (z  +  b)  cor¬ 
respond  to  the  image  vortex  system  at  z  =  -b. 
We  have  adopted  herein  the  convention  that  the 
lift  force  due  to  the  real  vortex  acts  upward  in 
the  positive  z  direction  for  p  >  0. 

Equations  (5)  and  (6)  require  a  large,  high 
aspect  ratio  wing,  whence  Af  A^.  «  1.  Conse¬ 
quently,  the  advection-diffusion  equation  (4) 
may  be  written  at  lowest  order: 

+  =  o  (7) 

Aw  \  Hz  Hz  J 

where  w^^’  =  kyH^IHz  is  the  0(1/A,,)  induced 
velocity.  Equation  (7)  is  subject  to  boundary 
conditions  which  require: 

p'®’  =  1  at  z  =  0  (8) 

p‘°’  =  Pt/pr  at  z  =  h/8  (9) 

With  (6)  a  solution  to  equation  (7)  can  be 
obtained  by  successive  integrations,  evaluating 


the  integration  constants  according  to  the 
boundary  conditions  (8)  and  (9)  to  yield  the  fol¬ 
lowing  result  in  dimensional  form: 


p  -  p,  exp|^[Ii(x,y,z)  -  li(0) 


-1.(0: 


-I.(h)] 


I  log^l  (10' 

h  PcJ 


where  1,(0)  =  Ii(x,y,z  =  0);  l,(h)  =  Ii(x,y,z  = 
h)  and  (see  equation  11). 

Here  I'o  is  the  circulation  of  the  bound  vortex 
of  the  wing  and  T  is  the  circulation  of  the  trail¬ 
ing  vortices  such  that, 


r„  =  UoC,c 

r  =  Tofix) 


(12) 

(13) 


C,  in  equation  (12)  is  the  lift  coefficient  which 
shall  be  positive  for  a  wing  above  the  bottom 
which  generates  lift  upward  in  the  positive  z- 
direction.  This  case  shall  be  referred  to  herein 
as  a  “downwashing”  wing.  When  C,  <  0,  the  lift 
is  directed  downward  toward  the  bottom  and 
these  shall  be  referred  to  as  “upwashing” 
wings.  The  decay  rate,  fix),  for  the  circulation 
of  the  trailing  vortices  in  (13)  is  due  to  the 
ground  effect  calculations  of  PEACE  and 
RILEY  (1983)  for  the  case  of  a  stress  free  bound¬ 
ary  at  z  =  0.  Their  results  by  numerical  inte¬ 
gration  give; 

.  ,1/2 

fe)  .  1  -  f  {«'»  +  (g) 


dx 


=  1  -t  (5.410  X  10-*) 


A.S 

xC, 

xC| 

,A„S 


(14) 


-  (6.670  X  10-*) 


+  (1.173  X  10'*) 


-  (9.241  X  10-®)h-5 


+  (2.729  X  10-'“) 


where  o  is  the  diffusivity  of  momentum. 
INNER  SOLUTION 

The  outer  solution  (10)  develops  singular 
behavior  as  x  -►  0  due  to  the  line  singularity  of 


Journal  of  Coastal  Research,  Vol.  6,  No.  4,  1990 


966 


Jenkins  and  Wasyl 


.  Tof.  r  (z  -  W(y  +  S)  1  _  (z  +  bKy-i-S) _ H 

l,(x.y.z.  = -jtan  tan  ^ 

_  h  L„„-.  [  ^  t»(y  -  S) _ 1  -  tan-  [  <^-b)(y-S)  .....1 

[  ‘ -  I  |x\  x2  ^  (y  -  R>2  -  (7  -  b)2| 

+  (y  ~  S)^  +  (z  +  b)2 


1  -  tan-  [ -  1! 

'  -  S)2  +  (Z  +  b)2j  L<y  -  S)  Vx2  +  (y  -  S)2  +  (Z  -  b)* J  J 


--!l|tan-[— — -= _ 

477  [  L(y  -  S)  Vx2  +  (y 

-  i  (tan-  - ^  ^  . -  -])  +  f  (“"“[— 

477  1  L(y  +  S)  Vx2  +  (y  +  s)2  +  (z  -  h)^ii  I  Lty  + 

Equation  11. 


x(z  +  b) 


S)  Vx2  +  (y  +  S)2  +  (z  +  b)2 


the  bound  vortex  which  approximates  the  wing 
itself  in  the  outer  domain.  This  non-uniformity 
in  the  lifting  line  approximation  is  well  known, 
(VAN  DYKE,  1964),  and  can  be  remedied  as  in 
aerodynamic  theory  by  introducing  minified 
dimensionless  inner  variables  (X,  Y,  Z)  such 
that: 


dZ" 


A. 


dZ 


+  p- 


dZ 


=  0  (18) 


where  =  Afr/A^di>/dZ  is  the  0(1)  induced 
velocity.  On  (18)  the  following  boundary  con¬ 
ditions  are  imposed: 


X  =  A  Jc 

Y  =  y  (15) 

Z  =  A^z 


The  local  flow  field  in  terms  of  (15)  then 
reduces  to  plane  flow  past  an  image  pair  of  flat 
plates  at  angle  of  attack^P  as  defined  by  equa¬ 
tion  16,  where  =  A;*d)  and  /„  is  the  imagi¬ 
nary  part  operator.  The  non-dimensional  form 
of  the  advection-diffusion  equation  in  the  inner 
domain  can  now  be  expressed  in  terms  of  inner 
variables  as: 


0  = 


a»p 

az* 


a„L\  uo/az 


Rp 

AwAp 

Rp 

A^P 

1 


aY 


.au 

.av 

p— *- 

dY 


^  dY^ 


(17) 


Inserting  (16)  in  (17),  the  leading  order  trans¬ 
port  equation  in  the  inner  domain  becomes, 


p<®’  -  1  (g  Z  =  0  (19) 

P<0)  =  Pt/pc  @  Z  =  A,^  h/8  (20) 

Upon  successive  integrations  the  solution  to 
(18)  subject  to  (19)  and  (20)  is  found  in  dimen¬ 
sional  form  to  be: 

p  =  p,  exp|^[l2(x,y,z)  -  l2(0) 

+  fyO)  -  fwh)]  +  ^log2i|  (21) 

h  h  h  pcj 

where  l2(0)  =  l2(x,y,z  =  0);  Mh)  =  l2(x,y,z  = 
h)  and  (see  equation  22). 

It  can  be  demonstrated  at  0(Rf/A,^)  that  the 
inner  solution  by  (21)  and  (22)  matches  asymp¬ 
totically  with  the  outer  solution  given  by  (10) 
and  (11)  asA^-^^  in  the  overlap  domain  where 
X  ~  0(1).  Even  for  relatively  small  aspect  ratios 
the  match  is  quite  good  as  shown  in  Figure  2  for 
the  cases  of  both  downwashing  and  upwashing 
wings.  The  match  remains  satisfactory  in  as 
close  as  x/S  =  2/3,  at  relatively  high  ambient 
velocities,  u^  =  50  cm/sec. 

We  find  in  Figure  3a  that  the  net  effect  in  the 
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4>  =  X  cobB  +  (sinpl'/ms 

'll  -  T*(Y) 

[  ^ 

V  [  \Ak  c 

/J 

^  TlYilog 

■  4  ■ ‘  J 

-  2 

-  rci'i 

[  XAy  Cl  V 

L  Ur  cl 

-1 

+  Xco8{-p)  +  l8in(-p)J ‘/ml  J 

2 

-  T^Y) 

-  T(Y)  log 

■  2 

-  T’^Y) 

Equation  16. 


MATCH  OF  INNER  AND  OUTER  SOLUTIONS 
at  X  =  S  and  Y  -  0 


350*1 

Cl  =  ±0.2  u  =  10  cm/sec 

300  ] 

\, 

downwashing  wing  Aw  =  5: 

•  INNER  SOLUTION 

200  ; 

'  -  OUTER  SOLUTION 

^  200- 

upwashing  wing  Aw  =  5: 

u 

\ 

\ 

\ 

\ 

INNER  SOLUTION 

^  150- 

—  OUTER  SOLUTION 

100- 

50 

(XOOO  0.025  0.050  0.075  0.100  0.125 

p  (  g/cm^ ) 

Figure  2.  Match  of  the  inner  and  outer  solutions  for  prototype 
conditions.  Calculations  are  based  on:  c  »  4.961  cm*/sec,  pi  ■= 
0.0125  g/cm^  and  Pc  =  0.10  g/cm^,  and  h  ■  S  «  304.8  cm. 


wake  of  an  upwashing  wing,  (C/,p)  <  0,  is  to  dis¬ 
place  suspended  sediment  upward  away  from 
the  bottom  creating  a  density  inversion  layer 
above  the  wing.  In  consequence  the  density  gra¬ 
dient  at  the  lutocline  becomes  more  strongly 
negative  relative  to  the  conditions  far  down¬ 
stream.  Conversely,  the  wake  of  the  downwash¬ 
ing  wing  in  Figure  3b  is  shown  to  drive  the  sed¬ 
iment  downward,  depleting  the  region  above  of 
the  floe  layer  and  elevating  the  suspended  con¬ 
centrations  beneath  the  wing  and  directly  adja¬ 
cent  to  the  lutocline.  This  results  in  an  inver¬ 
sion  layer  under  the  wing  which  can  reverse  the 
sign  of  the  density  gradient  at  the  lutocline. 
The  calculations  are  for  a  wing  moored  100  cm 
above  the  lutocline  in  a  10  cm/sec  ambient  cur¬ 


rent.  In  either  case  the  wing  induced  density 
perturbations  die  out  in  the  outer  solution  with 
increasing  distance  downstream  until  x  = 
where  the  density  profile  is  equivalent  to  the 
exponential  no-wing  case.  Figure  4  shows  cross¬ 
stream  variations  in  the  wing  induced  density 
perturbations  that  are  similar  to  those  found  in 
the  stream-wise  direction  in  Figure  3.  The 
cross-stream  variability  is  shown  for  a  tapered 
wing  whose  spanwise  circulation  distribution 
varies  according  to  VAN  DYKE  (1964)  as: 

T(y)  =  Re{a  -  y"/ST"}  (23) 

where  Re  is  the  real  part  operator.  We  find  that 
beyond  the  wing  tips,  lyl  >  S,  the  density  profile 
of  the  inner  solution  again  approaches  the  expo¬ 
nential  no-wing  profile  of  the  undisturbed 
ambient  flow. 

EROSION  AND  DEPOSITION  Shoaling 

of  the  bottom  in  a  time-averaged  sense  is  the 
net  resultant  of  two  processes,  erosion  and  dep¬ 
osition.  Therefore,  shoaling  can  be  retarded  or 
prevented  by  either  of  two  mechanisms:  (1)  by 
eroding  immobile,  partially  consolidated  mud 
below  z  =  -m,  and  (2)  by  reducing  or  prevent¬ 
ing  the  deposition  of  new  mud  above  z  =  -m. 
The  latter  is  accomplished  by  entraining  sus¬ 
pended  material  from  the  fluid  mud  layer  before 
it  can  become  incorporated  on  the  partially  con¬ 
solidated  bottom  at  z  =  m. 

Generally  erosion  results  when  an  applied 
stress  exceeds  the  cohesive  yield  stress,  t^,  of 
the  immobile  bed  below  z  =  —m  (MEHTA,  et 
al.y  1989).  In  our  formulation  of  the  bottom 
interface,  this  can  only  result  from  stress  trans- 
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Equation  22. 


fer  through  the  fluid  mud  layer.  The  velocity  at 
the  top  of  the  fluid  mud  layer  must  match  the 
fluid  velocity  along  the  stress  free  boundary  at 
z  =  0.  However  at  the  bottom  of  the  fluid  mud 
layer,  z  =  -  m,  the  velocity  must  vanish,  giving 
rise  to  shear  stresses  within  the  fluid  mud. 
Because  the  fluid  mud  is  in  motion  these  shear 
stresses,  t,  must  already  be  in  excess  of  the  coh¬ 
esive  yield  stress  of  the  fluid  mud.  Under  the 
assumption  that  the  fluid  mud  behaves  as  a  thin 
Bingham  fluid,  (VANONI,  1975,  FAAS,  1985, 
WRIGHT  and  KRONE,  1987),  this  excess  shear 
stress,  T  -  Tp,  will  vary  with  the  wing  distur¬ 
bance  as: 


where  is  the  apparent  viscosity  of  the  fluid 
mud. 

There  are  a  number  of  erosion  theories  which 
predict  that  the  rate  of  erosion  increases  with 
increasing  t  -  The  exponential  rate  theories 
of  PARTHENIADES  (1965),  CHRISTENSEN 
(1965),  CHRISTENSEN  and  DAS  (1973), 
RAUDKIVI  and  HUTCHISON  (1974)  and 
GULARTE  (1978)  were  found  to  grossly  under¬ 
estimate  the  wing-induced  erosion  fluxes  under 
the  prototype  conditions  described  in  Section 
7.0.  This  is  consistent  with  the  fact  that  these 
rate  equations  are  based  upon  data  from  arti¬ 
ficially  placed  beds  which  characteristically 
erode  more  slowly  than  naturally  deposited 
beds,  (MEHTA,  et  al.,  1982).  The  linear  theories 
of  KANDIAH  (1974),  ARULANANDAN  (1975), 
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DOWNSTREAM  DECAY  of  DENSITY  ANOMAUES 
UPWASHING  WING  at  Y  =  0 
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Figure  3  (A).  Downstream  decay  of  wing  induced  perturba¬ 
tions  to  the  sediment  component  density  as  defined  by  equa¬ 
tion  (2).  Calculations  are  based  on:  <  *=  4.951  cm^/sec,  pt  = 
0.0125  g/cm^  and  Pc  «  0.10  g/cm^,  and  h  S  =  304.8  cm. 

THORN  and  PARSONS  (1980)  and  MEHTA 
(1981)  can  be  calibrated  from  erosion  rates 
observed  over  deposited  beds  because  they  are 
parameterized  by  a  single  empirically  derived 
factor,  a,  according  to: 


Here  a  >  0  and  positive  values  of  erosion  flux, 
E,  correspond  to  erosion. 

A  contour  map  of  the  wing  induced  erosion 
flux  from  (24)  and  (25)  for  C/  =  +  0.2  is  calcu¬ 
lated  in  Figure  5a  for  prototype  conditions  typ¬ 


DOWNSTREAM  DECAY  of  DENSITY  ANOMAUES 
DOWNWASHING  WING  at  Y  *  0 
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Figure  3  (B).  Downstream  decay  of  wing  induced  perturba¬ 
tions  to  the  sediment  component  density  as  defined  by  equa¬ 
tion  (2).  Calculations  are  based  on:  e  =  4.951  cm^/scc,  pi  = 
0.0125  g/cm^  and  pc  =  0.10  g/cm^,  and  h  =  S  =  304.8  cm. 


SPANWISE  EVOLUTION  of  DENSITY  ANOMAUES 
INNER  SOLUTION:  UPWASHING  WING  at  X  =  S/4 


/)(g/cm®) 


Figure  4  (A).  Comparisons  between  theory  and  data  for  span- 
wise  variations  in  the  wing  induced  perturbations  to  the  sed¬ 
iment  component  density  in  the  near  wake  at  X  =  S/4.  Cal¬ 
culations  are  based  on:  e  =  4.951  cm^/sec,  pi  =  0.0126  g/cm* 
and  Pc  =  0.10  g/cm^  and  h  =  S  ~  304.8  cm. 

ical  of  the  experiment  described  in  Section  7.0. 
Horseshoe-like  erosion  features  are  predicted 
in  the  near  wake,  due  to  the  action  of  the  hor¬ 
izontal  induced  velocity  components.  Such  ero¬ 
sion  is  independent  of  the  sign  of  the  lift  coef- 
flcient.  However,  erosion  alone  fails  to  account 
for  shoaling  details  which  might  occur  further 
downstream  in  response  to  an  increased  level  of 
suspended  sediment  from  mobilized  bed  mater¬ 
ial. 

The  deposition  rate  at  z  =  m  is  directly  pro¬ 
portional  to  the  abundance  of  suspended  sedi- 

SPANWISE  EVOLUTION  of  DENSITY  ANOMAUES 
INNER  SOLUTION:  DOWNWASHING  WING  ol  X  =  S/4 


P  (  g/cm^) 


Figure  4  (B).  Comparisons  between  theory  and  data  for  span- 
wise  variations  in  the  wing  induced  perturbations  to  the  sed¬ 
iment  component  density  in  the  near  wake  at  X  =  S/4.  Cal¬ 
culations  are  based  on:  t  =  4.951  cm^/sec,  pt  =  0.0125  g/cm 
and  Pc  =  0.10  g/cm^,  and  h  =  S  ==  304.8  cm. 
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KKOSION  FLUX.  K  10  sec 


Figure  5.  Contour  maps  of  the  erosion  and  deposition  fluxes  induced  by  a  wing  with  Aw  =  5  and  C/  =  +  0.2  at  a  mooring  elevation 
of  b  =  100  cm  in  a  10  cm/sec  mean  flow  (A)  Erosion  flux  for  p  =  3.06  x  10“^  gm/cm  sec;  Tc  =  0.5  dynes/cm^  and  a  =  2.356  x  10  ^ 
g/cm^  sec. 


ment  near  the  bottom,  (ARIATHURAI  and 
KRONE,  1976,  and  COLE  and  MILES,  1983). 
Th  is  abundance  is  determined  by  the  mass 
fluxes  across  the  lutocline  at  z  =  0.  The  flux  of 
suspended  sediment  entering  the  fluid  mud 
layer  due  to  settling  under  gravity  is  -ptoo- 
These  downward  directed  settling  fluxes  are  the 
only  mass  fluxes  considered  by  EINSTEIN  and 
KRONE  (1962),  or  MEHTA  (1989)  when 
accounting  for  deposition  in  the  low  turbulence 
environment  of  a  settling  column  or  laboratory 
flume.  However  in  the  high  turbulence  environ¬ 
ment  of  a  vortical  wake  in  a  natural  estuary, 
vertical  transport  due  to  eddy  diffusion  must 
surely  become  important,  especially  along  the 
strong  density  gradients  in  the  neighborhood  of 
the  lutocline.  The  net  of  settling  and  diffusive 
fluxes  across  the  lutocline  determines  the  rate 
of  change  of  suspended  sediment  in  the  fluid 
mud  layer.  However,  there  is  only  a  finite  prob¬ 
ability,  p(0,l),  that  this  suspended  sediment 
will  actually  stick  to  the  bed  at  z  =  -m  and 
consolidate  to  some  ultimate  saturation  density 
p„  (KRONE,  1962).  Therefore  the  deposition 
flux,  D,  (including  vertical  diffusion  and  con¬ 
solidation),  may  be  written: 


”K.Peg  +  p(0, 


D  = 


dpi 

,1)  *  Pe(Do  - 


(1  -  Pe/pJ 


(26) 


Here  positive  values  for  D  correspond  to  depo¬ 
sition  while  negative  values  indicate  entrain¬ 
ment. 

With  equations  (10)  and  (21)  in  (26),  wing 
induced  deposition  fluxes  in  the  inner,  Di,  and 
outer  domains,  D^,  become: 


PeWo 

-•■<}] 


(27) 


A  contour  map  of  the  deposition  flux  induced 
by  a  downwashing  wing  at  C/  =  +0.2  is  given 
in  Figure  5b  for  the  prototype  conditions 
described  in  Section  7.0.  We  find  that  the  depo¬ 
sition  is  enhanced  in  the  wake  trail  when  C,  is 
positive.  In  this  case,  the  downwash  behind  the 
wing  has  reduced  or  even  reversed  the  density 
gradient  at  the  lutocline,  (see  Figure  3b),  so 
that  upward  diffusion  from  the  fluid  mud  layer 
is  either  greatly  diminished  or  turned  around  to 
augment  settling.  However,  when  C/  becomes 
negative  (upwashing),  downstream  entrain- 
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Figure  5  (B).  Deposition  flux  for  €  =  4.951  cm^/sec.  0.035  cm/sec.  pi  =  0.0125  g/cm^  p,  =  0.10  g/cm^  and  p,  =  1.2  g/cm^ 
h  =  S  =  304.8,  and  p(0,l)  =  1.  _ 


ment  occurs  in  both  the  inner  and  outer  domain 
as  a  consequence  of  enhanced  upward  diffusive 
fluxes  along  density  gradients  which  have  been 
rendered  more  negative  by  the  presence  of  the 
wing,  see  Figure  3a. 

Over  sufficiently  long  periods  of  time  (many 
tidal  cycles),  the  average  shoaling  flux,  J,  is  the 
net  resultant  of  the  mean  erosion  and  deposi¬ 
tion  fluxes,  or 

J  =  E  -  D  (28) 

where  overbars  denote  time  averages.  Here 
negative  values  of  J  correspond  to  shoaling  or 
losses  in  the  mean  depth  of  water  while  positive 
values  denote  scour  or  depth  gains.  Compari¬ 
sons  of  the  shoaling  performances  between 
upwashing  and  downwashing  wings  at  equal 
and  opposite  C/  are  shown  in  plots  of  shoaling 
flux  from  (28)  along  the  axis  of  the  wake  trail 
in  Figure  6  and  across  the  wake  at  x  =  S/3  in 
Figure  7.  We  find  that  the  upwashing  wing  (C/ 
<  0)  is  predicted  to  be  the  most  efficient  config¬ 
uration  for  scour  or  depth  maintenance.  The 
upwashing  wing  maintains  positive  shoaling 
fluxes  all  along  the  wake  trail  as  a  result  of  ero¬ 
sion  in  the  near  wake  and  entrainment  in  the 
far  wake.  This  action  extends  beyond  the  wing 
tips  in  consequence  of  roll-up  of  the  tip  vortices. 


SHOALING  FLUX  vs  DISTANCE  UP-DOWNSTREAM 
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Figure  6.  Comparison  of  shoaling  flux  theories  and  measure¬ 
ments  along  the  axis  of  the  wake  trail  for  both  downwashing 
0  =  +4^)  and  upwashing  wings  (p  -  4®).  Theory  curves 
based  on  same  parameterizations  as  in  Figure  6. 


The  downwashing  wing,  on  the  other  hand, 
appears  to  be  so  effective  in  advecting  sus¬ 
pended  sediments  toward  the  bottom  that  depo¬ 
sition  overwhelms  erosion  along  much  of  the 
wake  trail. 

The  induced  velocity  field  of  the  wing  disturbs 
future  bottom  profiles  by  fundamentally  differ¬ 
ent  mechanisms  when  causing  erosion  vs  pre- 
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Figure  7.  Comparison  of  shoaling  flux  theories  and  measure¬ 
ments  across  the  wake  trail  at  X  =  S/3.  Theory  curves  based 
on  a  spanwise  circulation  distribution  for  a  tapered  wing 
according  to  (23)  and  parameterizations  as  in  Figure  6. 
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Figure  8.  Variations  of  erosion  and  deposition  fluxes  with 
wing  aspect  ratio  and  mooring  elevation.  (A).  Erosion  flux  at 
X  -  S/4. 


venting  deposition.  In  the  case  of  erosion,  the 
tangential  velocities  of  the  wing  vortex  system 
disturb  the  bottom  directly  by  elevating  the 
shear  stresses  in  the  fluid  mud  layer  above  the 
cohesive  yield  stress  of  the  immobile  bed.  To 
prevent  deposition  the  normal  components  of 
the  induced  velocities  indirectly  disturb  future 
bottom  profiles  through  their  action  on  the  den¬ 
sity  gradients  above  the  lutocline.  These  nor¬ 
mal  induced  velocities  cannot  directly  cause 
advective  fluxes  of  sediment  across  the  lutoc¬ 
line  because  of  vanishing  normal  flow  along  the 
stress  free  boundary  at  z  =  0.  By  either  mech¬ 
anism  shoaling  is  minimized  by  any  given  flow 
Uo  by  maximizing  the  C/  which  the  wing  can 
generate  in  ground  effect.  Thus  the  sensitivity 
of  the  wing’s  lift-drag  polar  to  both  boundary 
proximity  and  low  Reynolds  numbers  is  critical, 
with  little  existing  airfoil  research  as  such  to 
support  optimal  design  selections. 

The  question  of  an  optimal  mooring  elevation 
and  wing  aspect  ratio  is  dependent  upon  the  rel¬ 
ative  importance  of  erosion  vs  deposition  in 
determining  future  bottom  profiles.  Erosion, 
which  prevails  near  the  wing,  is  maximized  by 
mooring  elevations  relatively  close  to  the  luto¬ 
cline,  b  =  0,25S,  and  increases  slightly  with 
increasing  wing  aspect  ratio  as  shown  in  Figure 
8a.  Very  little  improvement  in  maximum  E 
occurs  for  >  5.  This  behavior  is  dominated 
by  the  vortex  sheet  of  the  wing.  Lower  mooring 
elevations  increase  the  near  wake  tangential 
velocities  over  the  bed  and  hence  the  magnitude 


of  the  shear  stress  within  the  fluid  mud  layer. 
Higher  wing  aspect  ratios  diminish  tip  loss 
effects  near  Y  ==  0.  On  the  other  hand,  deposi¬ 
tion  is  minimized  at  mooring  elevations  some¬ 
what  further  above  the  lutocline,  at  b  =  0.8S, 
and  decreases  rapidly  with  decreasing  wing 
aspect  ratio,  as  in  Figure  8b.  Since  the  deposi¬ 
tion  flux  will  be  found  to  account  for  shoaled 
bottom  profiles  in  the  far  wake,  see  Experiment 
section,  such  dependence  on  b  and  is  expli¬ 
cable  in  terms  of  the  vertical  velocities  induced 
by  the  trailing  vortices.  Small  mooring  eleva¬ 
tions  and  large  wing  aspect  ratios  tend  to 
diminish  the  vertical  induced  velocities,  with 
corresponding  increases  in  the  minimum 
attainable  deposition  flux. 


DEPOSITION  FLUX  vs  WING  PARAMETERS 
at  X  «  3S  and  Y  =  0 


Figure  8  (B).  Deposition  flux  at  X  =  3S.  Parameterizations 
as  in  Figure  5,  with  h  5=  6S. 
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EXPERIMENT 

Three  buoyant  wings  were  moored  100  cm 
above  the  initial  bottom  elevation  of  the  Napa 
River  at  a  mean  depth  of  10  meters.  The  config- 
uration  of  the  wings  and  their  moorings  is 
shown  in  Figure  1.  Each  wing  was  constructed 
of  1.6  X  10“  ^  gm/cm^  urethane  foam  reinforced 
by  steel  tendons  and  ribs.  The  semi-span  was 
made  large  relative  to  the  fioc  layer  thickness, 
S  =  3.05  meters,  so  that  large  density  field 
aspect  ratios  are  compatible  with  (5).  A  large 
wing  aspect  ratio  was  required  to  remain 
within  the  purview  of  the  theory  and  to  maxi> 
mize  erosion.  However,  the  need  for  an  ade¬ 
quate  wing  cord  to  satisfy  (5)  and  subsequent 
ordering  arguments  dictated  a  compromise 
selection  for  the  wing  aspect  ratio  of  =  6, 
making  the  mean  semi-cord  c  =  0.61  meters. 

Because  low  wing  Reynolds  numbers  were 
anticipated  in  the  slow  moving,  highly  viscous 
fioc  layer,  a  stall  resistant  model  airplane  air¬ 
foil,  the  HQ-3.5/12,  was  employed  (QUABACK, 
1983).  This  airfoil  also  had  the  advantage  of  a 
nearly  constant  pitching  moment  coefficient 
over  a  wide  range  in  angle  of  attack,  allowing 
the  use  of  a  single  point  mooring  as  diagramed 
in  Figure  1.  The  mooring  point  incorporated  a 
swivel  joint  permitting  complete  rotational 
freedom  to  changes  in  current  direction.  The 
mooring  point  position  and  angle  of  attack  were 
also  fully  adjustable  by  means  of  a  pair  of 
opposed  turnbuckles  under  tension.  The  moor¬ 
ing  attachment  points  were  on  the  suction  side 
of  the  upwashing  wings  with  p  =  -4°,  corre¬ 
sponding  to  C/  =  -0.2  acting  downward 

against  the  buoyancy.  The  downwashing  wings 
were  moored  from  the  pressure  side  with  p  = 

4°,  so  that  C/  -  +  0.2  acting  against  the  moor¬ 
ing  tension.  To  maintain  pitch,  roll  and  yaw 
stability  about  this  single  point  mooring,  the 
wing  plan  form  was  configured  with  30°  of 
sweep  back,  5°  of  dihedral,  and  small  vertical 
stabilizers  on  each  wing  tip. 

Two  upwashing  wings  and  one  downwashing 
wing  were  installed  in  an  80  x  56  meter  test 
area  on  November  11,  1982.  The  test  area  and 
a  pair  of  75  x  50  memter  control  areas  were 
monitored  by  fathometer  soundings,  water  sam¬ 
pling  and  diver  inspection  until  November  6, 
1983.  The  initial  bottom  topography  in  the  test 
area  was  as  shown  in  Figure  9a  prior  to  wing 
installation.  These  contours  were  derived  from 


grid  surveys  using  a  40  kHz  fathometer.  The 
test  area  was  abreast  of  Dry  Dock  #1  at  the 
Mare  Island  Naval  Shipyard  (MINSY)  and  was 
bordered  by  a  sheet  pile  quay  wall  along  the 
shoreward  side.  The  two  upwashing  wings  were 
moored  at  points  A  and  B  in  Figures  9a  and  9b 
while  the  downwashing  wing  was  moored  at 
point  C. 

Two  Marsh-McBirney  electromagnetic  cur¬ 
rent  meters  were  placed  at  the  mooring  eleva¬ 
tion,  b  =  100  cm,  above  the  bottom  at  the  ebb 
and  flood  ends  of  the  test  area.  Seven  months  of 
current  records  indicated  that  the  near  bottom 
currents  were  modulated  by  tidal  action,  con¬ 
sistent  with  assumptions  of  time  variability 
made  in  (5).  Peak  current  amplitudes  were  0(20 
cm/sec)  while  mean  currents  were  u^  -  10  cm/ 
sec. 

To  measure  the  in  situ  values  of  (pt,  pc>  pJj 
bottom  cores  and  water  samples  of  the  lower 
portion  of  the  water  column  were  collected 
throughout  the  experiment  according  to  a 
schedule  regulated  by  barge  and  ship  move¬ 
ments  along  the  river.  Each  sample  volume, 
0(500  ml),  was  determined  to  0(0.5  ml)  accu¬ 
racy  and  Subsequently  suction  filtered  through 
#1  Whatman  filter  paper  with  an  11  micron 
pore  size.  After  ashing  and  drying  in  an  oven  for 
24  hours  to  remove  organics  and  absorbed 
water,  the  mass  of  sediment  retained  by  the  fil¬ 
ter  paper  was  determined  with  a  Mettler  ana¬ 
lytical  balance  to  0(10“®  g)  accuracy.  The  ratio 
of  sediment  mass  to  total  sample  volume  gave 
the  sediment  component  densities  equivalent  to 
equation  (2). 

The  test  period  coincided  with  the  1983  El 
Nino  when  regional  rainfall  and  sediment 
abundance  in  the  Napa  River  were  at  100  year 
record  levels,  (USGS,  1983).  During  the  high 
siltation  period  between  November  1982  and 
July  1983,  the  shoaling  rates  in  two  distinctly 
different  control  areas  were  monitored.  One  of 
these  control  areas  (Pier  21)  was  located  in  a 
quiet  water  cul-de-sac  where  there  were  no 
measurable  bottom  currents  and  hence  no  ero¬ 
sion.  Here  the  mean  shoaling  flux  was  found  to 
be: 


=  -1.185  X  10“'*  gm/cm^  sec  (29) 

where  t\  was  the  vertical  position  of  the  bed 
from  40  kHz  fathometer  soundings.  From  this 
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BOTTOn  TOPOGRAPHY:  BEGINNING  OP  EXPEPinENT  11  NOV  82 


Figure  9  (A).  Comparison  of  initial  and  final  bottom  configurations  during  the  prototype  wing  test.  Upwashing  wings  {Ci  <  0) 
were  moored  at  points  A  &  B.  A  downwashing  wing  (C/  >  0)  was  moored  at  point  C. 


eOTTOri  TOPC'GPAPH-:  END  C'F  EXPERIMENT  o  NOV  S3 


Figure  9  (B).  Comparison  of  initial  and  final  bottom  configurations  during  the  prototype  wing  test.  Upwashing  wings  (C/  <  0) 
were  moored  at  points  A  &  B.  A  downwashing  wing  (C?  >  0)  was  moored  at  point  C. 


observed  shoaling  flux  the  mass  diffusivity  was 
calibrated  by  solving  equation  (27)  for  the  case 
of  no-wings,  Ii  =  I2  =  0?  to  get: 


J(1  -  pc/pa)  -  K^Pcg  +  Pc^o 


=  4.951  cm^/sec 


Here,  pt  =  0.0125  g/cm®,  p,  =  0.100  g/cm^,  and 


p.  -  1.2  g/cm®  derived  from  averages  of  the  in 
situ  water  samples  collected  in  the  control 
areas  throughout  the  duration  of  the  experi¬ 
ment.  The  sedimentation  coefficient  was  taken 
as  K.,  =  4  X  10"  sec  after  the  work  of  FUJITA 
(1962).  A  characteristic  settling  velocity  in 
equation  (30)  of  Wq  =  0.035  cm/sec  was  based 
upon  a  25  micron  median  aggregate  size  deter¬ 
mined  by  scanning  electron  microscopy  size 
analysis  of  the  suspended  load  in  the  Napa 
River  water  samples. 
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The  fluid  mud  parameters  (pi,m,Tc)  were 
derived  from  grab  samples  recovered  by  divers. 
These  samples  were  placed  immediately  in  a 
Brookfield  Viscometer,  giving  average  values  of 
pi  =  3.06  X  10^^  gm/cm  sec  and  =  0.5  dynes/ 
cm^  The  thickness  of  the  fluid  mud  layer  was 
estimated  to  be  0(15  cm),  based  upon  the  band 
of  abrasion  appearing  on  the  upper  portions  of 
the  vertical  steel  members  of  the  shear  plate  as 
shown  in  Figure  1.  With  these  values  and  the 
calibrated  mass  diffusivity  (30),  we  have 
invoked  (27)  and  (28)  to  back  out  the  erosion 
flux  contribution  to  the  observed  mean  shoaling 
flux  at  the  remaining  control  area  (Berth  7), 
and  found  that  a  =  2.356  gm/cm*  sec. 

Sediment  component  densities  measured 
approximately  x  =  S/4  downstream  of  the 
upwashing  wing-B  are  compared  with  the  inner 
solution  calculations  (21)  in  Figure  4a.  The 
majority  of  these  measurements  are  found  to  lie 
in  the  envelope  of  the  theory  curves.  This  enve* 
lope  is  believed  to  be  the  correct  standard  for 
comparison  because  precise  cross  stream  posi¬ 
tioning  during  water  sample  collection  was  dif¬ 
ficult  while  surface  currents  were  running.  In 
general  the  data  show  high  levels  of  suspended 
sediments  raised  above  the  upwashing  wing 
and  locally  smaller  levels  at  the  mooring  ele¬ 
vations,  b  =  100  cm.  Such  a  near  bed  density 
distribution  is  not  typical  of  Nature.  Thus  an 
inversion  layer  has  been  created  in  the  region 
above  the  floe  layer  by  the  upward  advection  of 
the  near  wake  in  reaction  to  downward  lift. 

Comparisons  between  theory  and  measured 
sediment  component  densities  observed  at  x  = 
S/4  downstream  from  downwashing  wing-C  are 
shown  in  Figure  4b.  Samples  yielding  this  data 
were  gathered  on  the  same  days  as  those  in  Fig¬ 
ure  4a  and  show  comparable  values  above  and 
at  the  bottom  of  the  floe  layer.  However,  in  Fig¬ 
ure  4b  a  high  density  toe  of  suspended  sediment 
is  found  below  the  wing  and  above  the  pre¬ 
sumed  lutocline  at  z  =  0.  Some  of  the  sediment 
component  densities  in  this  toe  are  actually  in 
excess  of  those  found  at  the  lutocline.  This,  too, 
is  an  unnatural  finding.  Acoustic  evidence  of 
this  near-bottom  high  density  suspension  is 
found  downstream  (LHS)  of  the  echo  return  for 
wing-C  in  Figure  10a.  The  40  kHz  acoustic 
waves  are  known  to  back  scatter  from  strong 
density  gradients  like  those  appearing  above 
the  bottom  in  Figure  4b,  (KIRBY  and  PARKER, 
1974).  Hence,  the  downstream  acoustic  image 


in  Figure  10a  is  qualitatively  consistent  with 
the  measured  density  anomalies  near  the  bed  in 
Figure  4b.  Together,  these  observations  support 
theoretical  expectations  that  downwashing 
wings  tend  to  drive  suspended  sediment  down 
into  the  lower  portions  of  the  floe  layer,  causing 
an  inversion  layer  below  the  wing. 

Shoaling  flux  estimates  derived  irom  40  kHz 
fathograms  are  compared  with  theory  in  Fig¬ 
ures  6  and  7.  One  of  these  fathograms  is  shown 
in  Figure  10b,  taken  along  the  centerline  of  the 
wake,  (y  =  0),  behind  an  upwashing  wing,  B. 
The  upwashing  wings  (C/  <  0)  were  found  to 
diminish  downstream  shoaling  while  the  down¬ 
washing  wings  (C/  >  0)  were  found  to  cause 
downstream  accretion,  although  the  outer  solu¬ 
tion  for  the  shoaling  flux  tends  to  underesti¬ 
mate  both.  The  momentum  diffusivity  was  fixed 
at  O'  =  36  in  order  to  match  the  length  of  the 
scour  trail  with  the  data.  The  data  were  based 
upon  one  week  averages  shortly  after  the  wing 
installations  were  completed,  when  the  erosion 
fluxes  were  relatively  large. 

In  the  near  wake  of  the  wing,  the  inner  solu¬ 
tion  of  the  shoaling  flux  was  found  in  Figures  6 
and  7  to  overestimate  the  observed  erosion  due 
to  an  upwashing  wing  by  a  considerable 
amount.  This  is  believed  to  be  due  to  the 
assumption  of  a  constant  when  posing  the 
erosion  flux  according  to  (25).  As  erosion  cuts 
into  the  bed,  denser  more  consolidated  mud  is 
exposed  which  has  a  higher  cohesive  yield 
stress  and  is  thus  more  resistant  to  erosion  than 
predicted  by  the  calculations  in  Figures  6  and 
7.  This  argument  is  supported  by  the  fact  that 
the  smaller  positive  shoaling  fluxes  which  do 
not  entail  erosion  cuts  as  deeply  into  the  bed  are 
predicted  reasonably  well  in  the  near  wake  of  a 
downwashing  wing.  Unfortunately,  the  calcu¬ 
lation  does  not  account  for  the  unshoaled  inter¬ 
mediate  region  between  2S  <  x  <  8S  behind  a 
downwashing  wing.  This  may  again  be  due  to 
cross  stream  positioning  errors  if  the  fathome¬ 
ter  sweep  was  performed  along  one  of  the  ero¬ 
sive  tails  of  the  horseshoe  at  Y  =  S  in  Figure 
5a  rather  than  along  the  centerline  of  the  wake, 
Y  =  0. 

The  rate  of  erosion  of  naturally  deposited 
beds  typically  decreases  with  increasing  time 
and  may  eventually  arrest  altogether,  (MEHTA 
et  a/.,  1982).  Not  only  does  erosion  expose  den¬ 
ser,  more  erosion-resistent  material,  it  also 
increases  the  effective  mooring  elevation. 
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thereby  reducing  further  erosion  while  increas¬ 
ing  downstream  deposition  in  accordance  with 
the  results  of  Figures  8a  and  8b.  The  image 
method  calculations  also  become  increasingly 
inaccurate  as  the  bottom  deviates  from  a  plane 


configuration.  Flow  divergence  over  wing 
induced  scour  holes  will  surely  reduce  the  cir¬ 
culation  which  the  wing  can  generate.  These 
factors  all  work  together  over  time  to  slow  or 
arrest  erosion.  A  final  steady  state  bottom  con- 
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figuration  as  shown  in  Figure  9b  was  found  12 
months  later.  The  long-term  erosion  process 
has  produced  local  scour  holes  in  the  neighbor¬ 
hood  of  both  downwashing  and  upwashing 
wings  alike.  Such  features  are  qualitatively 
similar  to  the  erosion  simulation  depicted  in 


Figure  5a.  None  of  the  long  extended  scour 
trails  associated  with  entrainment  and  dimin¬ 
ished  deposition  behind  an  upwashing  wing  are 
still  in  evidence  in  the  final  steady  state  bottom 
configuration.  Thus  deviations  from  a  fiat  bot¬ 
tom  in  consequence  of  near-wake  erosion 
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appears  to  greatly  inhibit  vertical  mass  trans¬ 
port  in  the  wake.  Even  so,  it  is  significant  that 
none  of  the  wings  became  buried  after  one  year 
of  continual  submergence,  in  spite  of  the  higher 
than  normal  siltation  rates  that  persisted 
throughout  the  test  period. 

CONCLUSIONS 

Wing  induced  perturbations  to  the  density 
field  of  suspended  sediment  appear  to  be  qual¬ 
itatively  accounted  at  OCRp/Av^).  Wings  in 
ground  effect  which  generate  a  downward  lift 
force  give  rise  to  density  inversions  above  the 
wing.  Wings  producing  an  upward  lift  force 
cause  density  inversions  to  appear  below  the 
wing. 

Short  term  wing  induced  perturbations  to  the 
bottom  appear  to  be  qualitatively  explained  by 
entrainment  and  diminished  deposition  in  the 
far  wake  and  by  erosion  in  the  near  wake. 
Wings  which  generate  a  downward  lift  force 
induce  erosion  in  both  the  near  and  far  wakes. 
Wings  which  generate  an  upward  lift  force  give 
rise  to  erosion  in  the  near  wake  followed  by 
accretion  in  the  far  wake. 

Long  term  wing  induced  perturbations  to  the 
bottom  are  qualitatively  explicable  only  in 
terms  of  erosion. 

Both  erosion  and  entrainment  are  maximized 
by  maximizing  wing  section  lift  coefficients 
while  in  ground  effect  at  low  wing  Reynolds 
numbers. 

There  exists  no  unique  selection  of  mooring 
elevation  and  wing  aspect  ratio  which  simul¬ 
taneously  maximizes  erosion  and  minimizes 
downstream  deposition.  Erosion  is  maximized 
by  maximizing  wing  aspect  ratio  for  a  mooring 
elevation  of  b  =  0.25S.  Deposition  is  minimized 
by  minimizing  wing  aspect  ratio  with  b  =  0.8S. 

Wings  which  generate  a  downward  lift  force 
are  more  effective  in  causing  local  erosion  and 
retarding  downstream  deposition  than  wings 
generating  an  upward  lift  force.  Because  either 
effect  is  limited  to  a  wake  trail  whose  cross¬ 
stream  dimension  is  of  the  order  of  a  wing  span, 
large  areas  of  channel  bottom  will  require  mass 
arrays  of  wings  with  overlapping  wakes  (like  a 
flock  of  birds),  in  order  to  control  shoaling. 
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n  RESUMEN  D 

Se  muestra  en  este  trabajo  una  teorfa  con  correspondiente  expcrimentacibn  que  describe  la  resuspensibn  de  sedimentos  finos  debida 
a  la  accibn  de  alerones  anclados  cerca  del  fondo  de  un  estuario.  Los  cblculos  mediante  el  mbtodo  de  las  imaj^enes  indican  que  el 
sistema  de  vbrtices  de  un  alerbn  puede  ejercer  dos  efectos  diferentes  sobre  la  sedimentacibn  del  estuario:  1)  Puede  aumentar  la 
tensibn  tangencial  en  las  proximidades  del  alerbn  hasta  una  cantidad  suficiente  como  para  inducir  la  erosibn  de  las  capas  de 
scdimento  de  fondo  parcialmente  consolidadas;  y  2)  puede  promover  transporte  vertical  que  incremente  o  disminuya  la  sedimen¬ 
tacibn  corriente  abajo.  Se  ha  encontrado  que  la  sedimentacibn  en  la  cola  de  la  estela  depende  de  la  direccibn  de  la  fuerza  de 
sustentacibn  generada  por  el  alerbn.  Los  alerones  que  producen  una  fuerza  de  sustentacibn  descendente  limpian  la  regibn  ady- 
acente  del  fondo  de  sedimento  suspendido  y,  por  lo  tanto,  disminuye  la  sedimentacibn  aguas  abajo.  Reclprocamente,  los  alerones 
que  generan  una  fuerza  de  sustentacibn  ascendente  incrementan  la  conccntracibn  de  sedimento  suspendido  cerca  del  fondo.  incre- 
mentando  por  lo  tanto  la  tasa  de  sedimentacibn  aguas  abajo.  Se  formulan  las  configuraciones  de  alerones  bptimas  que  maximizan 
la  erosibn  y  minimizan  la  sedimentacibn.  Las  medidas  de  concentracibn  detrbs  de  prototipos  de  alerones  son  consistentes  con  las 
expectativas  tebricas.  Los  cambios  a  corto  plazo  de  la  batimetna  del  fondo  se  ajuatan  en  general  con  la  erosibn  predicha  en  e! 
dominio  interior  y  con  la  menor  acrecibn  en  el  dominio  exterior.  Department  of  Water  Sciences,  University  of  Cantabria,  Santander, 
Spain. 


□  RESUME  □ 

Prbsente  une  expbrimentation  dberivant  la  resuspension  des  sbdiments  fins  sous  Taction  d*ailes  amarrbes  prbs  du  fond  dans  un 
estuaire.  Les  calculs  rbalisbs  par  ajustement  k  des  expansion  d’images  internes  et  externes  indiquent  que  le  vortex  d’une  aile  peut 
exercer  deux  effets  distincts  sur  la  sbdimentation  estuarienne.  1)  11  peut  accroitre  par  une  quantitb  suffisante  la  force  de  cisaille- 
ment  du  voisinage  et  induire  une  brosion  des  couches  partiellement  consolidbes  du  fond.  2)  11  peut  promouvoir  un  transport  vertical 
qui  soit  augmente,  soit  diminue  le  dbpbt  aval.  Dans  la  tralnbe  du  sillagc,  le  dbpbt  est  dbpendant  de  la  force  d’blbvation  entrainbe 
par  Taile.  Les  ailes  qui  produisent  une  force  d’blbvation  en  aval  provoquent  le  dbpart  des  sbdiments  fins  en  suspension  prbs  du 
fond  et  diroinuent  le  dbpbt  aval.  A  Tinversb,  les  ailes  engendrant  une  force  vers  le  haut  accroissent  la  densitb  des  sbdiments  en 
suspension  prbs  du  fond  et  accroissent  done  le  dbpbt  en  aval.  Les  mesures  de  densitb  derribre  les  ailes  du  prototype  coincident 
avec  ce  qui  pouvait  btre  dbduit  de  la  thborie.  Les  changements  b  court  terme  des  formes  de  fond  coincident  gbnbralement  avec 
Tbrosion  prbdite  dans  le  domaine  interne  et  avec  une  acerbtion  diminube  dans  le  domaine  exteme. — Catherine  Bressolier,  Gio- 
morphologie  EPHE,  Montrouge,  France. 


a  ZUSAMMENFASSUNG  □ 

Bine  Theorie  und  ein  begleitendes  Experiment  werden  vorgestellt.  die  die  Wiederaufnahme  feinkbrniger  Sedimente  als  Suspen- 
sionsfracht  aufgrund  von  Fliigeln  beschreiben,  die  am  Grunde  eines  Astuars  verankert  sind.  Berechnungen  zeigen,  dafi  die  Wir- 
belbildungen  an  einem  solchen  Fliigel  zwei  verschiedene  Auswirkungen  auf  die  Sedimentation  ira  Astuar  haben  kbnnen:  1)  Sie 
kbnnen  die  Schubspannungen  in  der  Umgebung  des  Fliigels  um  ein  solches  AusmaB  erhbhen,  dafl  dadurch  die  Erosion  von  teilweise 
verfestigten  Schichten  am  Astuargrund  ausgelbst  wird;  und  2)  kbnnen  sie  den  vertikalen  Transport  fbrdern,  um  eine  stromab- 
wfirtige  Ablagerung  entweder  zu  steigern  oder  zu  hemmen.  Es  ergab  sich,  daO  die  Sedimentation  im  Kielwasserbereich  von  der 
Richtung  der  Hubkraft  abhftngig  ist,  die  von  dem  Flugel  erzeugt  wird.  Fliigel,  die  eine  nach  unten  gerichtete  Hubkraft  bewirken, 
beseitigen  Schwebgut  in  dem  am  Grunde  angrenzenden  Gebiet  und  vermindern  dadurch  die  stromabwftrtige  Ablagerung.  Umge- 
kehrt  vergrbfiern  Flugel,  die  eine  nach  oben  gerichtete  Hubkraft  bewirken,  die  Schwebgutdichte  nahe  dem  Grunde  und  dadurch 
auch  die  stromabwbrtige  Ablagerungsrate.  Bestmbgliche  Flugelgestaltungen  werden  vorgeschlagen,  die  die  Erosion  maximieren 
und  die  Ablagerung  minimieren,  Dichtemessungen  bei  Prototyp-Fliigeln  stimmen  mit  diesen  theoretisch  abgeleiteten  Erwartun- 
gen  iiberein.  Kurzzeitige  Veranderungen  des  Reliefs  am  Astuarboden  stehen  im  generellen  Einklang  mit  der  vorhergesagten 
Erosion  im  inneren  und  der  verminderten  Anlagerung  im  buBeren  Fliigelbereich, — Helmut  Bruckner,  Geographisches  Institut, 
U niversitat  Diisseldorf,  F.R.G. 
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Transport  of  Fine  Sediments  by 
Hydrostatic  Jets 

Scott  A.  Jenkins,  Saima  Aijaz  and  Joseph  Wasyl 


The  effects  of  hydrostatic  jets  (diose  without  propagating  horizontal  pressure 
gradi^ts)  on  the  transport  and  entrainment  of  fine  grain  estuarial  s^imftnfg 
from  the  near  bottom  fluid  mud  layers  are  considered.  Using  mixing  length 
argum^ts  it  is  proposed  that  the  elevated  shear  stresses  resulting  fiom  these  jets 
transport  su^nded  sediment  vertically  upward  across  the  lutocline,  thereby 
reducing  the  abimdance  of  suspended  sediment  directly  adjacent  to  the 
consolidated  bed.  Through  this  entrainment  the  jets  were  found  to  cause 
reductions  in  the  deposition  rates.  The  extinction  of  jet-induced  bottom  stress 
with  increasing  distance  from  the  jet  was  measured  on  laboratory  scales  and 
found  to  obey  a  power  law  over  several  orders  of  magnitude.  This  power  law 
was  invoked  in  a  vertical  advection-diffusion  model  to  calculate  the  variation  in 
jet  entrainment  flux  with  increasing  distance  away  from  the  jet  discharge  nozzle. 
These  calculations  were  conq>ared  with  field  measurements  conduct^  at  three 
different  sites  for  a  variety  of  jet  discharge  velocities,  diameters  and  flow  rates. 
Whai  low  salinity  surface  water  was  pumped  through  the  jets,  the  entrainment 
flux  at  great  distances  was  larger  than  expected.  Consequently,  the  calculations 
required  the  assumption  of  a  reduced  cohesive  yield  stress  to  matoh  the 
^easured  values.  This  was  supported  by  laboratory  viscometer  me^lirements 
using  native  suspensions  with  artificially  depressed  salinities. 

Introduction 

Water  jets  have  been  enqiloyed  since  the  late  nineteenth  century  as  a  for 
agitation  dredging.  This  application  requires  large  amounts  of  horsepower  in 
order  to  achieve  jet  induced  shear  stresses  in  excess  of  the  critical  erosion  stress 
of  consolidated  muds.  Consequently  the  jets  have  only  a  local  effect  and  require 
some  sort  of  mobile  platform  to  scour  and  erode  large  areas  of  muddy  estuarial 
waterways.  Furthermore  significant  increases  in  turbidity  result  from  this 
practice  which  ultimately  limit  its  application  in  environmentally  sensitive  areas. 

In  recent  years  automated  jetting  systems  have  been  attached  to  piers  and 
quay  walls  to  resuspend  fluid  mud  before  it  remains  immobile  long  enough  to 
form  partially  consolidated  mud,  see  Jenkins,  Bailardnnd  Tnman  (1980),  J^ikins 
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and  Bailard  (1989),  Heinz,  Bailard,  and  Jenkins  (1989).  In  this  application  the 
jetting  system  does  not  function  as  a  dredging  device  which  must  induce  erosion. 
Rather,  the  jetting  system  is  employed  as  a  depth  maintenance  device  which 
must  merely  produce  enough  bottom  stress  to  exceed  die  cohesive  yield  stress 
of  mobile  fluid  mud  in  order  to  effect  its  resuspension.  These  automated  jetting 
systems  fall  into  two  basic  classes;  (1)  pump  based  systems;  and  (2)  ducted 
propeller  based  systems.  The  pump  based  systems  consist  of  a  centrally  located 
pump  which  discharges  through  a  manifold  to  a  line  array  of  fixed  bottom 
mounted  jets,  see  Fig.  1.  Control  valves  at  each  branch  point  along  the  manifold 
regulate  the  entire  pump  discharge  sequentially  through  each  jet  for  a 
predetermined  duty  cycle.  The  bottom  mounted  jets  consist  of  a  convergent 
nozzle.  The  operation  of  the  jet  system  is  controlled  by  a  computer  which 
sequences  its  operations  to  coincide  with  early  ebbing  tide  immediately  after 
slack  water.  When  operated  on  a  daily  basis  only  a  relatively  thin  fluid  mud 
layer  is  reintroduced  into  the  water  column,  rather  than  the  deposition  of  an 
entire  season  as  in  die  case  of  agitation  dredging. 


Figure  1.  Section  view  of  the  pump  based  jet  system  at  the  Port  of  Grays 
Harbor,  Washingtcai. 
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The  propeller  based  syste^  operate  in  a  similar  manner,  but  employ  ducted 
fans  rather  than  centrifugal  impellers  to  drive  the  jet  flow,  see  Fig.  2  These 
systems  are  detailed  in  Bailard  and  Jenkins  (1990).  The  propeller  is  mounted 
^ide  a  vertically  ori^ted  90®  elbow  with  a  horizontal  discharge  at  the  bottom 
ms  configuration  protects  the  propeller  from  bottom  debris  or  slumping  mud 
banks,  as  commonly  found  beneath  piers  and  docks.  The  elbow  may  be  hinged 
to  a^frame  and  articulated  with  a  hydraulic  ram  to  allow  the  jet  to  scan  through 
180°  arcs,  thereby  increasing  the  area  of  influence  of  the  jet  discharge.  Whereas 
the  pump  based  systems  generally  produce  a  high  velocity  low  flow  rate  jet 
through  a  small  convergent  nozzle,  the  propeller  based  system  produces  a 
relatively  low  velocity  jet  at  a  very  high  flow  rate.  We  shall  show  subsequentiy 
that  either  high  velocity  or  high  flow  rate  jets  are  capable  of  elevating  bottom 

^r^ses  sufficiently  to  exceed  the  cohesive  yield  stress  of  fluid  mud  and  thereby 
induce  resu^nsion. 


Kgwe  2.  Schematic  of  the  propeller  based  jet  system  tested  at  the  Port  of  Grays 
Harbor  November  1989  -  February  1990.  Components  are  indicted  as  follows: 
1)  pump  chamber;  2)  400  h.p.  pump;  3)  water  intake  screens;  4)  access  catwalk; 
5)  24  mamfold  pipe;  6)  butterfly  control  valve;  7)  12"  nozzle  supply  hose*  8) 
nozzle  guide  pile;  9)  nozzle  support  assembly;  10)  nozzle  -  13,700  g  p'm 
11)  control  house;  12)  air  tank. 
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Theory 

Consider  a  jet  discharging  horizontally  across  a  partially  consolidated  muddy 
bottom  at  z  =  0,  as  ^own  in  Fig.  3.  Let  tihe  jet  diameter  at  die  discharge  point 
be,  d,  where  the  c^terline  axial  velocity  of  the  jet  discharge  is  Uq.  Above  the 
bottom  at  z  =  0  a  suspension  of  uniform  particles,  each  with  a  setding 

velocity  relative  to  the  fluid,  wq  and  a  solid  density,  Pq  ==  2.65  g/cm^ .  If  the 
density  of  the  fluid  is  pf,  then  the  density  of  the  fluid  sediment  mixture,  p„  may 
be  written: 


Pm  =  +  (1-N)/Jf  (1) 

where  N  is  the  volume  concentration  equal  to  the  volume  of  sediment  per 
volume  of  the  fluid-sedimait  mixture.  The  density  of  the  sediment  conqionent 
of  this  mixture,  p,  is  sometimes  referred  to  as  the  excess  density  and  is  defined: 

P=PqN  (2) 

The  upper  regions  of  the  water  column  are  relatively  deficient  of  suspended 
sediment,  where  the  density  of  the  sediment  component  at  the  free  surface,  z  = 
h,  is  p  =  Pt.  Because  of  setding  under  the  influence  of  gravity  the  abundance 
of  suspend^  sedim^t  bbgins  to  increase  in  die  neighborhood  of  the  bottom 
below  z  =  6,  forming  a  non-cohesive  floe  layer  as  shown  in  Fig.  3,  see  Mehta 
and  Partheniades  (1975)  and  Krone  (1978).  At  the  bottom  of  the  floe  layer,  at 
z  =  m,  there  is  a  rapid  increase  in  bulk  density  referred  to  in  the  recent 
literature  as  a  lutocline,  see  Dyer  (1985)  and  Mehta  (1989).  At  or  below  the 
lutocline  the  suspended  sediment  concentrations  are  sufficiently  high  for  the  fluid 
mixture  to  exhibit  a  non-zero  cohesive  yield  stress,  r^,  and  behave 
therefore  as  a  viscoplastic.  Above  the  lutocline  at  z  =  m,  the  cohesive  yield 
stress  vanishes  and  die  fluid  sediment  mixture  behaves  as  a  fluid,  i.e.  unable  to 
support  shear  stress  at  equilibrium.  The  sediment  component  density  at  z  =  m, 
where  die  fluid  sediment  mixture  first  begins  to  exhibit  a  cohesive  yield  stress, 
is  p  =  Pg.  The  region  between  the  lutocline  and  the  immobile  partially 
consolidated  bottom  at  z  =  0  is  referred  to  as  a  fluid  mud  layer  which  may 
exhibit  a  non-zero  velocity  when  applied  stresses  exceed  the  cohesive  yield 
stress.  Thus  the  lutocline  at  z  =  m,  between  the  fluid  and  the  viscoplastic 
material,  is  not  a  no-slip  boundary.  For  analytic  purposes  we  shall  assume  that 
the  lutocline  at  z  =  m  is  a  stress  free  boundary  much  like  a  free  surface,  with 
the  shear  stresses  resulting  from  fluid  motion  appearing  in  the  stress  field  of  the 
fluid  mud  layer  itself.  This  assunqrtion  relieves  the  equations  of  motion  for  the 
velocity  field  of  the  jet  on  the  fluid  side  of  the  interface  from  having  to  satisfy 
a  no-slip  condition  at  z  =  m.  The  no-slip  boundary  appears  at  the  irmnobile 
partially  consolidated  bottom  at  z  =  0.  We  shall  assume  there  are  minor 
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irregularities  along  the  surface  z  =  0  whose  elevations  are  characteristically,  Zq, 
where  Zq  «  m. 

Below  the  fluid  mud  layer,  z  <  0,  the  sediment  is  at  rest  and  there  is  a 
progressive  buildup  in  sediment  component  density  to  some  ultimate  saturation 
value,  p  =  Pg.  This  is  due  to  consolidation  and  compaction  under  the  weight  of 
die  overburden  which  proceeds  at  a  rate  characterized  by  p^K^g,  where  g  is  the 
acceleration  of  gravity  and  Kg  is  the  sedimentation  coefficient  in  seconds,  Fugita 
(1962).  The  portion  of  the  cohesive  bottom  below  the  fluid  mud  layer  for  which 
Pg  <  p  <  Pg  shall  be  referred  to  as  the  partially  consolidated  mud,  which  remains 
immobile. 

Now  superimpose  a  jet  discharge  across  the  lutocline  surface  at  z  =  m. 
Assume  this  jet  has  operated  a  sufficient  length  of  time  for  all  transient  pressure 
disturbances  to  propagate  into  the  far  field  (the  hydrostatic  assumption). 
Therefore,  along  the  axis  of  the  jet  we  may  assume  a  one-dimensional  steady 
state  in  which  upward  diffusive  mass  fluxes  due  to  the  turbulent  eddy  tran^ort 
of  the  jet  are  balanced  by  downward  advective  mass  fluxes  due  to  the  settling 
of  the  suspended  sediment  under  the  influence  of  gravity,  or; 


Here  €  is  the  eddy  mass  diffusivity  which  may  approximated  from  mixing  length 
theory  as: 


€ 


ako^/r 


(z  +  Zq) 


(4) 


where  Icq  is  Von  Karman’s  constant;  r  is  the  jet  induced  shear  stress;  and  a  is 
the  Peclet  number,  equal  to  a  ratio  of  the  mass  diffusivity  to  the  momentum 
diffusivity. 

Equation  (3)  is  solved  by  successive  integrations  subject  to  the  following 
boundary  conditions:  (1)  the  jet  will  not  entrain  sediment  from  the  fluid  mud 
layer  unless  the  applied  jet  stress,  t,  exceeds  the  cohesive  yield  stress,  t^,  or  p 
=  Pg,  at  z  =  m  for  t  <  t^;  and  (2)  p  =  p^  at  z  =  h  for  t  >  Tg.  A  solution  to 
Eq.  (3)  in  the  domain  between  the  lutocline  and  the  free  surface  subject  to  these 
boundary  conditions  is: 


P  = 


9a^kg(T  -  Tj 


+  PgCXp 


(h-m)  ^Pg 


(5) 


for  T  and  m  <  z  <  h 
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The  first  term  on  the  right  hand  side  of  Eq.  (5)  is  the  disturbance  to  the 
sediment  component  density  profile  induced  along  the  axis  of  the  jet,  while  the 
second  term  represents  the  ambient  density  profile  above  the  lutocline.  The 
sediment  component  density  above  the  lutocline  is  found  to  grow  in  direct 
proportion  to  the  excess  applied  jet  shear  stress  relative  to  the  cohesive  yield 
stress  of  the  fluid  mud  layer,  while  decreasing  with  increasing  particle  or 
increasing  settling  velocity. 

Fine  sediment  deposition  at  the  consolidated  muddy  bottom,  z  =  0  in  Fig.  3, 
is  directly  proportional  to  the  abundance  of  Suspended  sediment  in  the  fluid  iniiH 
layer,  Ariathursu  and  Krone  (1976)  and  Cole  and  Miles  (1983).  This  abundance 
is  determined  by  flie  net  mass  fluxes  across  the  lutocline  at  z  =  m.  The  flux  of 
suspended  sediment  entering  the  fluid  mud  layer  due  to  settling  under  gravity  is 
-pojQ.  These  downward  directed  settling  fluxes  are  the  predominant  mass  fluxes 
which  occur  across  the  lutocline  in  a  dredged  estuary.  However  in  the  presence 
of  a  turbulent  jet,  vertical  transport  due  to  eddy  diffusion  may  become 
important,  especially  along  strong  density  gradients  in  the  neighborhood  of  the 
lutocline.  The  net  of  settling  fluxes  and  diffusive  fluxes  across  the  lutocline 
determines  the  rate  of  change  of  suspended  sedimrat  in  the  fluid  mud  layer. 
However  there  is  only  a  finite  probability  p(0,l),  that  this  suspended  sediment 
will  actually  stick  to  the  bed  at  z  =  0  and  consolidate  to  some  ultimate 
saturation  density,  p^,  see  Krone  (1962).  Therefore  the  deposition  flux,  D,  may 
be  written: 


-KsPcS  +  P(0.1) 


D  = 


dp  I 


(1  -  PjPs) 


(6) 


In  Eq.  (6)  positive  values  of  D  correspond  to  deposition  while  negative  values 
indicate  entrainment  of  suspended  sediment  from  the  fluid  mud  layer  into  the 
floe  layer  and  fluid  above.  In  the  absence  of  jet  induced  shear  stresses,  (t  < 
as  in  control  areas)  the  diffusive  flux  term  in  Eq.  (6)  is  negligable  and  the 
deposition  flux  becomes  a  positive  value  controlled  by  local  settling  fluxes,  p^wq. 
If,  however,  jet  induced  shear  stresses  become  large,  r  >  t^.,  then  a  large 
negative  diffusive  flux  term  arises  in  Eq.  (6)  which  we  refer  to  as  the  jet 
entrainment  flux,  Ej,  or: 


-KsPeS  +  P(0,1) 


-  Re 


9a^ko(T-Tg)^'^(m+Zo) 


,3/2/ 
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log 
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m+ZQ 

2zo 

(7) 


wh^re  Re  is  the  real  part  operator.  Hence  there  is  no  jet  entrainment  flux  when 
^  Tp.  From  Eq.  (7)  it  is  clear  that  the  variation  of  the  jet  entrainment  flux 
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along  the  jet  axis  depends  decisively  on  the  decay  rate  of  the  applied  jet  shear 
stress  with  distance  from  the  jet  discharge,  and  upon  the  cohesive  yield  stress 
of  the  fluid  mud  layer. 

Applied  Jet  Shear  Stress 

A  laboratory  experiment  was  conducted  in  which  a  turbulent  jet  was  discharged 
across  a  plain  immobile  bottom  imbedded  with  a  flush  mounted  shear  stress 
probe.  The  position  of  the  flush  mounted  shear  stress  probe  was  varied  along  the 
axis  of  the  jet.  The  initial  discharge  velocity  of  the  jet,  Uq  was  varied  for 
several  different  jet  diameters,  d.  The  results  were  found  through  recession 
analysis  to  produce  a  systematic  power  law  relation  between  the  applied  jet  shear 
stress  T,  and  the  distance,  x,  along  the  jet  axis  from  the  jet  discharge.  These 
results  are  plotted  in  Fig.  4.  The  regression  line  for  the  data  was  found  to  obey 
the  following: 


T  =  (8) 

“  X 

We  find  that  the  applied  jet  shear  strength  decays  as  the  -2.4  power  of  the 
distance,  x,  from  the  jet  discharge  point.  The  flush  mounted  shear  stress  prote 
indjf-JitpH  that  lines  of  constant  stress  resemble  teardrops  as  diagrammed  in 


Figure  4.  Extinction  rate  of  applied  jet  shear  stress  measured  on  laboratory 
scales. 
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Fig.  3.  The  maximum  width  of  each  teardrop  occurred  at  roughly  2/3  the 
maximum  range  from  the  jet  for  any  given  contour  of  constant  stress.  The  width 
of  the  teardrop  was  found  to  be  approximately  1/3  die  maximum  range  of  the 
shear  stress  contour.  Our  focus  herein  is  along  the  jet  axis  for  which  the 
maximum  range,  X^,  for  any  given  contour  of  constant  applied  jet  shear  stress, 
T,  is  given  by: 


120o„U? 


-0.417 


(9) 


where  Rj  is  the  jet  Reynolds  number  =  U^d/v. 

In  Fig.  5,  plots  of  Eq.  (9)  show  that  the  maximum  axial  range  of  any  given 
stress  contour  increases  either  with  increasing  axial  jet  discharge  velocity  or 
increasing  jet  diameter.  In  otherwords,  X,„  grows  vsdth  increasing  jet  flow  rate, 
Uo'3rd^/4.  Equation  (7)  indicates  that  the  particular  value  of  applied  stress,  r, 
must  exceed  the  cohesive  yield  stress  ,  r^,  in  order  to  effect  deposition  out  to 
Consequ^tly  we  find  in  Fig.  5  that  X^  decreases  for  any  given  selection 
of  U  and  d  with  increasing  strength  of  the  fluid  mud. 


Cohesive  Yield  Stress  of  the  Fluid  Mud  Layer 

Because  the  steady  state  jet  plume  must  induce  motion  of  the  fluid  mud  layer  to 
effect  entrainment,  the  decisive  question  is  whether  the  applied  jet  shear  stress 
exceeds  the  cohesive  yield  stress  in  the  fluid  mud  layer.  To  assess  ttife  strength 
of  the  fluid  mud,  grab  samples  were  gathered  by  diver  in  flie  neigh^rhood  of 
the  consolidated  bottom.  Hydrometer  analysis  determined  that  the  sediments  in 
ffie  fluid  mud  layer  were  25  %  clay,  55  %  silt  and  20%  fine  sand  with  an  organic 
content  of  4  % .  X-ray  diffraction  revealed  that  the  predominant  minerals  present 
in  the  fluid  mud  layer  were  montmorillonite,  illite,  and  kaolinite.  The  cohesive 
yield  stress  was  determined  with  a  Brookfield  synchro-electric  viscometer 
equipped  with  a  special  UL  adapter,  (Brookfield  Engineering  Laboratories  Inc, 
Stoughton,  MA).  The  usefulness  of  the  instrument  in  predicting  this  behavior  of 
cohesive  sediments  both  in  the  field  and  in  the  laboratory  has  been  described  in 
detail  by  Faas  (1990).  The  coaxial  cylinder  geometry  of  the  UL  adapter  provides 
greater  sensitivity  for  low  viscosity  suspensions  and  the  instrument  is  accurate 
to  within  1  %  of  the  full  scale  reading. 

The  grab  samples  from  the  fluid  mud  layer  were  washed  with  deionized  water 
to  remove  sea  salts,  followed  by  screening  to  remove  pieces  of  debris,  shells, 
etc.  The  samples  were  than  sieved  through  a  62  micnon^sieve  to  separate  the  fine 
fis^tion  from  the  coarse  fraction.  Subsequently  the  suspension  was  agitated 
thoroughly  to  break  up  any  large  clumps  or  floes.  The  suspension  was  then 
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Figure  5.  Dependence  of  the  maximum  range  of  applied  jet  shear  stress  upon: 
a)  axial  discharge  velocity  (jet  diameter/  the  Mare  Island  pump  based  system); 
and  b)  jet  diameter  (axial  jet  discharge  velocity/  the  propeller  based  Grays 
Harbor  jet  system).  Symbols  indicate  design  points  for  prototype  tests  shown  in 
Fig.  8. 
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sheared  in  the  Brookfield  viscometer.  The  initial  reading  was  taken  at  zero 
salinity.  Subsequent  readings  were  taken  at  increasing  salinities  of  up  to  30  o/oo. 
Three  sq>arate  runs  were  made  on  the  same  sediment  suspension  and  3  or  4 
viscometer  readings  were  taken  for  each  increment  of  salinity.  There  were  no 
variations  in  the  measured  values  for  all  three  runs  thereby  assuring  the  accuracy 
of  die  viscometer  to  produce  data  of  good  quality. 

Two  sets  of  experim^ts  were  conducted.  One  consisted  of  analyzing  the  yield 
stress  behavior  of  the  natural  sediment  sample  and  the  other  measured  the  yield 
stress  after  removing  die  organic  matter.  The  organic  matter  was  removed  by 
treating  the  sediment  sample  with  30%  hydrogen  peroxide  solution.  The 
sediment  component  density  was  then  calculated  from  the  ratio  of  the  oven-dried 
weight  of  the  sample  to  die  total  volume  of  the  sediment  susp^sion. 

The  yield  stress  was  measured  at  a  fixed  sediment  concentration  while  the 
salinity  was  varied  in  increments  of  5  o/oo.  The  experiments  were  then  repeated 
for  varying  sediment  concentrations.  The  results  show  that  as  the  salinity  is 
increased  from  0  to  5  o/oo,  there  is  a  sharp  increase  in  cohesive  yield  stress 
followed  by  a  gradual  increase  beyond  5  o/oo,  see  Fig.  6.  This  result  is  in 
agreement  with  the  results  of  previous  investigators  in  the  lower  salinity  range, 
see  Edzviald  et  al.  (1974),  Gibbs  (1983).  However,  even  after  5  o/oo,  the  yield 
stress  does  not  level  off  but  continues  to  increase  at  higher  salinities.  This  result 
implies  that  the  pump  intakes  should  be  placed  as  close  to  the  surface  as 
possible,  to  draw  off  the  lens  of  low  salinity  surface  water.  The  resulting 
reduction  in  yield  stress  within  die  low  salinity  jet  plume  would  then,  by 
Eq.  (9),  maximize  the  range  of  deposition  reduction. 

Several  runs  of  experiments  were  conducted  in  the  sediment  concentration 
range  of  10  grams/liter  to  75  grams/liter.  This  range  is  representative  of  the 
seasonal  variations  in  fluid  mud  layers  at  the  prototype  test  sites  at  Mare  Island, 
CA  and  Grays  Harbor,  WA,  see  Jenkins  et  al.  (1980).  The  values  of  yield  stress 
obtained  were  0.85  dynes/cm^  at  10  grams/liter  and  zero  salinity  up  to  5 
djmes/cm*  at  75  grams/liter  and  30  o/oo.  The  yield  stress  in  Fig.  6  for  the 
sediment  concentration  of  75  grams/liter  rises  from  1.5  dynes/cm^  at  zero 
salinity  to  4.26  dynes/cm^  at  3  o/oo  and  finally  to  5.06  dynes  /cm^  at  30  o/oo. 
There  is  a  19%  increase  in  yield  stress  from  3  o/oo  to  30  o/oo.  Similarly,  other 
runs  show  a  percentage  increase  of  15  to  20  %  as  the  salinity  increases  from  3 
to  30  o/oo.  The  sediment  without  the  organic  matter  (Fig.  6)  shows  the  same 
trend  as  the  natural  sediments,  indicating  that  the  organic  matter  is  not  playing 
a  significant  role  in  affecting  the  yield  stress  increase  due  to  the  effects  of 
salinity. 


Jet  Experiments  in  the  Field 

Jet  experiments  were  conducted  at  three  west  coast  sites:  (1)  Mare  Island, 
-Califbmia;  (2)  Terminal  4,  Grays  Harbor,  Washington;  (3)  Terminal  2,  Grays 
Harbor  Washington.  The  Mare  Island  jet  system  was  a  pump  based  system 
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a.  Natural  Sediments 
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b.  Sediments  with  organics  removed 
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Figure  6.  Cohesive  yield  stress  dependence  upon  salinity  and  concentration. 
Right  hand  column  figures  are  based  upon  sediment  samples  with  the  organics 
removed.  Left  hand  column  figures  are  based  upon  natural  sediments  with  the 
organic  content  intact. 

comprised  of  ten  jets  at  914  cm  spacings,  each  with  a  diameter  of  7.3  cm  and 
a  discharge  velocity  of  2,879  cm/sec.  Two  distinctly  different  jetting  systems 
were  tested  at  Grays  Harbor,  Washington.  The  Terminal  4  system  was  pump 
based,  as  shown  in  Fig.  1,  comprised  of  24  jets  at  914  cm  spacings,  each  with 
a  diameter  of  22.2  cm  and  a  discharge  velocity  of  2,207  cm/sec.  The  Terminal 
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2  system  was  a  propeller  based  system  as  diagrammed  in  Fig.  2  conq)rised  of 
a  single  discharge  with  a  91.4  cm  diameter  and  a  307  cm/sec  discharge  velocity. 

Based  upon  cohesive  yield  stress  results  like  those  in  Fig.  6,  the  Mare  Island 
^stem  was  designed  to  produce  an  applied  shear  stress  of  5  dynes/cm*  out  to 
a  maximum  range  of  2,000  cm,  to  allow  for  seasonal  maximums  in  suspended 
sediment.  The  settling  velocity  of  the  floes  in  the  fluid  mud  layer  was 
determined  to  be  0.035  cm/sec  with  typical  excess  densities  in  the  fluid  mud 
layer  of  0.075  g/cm’.  The  saturation  density  of  the  consolidated  at  Mare 
Island  was  determined  to  be  1.1  g/cm^.  Dynamic  shear  stress  analyses  at  Grays 
Harbor  revealed  significantly  stronger  fluid  mud  layers  as  thick  as  150  cm  with 
cohesive  yield  stresses  as  high  as  10  dynes/cm^.  The  floes  of  these  fluid  mud 
layers  were  larger  with  median  settling  velocities  of  0.05  cm/sec.  Typical  excess 
deities  in  the  fluid  mud  layer  at  Grays  Harbor  were  found  to  be  0.1  g/cm^ 
with  saturated  densities  for  the  consolidated  material  on  the  order  of  1.1  g/cn^. 
Accordingly,  the  design  applied  shear  stress  for  both  the  Grays  Harbor  jetting 
systems  was  increased  to  10  dynes/cm*  out  to  a  maximum  range  of  3,000  cm. 
The  design  point  for  the  Mare  Island  system  is  indicated  in  Fig.  5a  while  the 
design  point  for  the  propeller  based  Grays  Harbor  system  is  shown  in  Fig.  5b. 

Fathometer  surveys  were  conducted  over  a  2  month  period  at  Mare  Island  and 
over  2  separate  3  month  periods  at  Grays  Harbor.  Comparisons  between  the 
local  depths  in  the  jetting  areas  to  those  in  control  areas  were  used  to  calculate 
the  jet  entrainmrat  flux  by  the  following: 


^jet  ~  ^control 
At 


-  D 


control 


(10) 


where  17  was  the  vertical  position  of  the  immobile  bed  as  determin^'^from  40 
.JdEIz  fathometer  soundings  in  the  test  areas  and  control  areas,  respectively.  The 
jet  entrainment  flux  is  a  measure  of  the  deposition  flux  which  did  not  occur  due 
to  the  action  of  the  jet  stress  over  a  period  of  time.  At.  Figure  7  shows  a 
bathymetry  contour  plot  from  the  survey  of  the  Terminal  4  pump  based  Grays 
Harbor  system  conducted  between  March  1987  and  June  1987.  These  soundings 
were  conducted  along  2,000  feet  of  waterfront,  one  half  of  which  (along  the 
upstream  portion)  served  as  a  control  area.  Seven  hundred  feet  of  the 
downstream  portion  was  influenced  by  a  700  ft  long  reach  of  24  jets.  It  is  clear 
from  these  bathymetry  surveys  that  the  jets  had  a  pronounced  effect  on 
preventing  deposition.  The  entire  610  meters  of  waterfront  was  dredged  to  -11.6 
meters  MLLW  just  before  the  activation  of  the  jetting  systems  in  March  1987. 
By  June  1987  the  control  area  had  shoaled  to  between  -7  and  -8  MLLW  while 
major  portions  of  the  jet  test  area  inside  30  meters  were  maintained  at  the 
original  dredged  depths.  Surprisingly,  deposition  was  retarded  as  far  out  as  42 
meters  firom  the  jet  discharge  nozzles. 

^  Measured  jet  entrainment  fluxes  from  the  avefage'depths  over  control  areas 
^  jet  test  areas  according  to  Eq.  (10)  are  plotted  in  Fig.  8  and  compared 
agamst  entrainment  flux  calculations  based  on  Eq.  (7)  for  the  respective  jet  and 
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figure  7.  Bathymetry  survey  at  the  end  of  the  test  period  of  the  pump  based 
Grays  Harbor  jet  system^  conducted  between  March  1987  and  June  1987. 
Bathymetry  contours  are  based  upon  echo  soundings  of  a  40  kHz  fathometer. 
Ebb  flow  proceeds  from  left  to  right  with  the  control  area  situated  on  the 
upstream  side  during  jet  operations  in  ebbing  flow  arrows  indicate  the  position 
of  individual  jets.  _ 

sediment  properties  listed.  The  value  for  the  probability  factor,  p(0,l),  was 
adjusted  to  achieve  best  fit  between  the  calculations  and  the  data.  Even  then,  it 
was  necessary  to  base  the  calculations  on  cohesive  yield  stress  values  less  then 
die  design  values  in  order  to  account  for  the  non  zero  entrainment  fluxes 
measured  at  great  distances  from  the  discharge  points.  The  entrainment  flux  data 
suggests  that  the  cohesive  yield  stress  values  during  the  tests  were  effectively 
only  about  1  and  4  dynes/  cm^  at  MINSY  and  Grays  Harbor  respectively. 
However,  these  values  are  consistent  with  the  results  of  Fig.  6,  for  artificially 
depressed  salinities  which  would  arise  in  the  field  from  the  jet  discharges  which 
derive  their  drive  water  from  intakes  near  the  surface. 


Conclusion 

The  results  indicate  that  deposition  of  fine  grain  sediments  may  be  controlled 
through  the  semi-diurnal  application  of  hydrostatic  jets  producing  applied  shear 
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Figure  8.  Measured  deposition  fluxes  derived  from  spatial  averages  of 
fathometer  surveys  as  compared  with  calculated  deposition  fluxes  based  on 
Eq.  (7).  The  dotted  curve  represents  the  theoretical  entrainment  expected  from 
the  pump  based  jet  array  system  at  Mare  Island  Naval  Shipyard  using  the 
following  parameterization:  p  =  0.001,  Zq  =  1  cm,  m  =  10  cm,  and  a  = 
0.002.  The  solid  curve  is  the  theoretical  entrainment  expected  from  the  turbine 
chimney  unit  at  Terminal  2,  of  the  Port  of  Grays  Harbor,  WA,  using  the 
following  values:  p  =  0.0000005,  Zq  =  10  cm,  m  =  150  cm,  and  a"  =  0.02. 
The  dashed-dotted  curve  corresponds  to  the  theoretical  predictions  of  the 
j)eiformance  of  the  pump  based  jet  array  system  at  Terminal  4,  jOje  Port  of 
Grays  Harbor,  WA,  with  parameterization  as  follows:  p  =  0.0002,  Zq  =  10  cm, 
m  =  150  cm,  and  a  =  0.0036. 


stresses  in  excess  of  the  cohesive  yield  stress  of  the  fluid  mud  layer.  This  may 
be  accomplished  by  either  of  two  means:  (1)  application  of  high  velocity  jets  or 
(2)  application  of  large  diameter  jets.  By  either  approach  the  effect  of  controlling 
^position  thru  jet  entrainment  in  the  far  field  is  maximized  by  maximizing  the 
total  jet  flow  rate.  This  suggests  that  the  total  eddy  momentum  content  of  the  jet 
discharp  is  the  limiting  factor  in  facilitating  entrainment  and  hence  controlling 
deposition.  In  other  words  the  turbulent  momentum  of  the  jet  discharge  is 
progressively  diluted  with  increasing  distance  from  the  jet  discharge  point.  At 
some  maximum  range,  the  eddy  momentum  of  the  jet  discharge  is  diluted  to 
such  an  extent  that  it  is  incapable  of  entraining  sediment  from  the  fluid  mud 
layer  across  the  lutocline  interface.  The  unexpected  persistence  of  this 
^ramment  in  the  far  field  of  low  salinity  jets  also  suggests  that  the  jet  salinity 
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may  exert  a  modifying  effect  on  the  cohesive  yield  stress  of  the  fluid  mud, 
i-ftnHftring  it  more  susceptible  to  eddy  transport  across  the  lutocline. 
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[57]  ABSTRACT 

An  apparatus  for  preventing  sedimentation  buildup  in 
waterways  having  a  current,  the  apparatus  being  sub¬ 
mergible  at  a  location  proximate  the  edge  of  the  water¬ 
way  and  adapted  to  produce  a  scour  pattern  of  prede¬ 
termined  configuration.  The  apparatus  includes  an  elon¬ 
gated,  substantially  vertical  tubular  member  having  an 
inlet  assembly  connected  proximate  the  upper  end  of 
the  tubular  member;  a  turbine  fan  unit  mounted  interi¬ 
orly  of  the  inlet  assembly  for  coniroliably  drawing 
water  into  the  apparatus;  and  a  discharge  assembly 
connected  to  the  tubular  member  proximate  its  lower 
end  for  directing  the  flow  of  water  outwardly  from  the 
tubular  member  in  generally  horizontal  direction  at  an 
optimum  flow  rate. 

3  Claims,  4  Drawing  Sheets 
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method  and 

prevention  of  SEDIME^A-nON  IN 

harbors  and  waterways  ^ 

background  of  the  INVaNTlOK 

I  Fidd  of  the  Invention  , 

The  present  invention  relates  generally  to  a  method 

waterways  to  accommodate  aeep  oraft 
vessels.  More  particularly,  the  invention  concerns  a 
method  and  app^tus  for  accomplishmg  daily  resus- 
of  newly  deposited  sediments  using  a  suo- 
mergetl.  turbulent  water  discharge,  thereby  eluniMtmg 
for  periodic  dredging  and  its  associated  ad- 
verse  enviromncatal  impacts. 

2.  Dbcuaaion  of  the  Invention 
Introdocnon  -  Most  present  day  haroore  are 
in  bays  and  estuaries  which  formed  at  the  mouths  of 
mawTriveis.  These  harbors  have  natural  water  depths  - 
of  about  20  feet.  At  this  depth,  ebbmg  water 
are  sufficiently  intense  to  carry  sediments,  mtroduced 
by  the  nvers  and  streams,  through  the  harbor  and  out  to 

the  ,  r 

Although  20  foot  water  depths  were  adequate  tor 

carlv  sailing  vessels,  the  introduction  of  deeper  draft 
vessels  around  the  turn  of  the  century  created  tne  neca 
for  greater  water  depths.  Bucket  and  suction  arwiges 
were  introduced  to  deepen  the  navigation  channels  and 
berthing  areas.  With  these  increased  depths  came  the  J 
need  for  cononual  maintenance  dredging.  This  process 
of  progressive  deepening  and  subsequent  maintenance 
dredging  continues  to  this  day.  For  example,  tray’s 
major  ports  and  harbors  have  water  depths  of  40  feet  or 
more,  reflectmg  the  massive  size  of  modern  deep  draft  3 
tankers  and  bulk  carriers. 

Increasing  the  depth  of  channels  and  berthing  areas 
reduces  the  scouring  action  of  natural  water  currents. 
AS  a  result  these  quiet  water  areas  have  become  natural 
collection  points  for  sediments  which  in  the  pas:  would  • 
be  transported  through  the  area  and  out  to  sea.  Further¬ 
more,  the  deeper  one  dredges  these  areas,  the  faster 
they  fill  ui.  As  a  result,  maintenance  dredging  has  be¬ 
come  an  increasingly  burdensome  probiem. 

After  150  years  of  dredging,  dredged  material  dis-  - 
posaJ  areas  are  becoming  scarce.  The  problem  has  been 
compounded  by  recent  environmental  regulations 
which  severely  limit  available  options  for  dredged  ma¬ 
terial  disposal.  These  regulations  came  about  through 
the  rccogninon  that  harbor  silts  and  clays  have  a  pro- 
cliviry  for  absorbing  heavy  metals  and  other  toxins.  As 
a  resulL  the  cost  of  dredging  has  risen  sharply  in  recent 
years- 

Convcnoonal  dredging  techniques  focus  on  the  effi¬ 
cient  removal  of  consolidated  bottom  sediments.  Con- 
ventsonai  dredges  achieve  their  economies  of  scale  by 
removing  large  volumes  of  sediment  relatively  infre- 
quentiy.  As  a  result,  harbor  channels  and  berthing  areas 
are  typically  dredged  every  year  or  two.  In  between 
oredgings,  sedimentauon  continues  to  occur,  so  that 
towards  the  end  of  a  dredging  cycle,  ship  operators  arc 
forced  to  carry  only  a  panial  load  to  avoid  grounding. 

The  ihoncommgs  associated  with  conveniionai 
dredging  methods  have  led  to  the  search  for  alternative 
methoas  of  maintaining  adequate  water  oeothi.  One 
approach  has  been  to  miemip:  the  sedimentation  pro¬ 
cess,  other  by  preventing  the  influx  of  sediment  to,  or 
the  deposition  of  sediment  within  a  harbor  or  waterway 


area.  This  approach,  termed  sedimentation  control,  has 
proven  most  cffecrivc  in  harbor  berthing  areas.  An 
important  benefit  of  the  sedimentation  control  approach 
is  that  the  costly  dredge  material  disposal  problem  is 
avoided. 

Another  prior  an  sedimentation  control  methoG  in¬ 
volves  the  use  of  the  scour  jet  array.  This  system  typ»- 
callv  consists  of  a  series  of  near-bottom  water  jets  posi¬ 
tioned  alone  the  face  of  a  wharf  or  dock.  The  jet  array 
>  is  powered  by  a  high  pressure  water  pump,  which  dis- 
charees  to  a  common  manifold  pipe.  The  jets  arc  con¬ 
nected  to  the  manifold  pipe  by  a  scries  of  valves  and 

hoses.  n 

Although  scour  jet  array  systems  are  generally  ctiec- 

5  tive  in  preventing  unwarned  sedimentation,  they  are 
Quite  expensive  to  build  and  operate.  Tests  have  shown 
that  the  roost  efficient  method  of  producing  a  water  jet 
with  a  sizable  range  of  scour  is  to  use  a  high  flow  rate, 
low  velocity  discharge.  High  flow  rates,  however,  re- 
°  quire  large  pumps  and  pipes,  thereby  increasing  the  cost 
of  the  scour  jet  array  system. 

An  alternative  method  of  generating  a  high  flow  rate, 
low  velocity  discharge  is  to  use  a  motor  driven  proi«l- 
ler.  A  number  of  agitation  dredging  and  sedimentation 
control  systems  have  been  built  using  open  mounted, 
near-bottom  propellers  to  produce  the  required  scour 
flow.  One  type  of  prior  an  propeller  system  is  described 
in  U.S.  Pat.  No.  3,449,915.  issued  to  T.  E.  Cummings. 

The  prior  an  propeller  type  systems  typically  suffer 
from  several  fundamental  problems:  (1)  the  open  pro¬ 
peller  design  is  subject  to  fouling  or  damage  by  water¬ 
borne  debris:  (2)  the  unprotected  propeller  can  be  bur¬ 
ied  bv  sediment  avalanching  from  adjacent  undredged 
„  banks;  (3)  the  propeUer  type  systems  utilize  high  salinity 
"  (and.  therefore.  less  efficient)  bonom  water  m  produc- 
inc  a  turbulent  scour  plume;  and.  (4)  the  exposed  pro¬ 
peller  poses  a  hazard  to  divers  and  manne  mammals. 
The  method  and  apparatus  of  the  present  invention 
^  overcomes  the  drawbacks  of  prior  an  systems  through 
the  use  of  an  array  of  turbine  chimney  scour  units  of 
hichlv  novel  desien.  Compared  to  conventional  prior 
art  dredging  methods,  the  turbine  chimney  scour  sys¬ 
tem  of  the  present  invention  is  less  costly  and  considera- 
4^  bly  less  damaging  to  environmental  wa.ter  quality.  By 
continually  preventing  sediment  deposition,  the  turbine 
chimney  scour  svstem  avoids  the  costly  problem  of 
dredged  material' disposal.  AddttionaUy,  reducing  the 
residence  time  of  the  resuspended  sediments  minimizes 
50  the  very  serious  environmental  problem  of  adsorption 
and  release  of  heavy  metals  and  other  toxins.  The  tur¬ 
bine  chimney  scour  system  also  minimizes  the  impact  to 
benthic  organisms  which  are  periodically  destroyed  by 
conventional  dredging  methods. 

<5  The  turbine  chimney  scour  system  of  the  present 
invention  is  also  considerably  more  effective  and  less 
environmentallv  damaging  than  the  prior  an  agitation 
dredging  methods.  The  prior  an  agitation  dredging 
methods  typically  function  by  remobihzmg  consoli- 
60  dated  bottom  sediments,  allowing  ambient  currents  to 
carrv  them  from  the  area.  For  example,  one  common 
prior  an  method  utilizes  a  specially-designed  vessel  to 
iuck  sediment  off  the  bottom  and  discharge  it  honzon- 
tally  across  the  water  surface.  Tne  objective  is  to  raise 
:  65  me' materia!  high  enough  within  the  water  column  so 
that  it  is  carried  a  significant  distance  downstream  bj 
ebbing  currents.  .Agitation  dredging  is  typically  per- 
forme'd  only  once  or  twice  a  year.  The  lengthy  penod 
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of  nme  between  dredging  is  critical,  for  it  allows  the 
sediments  to  consolidate,  greatly  increasing  their  shear 
strength  and  resistance  to  resuspension.  As  a  result,  ^ 
agitation  dredging  methods  are  not  particularly  effec-  «« 
tive,  since  they  fail  to  fully  resuspend  the  consolidat^  ’  me 
bottom  sediments,  causing  them  to  deposit  a  shon  db* 

tance  downcurrent.  /■u  rn* 

In  contrast  to  agitation  dredging,  the  apparatus  of  the 
present  invention  b  typically  activated  on  a  daily  basis.  “ 
fully  resuspending  the  thin  layer  of  newly  deposited  ^ 
sediment.  Thb  provides  three  important  benefits:  (1)  . 

water  turbidity  is  less  because  less  sediment  is  resus¬ 
pended  during  each  operation;  (2)  transport  distances  ^ 
arc  greater  because  the  sediment  is  fully  resuspended;  ^5 
and,  (3)  fewer  toxins  arc  released  because  the  bottom 
sediments  have  a  negligible  residence  time  before  being  ^ 
resuspended- 

Compared  with  the  pump  based  prior  art  scour  jet  ^ 
array  systems,  the  turbine  chimney  scour  system  of  the  20  ^ 
present  invention  is  less  costly  and  considerably  more 
cfTicicnL  Cost  savings  accrue  primarily  through  the  use 
of  individual  turbine  fan  units  to  produce  the  required  ^ 

jet  discharge.  This  eliminates  the  need  for  a  long  and  ci 

costly  laree  diameter  piping  system.  Moreover,  by  uti-  ^ 
lizing  a  high  flow  rate/low  velocity  discharge,  the  tur¬ 
bine  chimney  scour  system  requires  about  1/lOth  the  p 
input  power  to  achieve  the  same  scour  performance.  ii 

The  turbine  chimney  scour  system  of  the  present  v 

invention  is  readily  distinguishable  from  and  clearly 
superior  to  the  prior  an  open  propeller  scouring  system  c 
for  several  reasons.  In  the  first  place,  the  turbine  fan  is  t 
a  different  and  intrinsically  more  efficient  propulsion 
method  than  an  exposed  propeller.  The  ducting  around  35 
the  turbine  fan  acts  to  produce  a  hi^er  discharge  ve¬ 
locity  per  unit  j>ower  input.  Hnclosing  the  propulsion  ^ 
device,  also,  protects  divers  and  unsuspecting  marine 
TT^ammals  from  accidental  injury. 

Secondly,  the  abUity  to  readily  adji^t  the  intake  loca-  40 
don  of  the  apparatus  allows  low  salinity,  oxygen  nch 
surface  water  to  be  used  for  scouring.  Tests  have  shown 
that  low  salinity  water  reduces  the  cohesive  strength  of 
the  bonom  sediments  and  b  more  effective  in  res^ 
pending  fine  bonom  sediments.  Additionally,  the  higher 
oxygen  content  of  the  ingested  surface  water  provides 
improved  water  quality  near  the  bottom,  which  tends  to 
be  anaerobic  with  consequently  high  bacterial  levels. 

The  combination  chimney/elbow  configuration  of  the 
apparatus,  also,  provides  excellent  protection  of  the 
turbine  fan  unit  from  sediment  avalanching  of  ^jacent 
mud  banks.  Field  tests  of  propeller  based  scour  jet  array 
systems  have  shown  thb  to  be  a  serious  problem. 

Ajiother  advantage  of  the  apparatus  of  the  invention  jj 
resides  in  the  novel  design  of  the  screened  intake  wmch 
helps  to  prevent  damage  to  the  propubion  unit  from 
waterborne  debris.  The  design,  abo.  prevents  the  mges- 
tion  of  fish  during  operation.  The  latter  b  an  important 
permitting  consideration  in  areas  having  sigmfiCMt  go 
marine  fisheries  (e.g.  the  salmon  fishery  in  the  Pacific 
Northwest).  As  an  added  feature,  the  flow  from  the 
turbine  fan  can  be  briefly  reversed  following  operation, 
helping  to  clean  the  screen  of  any  accumulated  deons. 
Still  another  advantage  of  the  apparatus  of  the  invention  65 
b  the  gravity  support  frame  which  permits  rapid  instm- 
lation  and  recovery  of  the  turbine  chimney  scour  umts 
for  routine  maintenance. 


SUMMARY  OF  THE  INVENTION 
It  b  an  object  of  the  present  invention  to  provide  :m 
efficient,  cost  effective  and  environme^ly 
method  and  apparatus  for  preventing  sediment  build-up 

in  harbor  and  waterway  areas. 

Another  object  of  the  invennon  is  to  provide  a 
method  and  apparatus  for  reducing  adverse  enviro^ 
mental  impacts  associated  with  convenuonal  pnor  an 
methods,  including  release  of  sediment  borne  toxins  into 
the  water  column  release  of  excessive  sedimrat  turbid¬ 
ity  and  the  periodic  killing  of  benthic  orgamsms.,  .^ 
Another  object  of  the  invennon  is  to  pro^de  a 
method  and  apparatus  for  preventing  sedunent  deposi¬ 
tion  which  uniquely  uses  a  safe  and  efficient  honzont^ 
directed,  turbulent,  ne^bonom  water  discharge  typi¬ 
cally  operated  on  a  daily  basb. 

Still  another  object  of  the  invention  is  to  provide  a 
safe  and  easy  to  use  method  and  apparatus  for  prevrat- 
ing  sediment  deposition  in  harbors  and  waterways  with¬ 
out  interfering  with  ship  navigation. 

Another  object  of  the  present  invention  is  to  provide 
a  method  and  apparatus  for  significantly  reducing  the 
cost  of  maintaining  adequate  water  depths  m  harbors 
'  and  waterways. 

Still  another  object  of  the  present  invention  is  to 
provide  a  method  and  apparatus  which  effectively  clim¬ 
ates  the  need  for  periodically  dredging  harbor  and 
waterway  areas  and  disposing  of  the  dredged  material. 

^  A  further  object  of  the  invention  is  to  provide  a 
method  and  apparatus  in  which  low  salinity  water  can 
be  used  for  scouring. 

BRIEF  description  OF  THE  DRAWINGS 
^  FIG.  1  b  a  front  view  of  one  form  of  the  apparatus  of 
the  invention. 

FIG.  2  is  a  side  clevarional  view  of  the  app^atus. 
fig.  3  is  a  cross-sectional  view  taken  along  lines  3—3 

r\  of  FIG.  1.  r  e 

FIG.  Aba  front  view  of  an  alternate  form  of  appara¬ 
tus  of  the  present  invention.  . 

FIG.  5  is  a  side  view  of  the  apparatus  shown  in  riu. 

4 

,,  ’  FIG.  6  b  a  generally  diagrammatic  view  illustrating 

alternate  scour  patterns  obtainable  through  the  use  of 
the  apparatus  of  the  present  invention. 

FIG.  7  is  a  front  view  of  another  alternate  form  of  the 
apparatus  of  the  invention  comprising  a  gravity  support 
50  f^e  which  rests  against,  but  b  not  attached  to  the  pier 

^  PIG.  g  is  a  side  view  of  the  apparatus  shown  in  FIG. 
7. 

FIG.  9  is  a  cross-sectional  view  taken  along  lines  9  9 

PIG.  10  is  an  enlarged  schematic  view  of  the  turbine 
assembly  of  the  apparatus. 

description  of  the  INVENTION 

60  Referring  to  the  drawings  and  particularly  to  FIGS. 
1,  2  and  3,  one  form  of  the  apparatus  of  the  present 
invention  is  there  illustrated.  Here,  the  apparatus 
shown  being  used  for  preventing  sedimentation  buud- 
up  in  a  harbor  berthing  area  of  the  character  having  a 
65  wharf  W  fronting  a  berthing  area  B.'  In  this  embodiment 
of  the  invention,  the  apparatus  compnses  a  plurality  ot 
scour  units  10  and  connected  to  the  wharf  W  m  a  side- 
by-side  relationship  (FIG.  6). 
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jcoor  unit  10  compraes  i  vertical  tubular  mem¬ 
ber  12  Itf  wing  a  vertical  azn»  a  **ppcr  cod  14  and  a  lower 
cad  14^  C4:eiocctrd  proximaie  the  upper  end  of  each 
tubular  member  12  ia  an  inlet  means  for  directing  the 
flow  of  wate*  from  a  location  near  the  surface  of  the  5 
water  of  the  berthing  area  into  mbuiar  member  IZ  This 
inlet  b  provided  in  the  form  of  a  bcD-mouth 

inlet  aa«ibly  18.  As  best  seen  in  FIGS.  1  and 
Z  amernWy  18  is  connected  with  tubular  member  12  by 
of  mating  flanges  20  and  22  provided  respec- 
nvdy  on  members  18  and  IZ  Flanges  20  and  22  are 
iniereocmectod  together  by  suitable  fasteners,  such  as 
bolts  24.  With  this  arrangement  inlet  members  of  van- 
ow  heights  can  be  readDy  interconnected  with  the  tubu¬ 
lar  12  so  IS  to  pOBtioo  the  upper  end  of  the 

beH-shaped  inlet  assembly  18  at  an  optimum  location 
below  the  surface  of  the  water  in  the  particuiar  berthing 
area  bemg  scoured. 

Mounted  imerioriy  of  each  inlet  member  18  is  means 
for  cootroUably  drawing  water  from  the  near  surface  of  * 
the  berthing  area  into  the  inlet  means.  As  best  seen  in 
RG.  IC,  in  the  present  embodiment  of  the  invention, 
this  TTV^"^  is  provided  as  a  turbine  assembly  26,  which 
includes  impeller  blades  28.  stator  blades  28a,  a  nose 
fainng  28i.  and  a  tail  fairing  28c.  Stator  blades  28c  are 
attached  to  the  inlet  member  or  chimney  (RG.  10).  The 
impeUer  blades,  or  fan.  are  controUably  rotated  by  a 
motor  means*  such  as  a  submersible  electric  motor  30 
mounted  inside  motor  housing  30a.  Electric  or  hydrau- 
lie  motors,  suitable  for  use  in  the  present  application,  are 
wcil-kiK>wn  to  those  skilled  in  the  art  and  are  readily 
commercially  available.  However  the  fan  stator  blades 
and  fainng  are  custom  designed  for  the  defined  use. 

Coimected  proxiniate  the  lower  end  of  each  tubular 
member  12  is  a  discharge  means,  shown  here  as  an 
elbow  shaped  member  48.  for  directing  the  flow  of 
water  outwardly  therefrom  in  a  generally  horizontal 
direcnoD  lo  produce  a  near  bottom  scour  pattern  of  a 
prtrieTer»^^<*d  configurarioiL  As  best  seen  in  FIG.  6.  4^ 
the  scour  units  10  are  strategically  positioned  along  the 
wharf  W  so  that  the  scour  patterns  P  produced  by  the 
discharge  means  overlap  one  another  in  an  optimized 
maziDer. 

R.eferTing  particularly  to  RG.  3,  actuating  means,  45 
gcncraily  designated  by  the  numeral  3Z  are  operably 
conDccted  to  each  tubular  member  12  to  affect  a  desired 
degree  of  rotatxjn  of  tubular  member  as,  for  example, 
sciecQveJy  through  angles  of  90  or  180  degrees.  As  best 
by  referring  to  RGS.  2  and  3,  each  of  the  tubular  50 
12  is  rotatably  mounted  within  vertically 
spaced  apart  collar  assemblies  34a  and  34^.  These  collar 
assemblies  are.  in  turn,  fixedly  interconnected  to  adja¬ 
cent  wharf  pilings  36  by  brackets  37a  and  31b  in  the 
manner  shown  in  RG.  1.  55 


In  the  present  embodiment  of  the  invention,  the  hy¬ 
draulically  actuated  means  comprises  a  hydraulic  cylin¬ 
der  42  within  which  a  connector  arm  and  piston  assem¬ 
bly  44  is  carried.  Hydraulic  cylinder  42  is  pivotally 
interconnected  with  mounting  bracket  31b  and  co^ec- 
tor  arm  44  is  pivotally  inlcrconnccied  with  car  40  in  the 
manner  shown  in  FIG.  3.  Reciprocal  movement  of  the 
piston  interiorly  of  cylinder  42  urges  telescopic  move¬ 
ment  of  connector  arm  44  inwardly  and  outwardly  of 
the  cylinder  42,  which,  in  turn  imparts  rotational  move¬ 
ment' to  ring  39  and  to  tubular  member  12  to  which  it  is 
connected.  With  this  construction,  as  tubular  member 
12  is  rotated,  the  discharge  means,  shown  here  as  an 
elbow  48.  which  is  interconnected  to  tubular  member  12 
proximate  its  lower  end  16,  is  also  rotated  through  a 
predetermined  angle.  As  the  tubular  member  12  and 
discharge  elbow  48  are  rotated  by  the  actuating  means, 
a  scour  pattern  of  the  character  illustrated  in  the  lower 
part  of  FIG.  6,  is  achieved.  In  the  drawings,  this  scour 
pattern  is  shown  as  semicircular  in  shape.  However,  the 
scour  pattern  may  be  of  any  angle  desired  for  the  opti¬ 
mum  displacement  of  sedimentation  within  any  particu¬ 
lar  berthing  area  or  waterway. 

An  important  aspect  of  the  apparatus  of  the  present 
invention,  is  means  for  preventing  ingestion  into  the 
scour  units  of  fish,  debris  and  other  articles  larger  than 
a  predetermined  size.  In  the  embodiment  of  the  inven¬ 
tion  shown  in  FIGS.  1  and  2.  this  means  is  provided  in 
the  form  of  a  mesh  screen  50  which  covers  the  open  end 
of  the  bell-mouthed  shaped  inlet  18. 

Referring  particularly  to  the  upper  portion  of  FIG.  6, 
the  scour  units  10  are  arranged  along  the  face  of  the 
wharf  in  the  manner  shown  to  provide  near  uniform 
scour  coverage  throughout  the  protected  area.  Spacing 
between  the  assemblies  is  governed  by  the  lateral  di¬ 
mensions  of  the  individual  scour  patterns.  In  the  ibrm  of 
the  apparatus  shown  in  the  drawings,  the  discharge 
means  are  strategically  configured  to  produce  fixed 
scour  patterns  of  a  generally  tear-dropped  conngura- 
tion  having  a  length  to  width  ratio  of  about  3.  To  mini¬ 
mize  zones  of  incomplete  scour,  the  assemblies  arc  ar¬ 
ranged  so  as  to  produce  approximately  20  percent  over¬ 
lap  ^iween  adjacent  scour  patterns.  As  indicated  in  the 
lower  portion  of  FIG.  6,  rotating  the  discharge  means  in 
the  manner  previously  discussed,  produces  arcuate 
shaped  scour  patterns  which  span  up  to  180  degree 
sectors  depending  on  whether  the  jets  are  operated  only 
during  ebb  and  flow  conditions  or  during  both  ebb  and 
flood  flow  conditions. 

Turning  now  to  FIGS.  4  and  5,  another  form  of  the 
apparatus  of  the  present  invention  is  there  illustrated. 
The  apparatus  of  this  form  of  the  invention  is  similar  in 
many  respects  to  the  earlier  described  apparatus  and 
like  numbers  arc  used  to  identify  like  components.  Each 
scour  unit  54  of  this  embodiment  of  the  invention  com- 


Collar  aaaeanblies  34c  and  34^  are  of  identical  con-  prises  a  vertical  tubular  member  56  having  a  vertical 

smicoon  each  including  angularly  inwardly  extending  axis,  an  upper  end  58  and  a  lower  end  60.  Connected 

arms  35c  and  35^  which  are  connected  to  brackets  37a  proximate  the  upper  end  of  each  tubular  member  56  is 

and  31b.  Collar  assemblies  34a  and  34^  also  include  an  inlet  means  for  directing  the  flow  of  water  from  a 

transverse  members  38  which  engage  the  adjacent  pil-  60  location  near  the  surface  of  the  water  of  the  berthing 


in p  36  to  counteract  the  forces  generated  by  the  dis¬ 
charge  of  water  through  elbows  48.  A  ring  39  having  an 
car  40  a  connected  to  each  tubular  member  12  interme¬ 
diate  collars  34.  E>isposed  intermediate  ear  40  and  a 
mounung  bracket  31b  which  interconnects  the  lower 
collar  assembly  34^  with  the  piling  36  is  a  hydraulically 
actuated  means  for  urging  rotational  movement  of  ring 
39  and  m  turn  tubular  member  12, 


area  into  tubular  member  56.  This  inlet  means  is  pro¬ 
vided  in  the  form  of  a  generally  “T”shaped,  inlet  assem¬ 
bly  62  integrally  formed  with  tubular  member  56. 

Mounted  interiorly  of  each  tubular  member  56  is 
ti  means  for  conirollably  drawing  water  from  the  near 
surface  of  the  berthing  area  into  the  inlet  means.  In  this 
second  embodiment  of  the  invention,  this  means  is  pro¬ 
vided  as  a  turbine  assembly  64,  which  includes  a  fan  66, 
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a  stator  66a  and  a  fairing  66b,  The  fan  is  controUably 
rotated  by  a  pneumatic  motor  means,  such  as  a  submers- 
ibie  hydraulic  motor  68.  Hydraulic  motor  68  is  con¬ 
nected  to  a  central,  deck-mounted  hydraulic  power  unit 
70,  via  a  pair  of  supply  hoses  72  and  74  (FIG.  5).  Power  5 
unit  70  feeds  a  supply  manifold  loop  (not  shown)  with 
branch  lines  to  each  scour  unit.  Flow  to  each  branch 
line  is  controlled  by  a  control  device  76  operably  assc^i- 
ated  with  power  unit  70.  Power  units  and  manifolding 
and  control  devices  of  the  character  described  ^e  well-  10 
known  to  those  skilled  in  the  art  and  are  readily  com¬ 
mercially  available. 

Routably  connected  proximate  the  lower  end  of  each 
tubular  member  56  is  a  discharge  means,  shown  here  as 
an  elbow  shaped  member  57,  for  directing  the  flow  of  15 
water  outwardly  therefrom  in  a  generally  horizontal 
direction  to  produce  a  near  bonom  scour  pattern  of  a 
predetermined  configuration.  As  was  the  case  with  the 
earlier  described  embodiment,  the  scour  units  54  are 
strategically  positioned  along  the  wharf  W  so  that  the  20 
scour  panems  P  produced  by  the  discharge  means  over¬ 
lap  one  another  in  an  optimized  manner.  In  this  ^ond 
form  of  the  invention,  the  tubular  members  56  are 
fixedly  mounted  and  only  the  discharge  means  is  rotat¬ 
able  to  produce  the  scour  patterns  shown  in  the  lower  25 
portion  of  FIG.  6. 

Referring  particularly  to  FIG.  5,  actuating  means, 
generally  designated  by  the  numeral  JO.  are  operably 
connected  to  each  discharge  elbow  57  to  affect  a  de-' 
sired  deirrec  of  rotation  of  the  elbow,  for  example,  sciec-  30 
tively  though  angles  of  90  or  180  degrees.  The  actual 
ing  means  80  is  of  similar  construction  to  that  previ¬ 
ously  described  and  shown  in  FIG.  3.  More  pameu- 
lariy,  means  80  comprises  a  hydraulic  cylinder  42 
within  which  a  connector  arm  and  piston  assembly  44  is  35 
carried.  HydrauHc  cylinder  42  is  pivotally  intercon¬ 
nected  with  mounting  bracket  82  (FIG.  5)  and  connec¬ 
tor  arm  44  is  pivotally  interconnected  with  a  ring  84 
having  car  86.  In  the  present  instance,  however,  ring  84 
is  connected  to  discharge  elbow  57  which,  in  turn,  is  40 
rotatably  mounted  to  tubular  member  56  proximate  its 
lower  end  60.  Reciprocal  movement  of  the  piston  interi¬ 
orly  of  cylinder  42  urges  telescopic  movement  of  con¬ 
nector  arm  44  inwardly  and  outwardly  of  the  cylinder 
42.  which,  in  turn  imparts  rotational  movement  to  ring  45 
84  and  to  elbow  57  to  which  it  is  connected.  With  this 
construction,  as  elbow  57  is  rotated  by  the  acn^tion 
means,  a  scour  pattern  of  the  character  illustrated  in  the 
lower  pan  of  FIG.  6,  is  achieved.  As  before,  while  the 
scour  partem  is  shown  as  semicircular  in  shape,  the  50 
scour  pattern  may  be  of  any  angle  desired  for  the  opti¬ 
mum  displacement  of  seduncniation  within  any  particu¬ 
lar  berthing  area  or  waterway. 

The  apparatus  of  the  second  form  of  the  invention, 
also,  includes  means  for  preventing  ingestion  into  the  55 
scour  units  of  fish,  debris  and  other  articles  larger  than 
a  predetermined  size.  This  means  is  here  provided  in  the 
form  of  a  cylindrically  shaped,  closed  end  mesh  screen 
90  which  forms  a  part  of  the  “T**  shaped  inlet  62.  The 
“T’  shaped  screen  provides  an  increased  surface  area  60 
which  reduces  screen  intake  velocities  sometimes. 

Referring  particularly  to  the  upper  portion  of  FIG.  6, 
the  scour  units  54  are  arranged  along  the  face  of  the 
wharf  in  the  manner  previously  described  to  provide 
near  uniform  scour  coverage  of  the  character  shown  65 
throughout  the  protected  area.  ■  r-tr-c 

In  the  embodiment  of  the  invention  shown  in 
1  through  3,  as  well  as  that  shown  in  HGS.  4  and  5,  the 


support  frame  is  designed  to  span  two  adjacent  pilm^ 

36.  Chains  G  (91)  which  are  connected  to  the  wharf, 
help  to  support  the  weight  of  the  scour  units  and  to  hold 
them  in  place  relative  to  the  pilings.  The  frame  allows 
the  thrust  from  the  jet  discharge  to  be  carried  back 
against  the  face  of  the  pilings.  The  use  of  simple  chain 
hoist  and  flexible  supply  hoses  72  and  74  allows  m 
entire  scour  unit  to  be  raised  out  of  the  water  for  servic¬ 
ing.  Hoist  arrangements  of  this  character  are  well 
known  to  those  skilled  in  the  an  and  the  details  thereof 
will  not  be  described  herein. 

Referring  to  RGS.  7  through  9,  still  another  form  of 
the  apparatus  of  the  invention  is  shown.  In  this  form  of 
the  invention,  an  alternative  suppon  method  is  illus¬ 
trated.  .  ,  . 

The  apparatus  of  this  third  form  of  the  mvennon  is 
similar  in  many  respects  to  the  earlier  described  appara¬ 
tus  and  like  numbers  arc  used  to  identify  like  compo¬ 
nents.  Each  scour  unit  80  of  this  embodiment  of  the 
invention  comprises  a  vertical  tubular  member  93  hav- 
ing  a  vertical  axis,  an  upper  end  94  and  a  lower  end  95. 
Connected  proximate  the  .  upper  end  of  each  tubular 
member  93  is  an  inlet  means  for  directing  the  flow  of 
water  from  a  location  near  the  surface  of  the  water  o 
the  berthing  area  into  tubular  member  93.  This  inlet 
means  is  provided  in  the  form  of  an  inlet  assembly  mie- 
erally  formed  with  tubular  member  93  and  having  a 
convex  shaped  inlet  screen  90  covering  the  inlet  port  92. 

Mounted  interiorly  of  each  tubular  member  93  is 
means  for  controUably  drawing  water  from  the  nw 
surface  of  the  berthing  area  into  the  mlct  means  (FIG. 
8).  In  this  embodiment  of  the  invention,  this  means  is 
provided  as  a  turbine  assembly  26,  of  the  character 

shown  in  HG.  1.  , 

Connected  proximate  the  lower  end  of  each  tubular 
member  93  is  a  discharge  means,  shown  here  as  an 
elbow  shaped  member  97,  for  directing  the  flow  of 
water  outwardly  therefrom  in  a  generally  horizontal 
direction  to  produce  a  near  bottom  scour  pattern  of  a 
predetermined  configuration.  As  was  the  case  with  the 
^lier  described  embodiment,  the  scour  units  80  are 
strategically  positioned  along  the  wharf  W  so  that  the 
scour  patterns  P  produced  by  the  discharge  means  over¬ 
lap  one  another  in  an  optimized  manner. 

An  important  aspect  of  the  form  of  the  invention 
shown  in  FIGS.  7  through  9,  is  the  design  of  the  suppon 
frame  96  which  comprises  upp>er  and  lower  assemblies 
96c  and  96b.  Importantly,  frame  96,  which  is  connected 
to  and  supports  tubular  mem^r  93,  rests  against,  but  is 
not  connected  to  the  pier  pUings  36. 

Upper  frame  assembly  96a  comprises  a  generally 
rectangularly  shaped  face  member  98  which  is  con¬ 
nected  to  upper  and  lower  braces  100  (FIG.  8),  which 
arc,  in  turn,  connected  to  support  rings  102  mounted  on 
member  93.  Lower  frame  assembly  96^  also  comprises  a 
generally  rectangularly  sha^  face  member  104  having 
a  cross  member  105  which  is  connected  to  a  brace  106, 
which  is,  in  turn,  connected  to  a  support  ring  108 
mounted  on  member  93  {FIG.  8). 

Turning  to  FIG.  9.  it  can  be  seen  that  lower  frame 
assembly  96b  further  includes  a  pair  of  horizontally 
extending  arms  110  which  engage  pilings  36.  These 
arms  function  to  carry  the  lateral  thrust  of  the  discharge 
means  as  the  scour  unit  is  rotated  in  a  manner  presently 
to  be  described.  Additionally,  face  member  104  (FIG.  T) 
is  provided  with  a  lower,  transversely  extending  mem¬ 
ber  104<2  which  is  disposed  below  elbow  97  to  eliminate 
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the  rooDOMl  movement  indoced  by  the  discbarge 


10 


15 


20 


30 


A.  best  ssa  by  referring  to  FIGS.  8  and  9.  actuating 
generaUy  designated  by  the  numer^lli  ^ 

opciaWy  connected  to  each  tubular  member  93  to  afto 

degree  of  roobon  of  the 
.bout  a  fixed  point  P  on  the  support 
.ctuatmg  mema  lU  comprises  a  hydraulic  cylmdw  m 
within  which  a  connector  arm  and  piston  a^bly  116 
is  earned.  Hydraulic  cylinder  lU  is  pivotally  mtercon- 
Dccted  with  frame  96fi  at  a  pomt  117  and  connector  ^ 

Xl«  is  pivotally  interconnected  with  a  broket  11  • 
Bracket  118  is.  in  turn,  connected  to  nng  120  whi^  is 
connected  to  tubular  member  93  at  a  locanon  above 
nng  10*.  The  inboard  end  of  bracket  118  is 
competed  to  frame  9«.  by  a  pivot  pm  122  located  at 
pomt  P.  Reciprocal  movement  of  the  piston  mtenoriy 
blinder  11*  urges  telescopic  movement  of  connector 
Mm  116  inwardly  and  outwardly  of  the  cylinder  11  . 
which,  m  mm.  imparts  pivotal  movement  to  bracket 
118  about  point  P.  This  results  in  an  arcuate  movement 
of  mbular  member  93  in  the  manner  illustrated  m  rlU. 

^  The  scouring  units  80  of  this  laner  embodiment  are 
typically  held  in  place  against  pilings  36  by  cnains  C 
(91)  of  the  same  character  shown  in  FIGS.  1  anC  - 
Similarlv.  the  units  are  raised  and  lowered  oy  chmn 
hoists  of  a  character  well  known  to  those  skilled  in  the 

In  operating  the  various  embodiments  of  the  inven¬ 
tion  in  accordance  with  the  method  of  the  present  in- 
vennor.  the  scour  units  are  typically  activated  twice 
daily  bv  an  automatic  control  system.  The  control  sys¬ 
tem  is  programmed  to  initiate  system  operation  at  the  -5 
optimum  time  to  ensure  maximum  sediment  flushing  by 
ambient  currents.  Typically,  the  individual  scour  umts 
are  seoucntially  activated,  beginning  with  the  further¬ 
most  upstream  unit.  Tne  duty  cycle  for  each  scour  unit 
can  be  adjusted  to  accommodate  differences  in  sediment  40 
scour  properties  and  the  arc  of  scan  rotation.  In  addi- 
oon.  the  rate  of  scan  of  each  scour  unit  can  be  adjusted. 
Tests  have  shown  that  highly  cohesive  sediments  re¬ 
quire  longer  duty  cycle  times  and  slower  scan  rates  to 
ensure  proper  scouring  action.-  As  previously  men-  45 
tioned.  by  adjusting  the  location  of  the  intake  means, 
low  salinttv  water  can  be  used  for  scouring. 

"gbe  scouring  action  of  the  turbine  fan  discharge 
serves  to  resuspend  the  thin  layer  of  sediment  newly 
deposited  during  the  previous  flood  tide.  The  frequent  50 
operaoon  of  the  scour  system  reduces  input  power  re- 
quinancnts  and  mmimizes  the  turbidity  produced  by  the 
discharge  plume.  In  addition,  the  shon  residence  time 
of  the  newly  deposited  sediments  minimizes  adsorption 
of  undesirable  toxins.  Both  features  serve  to  minimize 
poienual  adverse  environmental  impacts  to  water  quai- 
ity. 

Having  now  described  the  invennon  in  deiad  in  ac¬ 
cordance  with  the  requirctncnis  of  the  patent  sututes. 
those  skilled  in  this  an  will  have  no  difficulty  in  making 
changes  and  tnodificanona  in  the  individual  pans  or 
ihctf  relauve  assembly  in  order  to  meet  specific  require¬ 
ments  or  conditions.  Such  changes  and  modifications 
mav  be  made  without  departing  from  the  scope  and 


spirit  of  die  invention,  as  set  forth  in  the  following 

claims. 

Wc  claim: 

1.  An  apparatus  for  preventmg  sedimentation  build¬ 
up  in  harbor  berthing  area  of  the  character  having 
wharf  fronting  the  berthing  area,  comprising  a  pluW 
of  submereerd,  scour  units  disposed  adjacent  “‘d 
in  side-by-side  relationship,  each  said  scour  unit  com- 

vertical  mbular  member  having  a  venical  axis 
including  an  upper  end  and  a  lower  end: 

(b)  inlet  means  connected  proximate  said  up^r  end 

of  said  tubular  member  for  directing  the  flow 
water  into  said  tubular  member;  • .  •  i 

(c)  turbine  means  mounted  interiorly  of  »d  inlet 
means  for  controUably  drawing  water  mto  said 
inlet  means,  said  mrbine  means  includmg  a  routing 
fan  and  hydraulic  motor  means  for  controUably 
rotating  said  fan; 

(d)  discharge  means  connected  proamaie  the  lower 
end  of  said  tubular  member  for  directing  *e  flow 
of  water  outwardly  therefrom  in  a  generally  hon- 
zomal  direction  to  produce  a  near  bottom,  scour 
panem  of  a  predetermined  configuratim.  smd 
scour  units  being  arranged  along  said  wharf  so  that 
said  scour  patterns  produced  by  said  discharge 
means  overlap  one  another; 

(e)  means  for  rotating  each  said  discharge  means 
about  the  venical  axis  of  each  said  tubular  mcmocr. 
whereby  each  said  discharge  means  produces  a 
scour  pattern  spanning  about  180  degre«; 

(0  operating  means  mounted  on  said  whm  oper- 
abiy  iniercomiecied  with  each  said  hydraulic 
motor  means  for  providing  hydraulic  fluid  for  op¬ 
erating  each  said  hydraulic  motor  means;  and 
(g)  means  mounted  on  said  wharf  for  raising  and 
lowering  each  said  tubular  member.  . 

2.  .A.n  apparatus  for  preventing  sedimentation  build- 
UP  in  a  harbor  berthing  area  of  the  character  having  a 
wharf  fronting  the  berthing  area,  compnsmg: 

(a)  a  venical  tubular  member  having  a  venical  axis 
including  an  upper  end  disposed  below  the  surface 
of  the  benhing  area  and  a  lower  end; 

(b)  iniei  means  connected  proximate  said 
of  said  tubular  member  for  directing  the  flow  ot 
subsurface  water  into  said  tubular  member; 

(c)  mrbine  means  mounted  interiorly  of  said  inlet 
means  for  controUably  drawing  water  mto  said 

inlet  means;  .  , 

(d)  discharge  means  connected  by  means  of  a  single 
45  degree  to  125  degree  turning  elbow  proximate 
the  lower  end  of  said  tubular  member  for  directing 
the  flow  of  water  outwardly  therefrom  in  a  gener¬ 
ally.  horizontal  direction  to  produce  a  near  bottom 
scour  pattern  of  a  predetermined  configuration. 

(e)  fully  submerged  means  for  rotating  said  discharge 

means  about  the  vertical  axis  of  said  mbular  mem¬ 
ber,  and  . 

(f)  means  mounted  on  said  wharf  for  raising  and  low¬ 
ering  said  tubular  member.  ^  •  i  ^ 

3.  An  apparatus  as  defined  in  claim  2  further  includ¬ 
ing  operating  means  mounted  on  said  wharf  for  control- 
lab  I  y  operating  said  turbine  means. 
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Crater-sink  fluidization  has  been  used  to  recover  sediment  from  harbors  for  the  past  couple  of 
decades.  Can  this  method  bypass  sediment  around  dams  and  return  it  to  eroded  beaches,  solv¬ 
ing  two  major  problems?  Here's  how  it  might  happen  in  Southern  California. 

Tum-of‘the-century  photographs  show  an 
immense  beach  near  Oceanside,  Calif., 
about  1000  ft  wide.  Today,  no  beach  in  the 
area  has  a  width  of  more  than  a  few  hun¬ 
dred  feet  Many,  in  fact  are  totally  denuded 
of  sand  during  the  winter  months. 

An  eight-year,  $6  million  study  of  San 
Diego  County’s  coastline,  completed  by  the 
U.S.  Army  Corps  of  Engineers  in  1990,  con¬ 
firmed  that  the  beaches  are  in  dire  need  of 
nourishment  One  of  the  Corps’  recommen¬ 
dations  involves  restoring  the  equilibrium 
in  the  area’s  littoral  cells  (the  coastal  re¬ 
gions  defined  by  the  limits  of  beach  sedi¬ 
ment  movement)  by  augmenting  the  sand 
supply.  One  possible  source  of  this  sand 
would  be  to  bypass,  or  transfer,  sediment 
trapped  behind  area  dams  to  the  beaches. 

Due  to  the  steep  topography  of  the  region, 
dams  are  situated  near  the  coast,  making 
the  suggestion  feasible. 

At  the  Sciipps  Institution  of  Oceanogra¬ 
phy,  La  Jolla,  Calif.,  research  is  under  way 
with  the  California  Water  Authority  and  the 
San  Diego  Association  of  Governments  to 
discern  the  best  method  of  dam  bypassing 
for  the  area.  One  method  under  considera¬ 
tion  is  crater-sink  fluidization.  The  key  ad¬ 
vantage  of  this  type  of  hydraulic-extraction 
system  is  its  fixed-place  installation,  which 
allows  automation.  This  system  would  col- 


lect  coarse-grmn  sediments  from  a  top-set 
delta  by  means  of  mobile  fluidizer  pipes 
that  inject  high-pressure  water  into  the  sed¬ 
imentary  mass. 

During  fluidization,  the  weight  of  sand 
particles  is  balanced  by  the  upward  flux  of 
fluid  momentum  via  weight  drag.  As  the 
body  of  sand  dilates,  there  is  an  almost  to¬ 
tal  loss  of  static  shear  strength.  This  loss  al¬ 
lows  the  slurry  to  readily  flow  down  a  shal¬ 
low  slope  to  a  centrally  located  eductor 
crater,  formed  by  excavation  of  sediments 
with  a  jet  pump  (see  Fig.  1).  From  here, 
the  pipes  would  pump  sediments  through  a 
hydraulic  pipeline  to  the  eroding 
streambed  directly  downstream  of  the  dam, 
or  to  a  collection  point  for  overland  trans¬ 
port  to  eroding  beaches. 

To  understand  the  efficiency  of  this 
method,  it  helps  to  take  a  look  at  the  prob¬ 
lems  facing  the  coastal  area  between  Point 
La  Jolla  and  Dana  Point,  Calif.,  which  com¬ 
prises  the  Oceanside  littoral  cell  (Fig.  2). 
Dam  building  in  this  region  began  in  1922. 
The  most  recent  dam  was  completed  in 
1977.  Altogether,  there  are  21  catch  basins 
and  nine  major  dams  with  reservoirs. 

A  number  of  streams  and  rivers  in  this 
littoral  cell  empty  into  the  ocean  in  the 
neighborhood  of  Oceanside,  Calif.  The 
delta  deposits  from  the  Santa  Margarita 


and  the  San  Luis  Rey  rivers,  as  well  as  nu¬ 
merous  creeks,  once  contributed  an  aver¬ 
age  sediment  load  estimated  at  250,000  cu 
yd  per  year  to  the  foreshore. 

The  entire  coastal  area  in  Southern  Cali¬ 
fornia  is  typical  of  a  so-called  collision 
coastline.  The  steep  topography  of  the 
nearby  mountain  range  and  narrow  conti¬ 
nental  shelf  are  the  result  of  tectonic  pro¬ 
cesses  involving  collisions .  between  the 
oceanic  and  continental  plates.  The  semi- 
arid  climate  combined  with  the  steep  hy¬ 
draulic  grading  results  in  high  sediment 
transport  rates  down  streams  and  rivers 
during  seasonal  rainfall. 

Construction  of  the  Oceanside  Harbor- 
Del  Mar  boat  basin  in  1942  exacerbated 
beach  erosion  by  intercepting  the  souther¬ 
ly  littoral  drift.  Sediment  trapped  by  the 
harbor  is  now  dredged  on  a  regular  basis 
and  returned  to  littoral  cell  south  of  the 
harbor  entrance.  This  results  in  no  net 
long-term  effect  on  the  sediment  supply  to 
the  littoral  cell,  but  the  episodic  nature  of 
the  dredge  cycle  does  contribute  to  cycles 
of  beach  erosion,  called  erosion  waves. 

Bypassing  is  better _ 

The  principal  cause  of  beach  erosion  in  this 
area  lies  in  the  insufficient  supply  of  beach¬ 
sized  sand  to  replace  sand  transported 
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south  along  the  coast  by  waves  and  wave* 
driven  currents.  From  1922  to  1976,  an  esti¬ 
mated  14.5  million  cu  yd  of  sediment  was 
lost  down  the  La  Jolla  submarine  canyon. 

As  the  foreshore  recedes  landward, 
waves  act  directly  on  the  base  of  the  sea 
cliffs.  The  beach  communities  south  of 
Oceanside  suffered  millions  of  dollars  of 
property  damage  from  cliff  failures  during 
the  1978,  1980  and  1983  winter  storms. 
When  large  blocks  of  sea  cliff  collapse,  tens 
of  thousands  of  cubic  yards  of  potentially 
new  beach  sediments  enter  the  near-shore 
zone.  But  intense  development  of  the  sea 
cMs  precludes  this  as  an  alternative  for 
beach  nourishment 

Sandy  beaches  in  equilibrium  with  in¬ 
coming  waves  are  the  most  desirable 
means  of  wave  energy  dissipation  and 
shoreline  protection.  Structural  shore  pro¬ 
tection  and  other  remedial  measures  taken 
so  far  treat  only  the  symptoms  of  beach 
erosion.  Building  seawalls  at  the  base  of 
the  sea  cliffs  protects  only  the  property  di¬ 
rectly  behind  them.  In  addition  to  being 
costly,  conventional  seawalls  accelerate 
erosion  rates  by  increasing  the  reflectivity 
of  the  shoreline.  This  causes  an  increase  m 
the  rate  at  which  sand  travels  cross  shore 
and  down  slope. 

The  building  of  groins  and  jetties  in¬ 
hibits  littoral  drift  by  creating  pocket 
beaches.  These  beaches  don't  stop  regional 
erosion,  they  only  defer  it  to  the  beaches 
down  drift  of  the  structures.  No  new  sedi¬ 
ment  is  added  to  the  foreshore.  The  dwin¬ 


dling  supply  of  available  sand  is  merely  re¬ 
arranged. 

Bypassing  sediment  around  area  dams 
would  not  only  nourish  beaches,  but  would 
help  mitigate  many  of  the  other  adverse  en¬ 
vironmental  impacts  of  dam  building.  For 
instance,  the  streambeds  below  the  dams 
were  in  equilibrium  with  the  seasonal  flows 
of  sediment-laden  water.  With  no  new  sedi¬ 
ment  introduced  from  upstream  sources, 


The  three  systems 


could  completely 


replenish  the 


14.5  million 


sand  deficit 


Z4  years 


the  streambeds  can't  resist  erosion  caused 
by  elevated  bottom  stresses  arising  from 
the  excess  kinetic  energy  of  strong,  clear 
water  flows.  High  river  flow  rates  occur  ei¬ 
ther  during  periods  of  heavy  rainfall  or 
high  water  discharge  rates  from  the  dams. 

Another  adverse  environmental  impact 
in  the  Oceanside  littoral  cell  can  be  found 
along  the  San  Luis  Key  River  below  the 
Lake  Henshaw  Dam.  Bank  erosion  has 
damaged  numerous  trestles  and  utility 
poles  along  state  Highway  76  and  has  un¬ 
dermined  the  roadbed  in  several  places. 
Reintroducing  bypassed  sediment  directly 
to  the  streambeds  below  the  dams  would 


mitigate  much  of  this  erosion. 

Also,  deposition  of  sediment  in  the 
reservoirs  has  decreased  reservoir  capaci¬ 
ty,  Lake  Mead,  on  the  Colorado  River,  traps 
an  estimated  3  cu  mi  of  sediment  every  100 
years.  Early  dam  designers  recognized  this 
problem  and  specified  sediment  traps  up¬ 
stream  of  the  dams.  Nevertheless,  large 
sediment  deposits  still  form  in  the  upper 
reaches  of  the  reservoirs. 

These  deposits  can  cause  serious  prob¬ 
lems,  even  before  the  useful  capacity  of  the 
reservoir  is  seriously  threatened.  Newly 
formed  deposits  place  tremendous  static 
loads  on  the  native  soils  of  the  original 
riverbed,  resulting  in  soil  subsidence  with 
costly  consequences.  The  static  loading  of 
the  sedimentary  deposits  behind  Lake  Hen¬ 
shaw  has  caused  parts  of  the  Warner  Val¬ 
ley  to  subside  as  much  as  0.5  ft.  This  subsi¬ 
dence  produced  large  cracks  in  the  dam's 
earthen  core,  reducing  the  reservoir's  wa¬ 
ter  storage  capacity  to  only  one-third  of  its 
original  design  value. 

Crater-sink  fluidization _ 

In  choosing  the  best  method  of  dam  by¬ 
passing  for  the  area,  there  are  four  aspects 
to  consider:  (1)  the  extraction  of  the  sedi¬ 
ments  from  their  at-rest  position  at  the 
reservoirs;  (2)  the  transportation  of  sedi¬ 
ments  around  the  dams;  (3)  the  transporta¬ 
tion  of  sediments  from  the  dam  site  to  the 
beaches  or  eroding  areas;  (4)  the  disper¬ 
sion  of  sediments  within  the  erosion  areas. 

To  effectively  extract  sediments  from 
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the  reservoir,  the  system  must  adapt  to  the 
depositional  pattern  within  the  reservoir. 
Much  evidence  suggests  that  in  most 
reservoirs  this  pattern  follows  a  character¬ 
istic  delta  pattern  (see  Fig.  1).  Surveys  of' 
Lake  San  Marcos  in  San  Diego  County,  for 
instance,  demonstrate  such  a  pattern. 

Typically,  these  deltas  consist  of  three 
distinct  depositional  beds:  the  top  set,  fore 
set  and  bottom  set  These  deposits  eictend 
along  the  old  river  channel  and  may  reach 
the  base  of  the  dam.  Grain-size  matching, 
the  main  criterion  in  selecting  a  bypass 
site,  is  helped  by  the  natural  horizontal 
sorting  within  a  depositional  delta. 

Top-set  beds  consist  of  the  coarsest 
grain-size  sediments,  and  their  packing 
yields  the  lowest  percentage  by  weight  of 
water.  Generally,  these  beds  most  nearly 
match  the  grain  sizes  of  the  beach  and 
streambed  sediments  and  are  the  preferred 
nourishment  material.  Fore-set  beds  are  in¬ 
termediate  in  grain  size.  Core  analysis  must 
determine  their  suitability  as  beach  nourish¬ 
ment  material.  Bottom-set  deposits  consist 
primarily  of  clays  and  fine  sand  and  are  un¬ 
suitable  for  beach  nourishment,  though  they 
have  value  as  topsoil  replacements. 

This  depositional  pattern  exhibits 
some  variation  with  annual  water-level 
changes.  However,  the  pattern  is  suffi¬ 
ciently  predictable  to  allow  consideration 
of  a  fixed-place  hydraulic-extraction  tech¬ 
nique,  known  as  crater-sink  fluidization,  to 
bypass  sediment  from  the  top-set  deltas. 
Where  road  and  railroad  networks  render 
dam  sites  accessible,  a  combination  of 
conventional  earth-moving  technologies 
(dredges,  bulldozers,  clam  shells,  dump 
trucks  and  conveyor  belts)  must  also  be 
considered. 

The  crater-sink  fluidization  system  has 
several  advantages.  It  can  be  run  more  or 
less  continuously  in  order  to  bypass  the  ex¬ 
isting  top-set  delta,  then  operated  as  need¬ 
ed  to  bypass  sediments  that  accumulate  on 
the  top-set  portion  of  the  delta  during  sub¬ 
sequent  depositional  events.  Manpower 
costs  are  reduced  since  the  natural  dynam¬ 
ics  of  deposition  within  the  reservoir  trans¬ 
port  the  sand  to  the  system.  In  addition, 
sediment  can’t  accumulate  undesirable 
anaerobic  bacteria  or  toxicity  when  ifs  by¬ 
passed  at  the  same  rate  it  is  deposited. 

The  dimensioning  of  a  crater-sink  flu¬ 
idization  is  primarily  dependent  on  the  by¬ 
passing  transport  rate,  the  capture  volume 
and  the  median  grain  size.  The  system  is 
sized  so  that  bypassing  transport  rates  are 


sufficient  to  match  the  deposition  rate  in 
the  top-set  delta.  Additionally,  the  com¬ 
bined  capture  volume  of  the  crater  and  flu¬ 
idization  trenches  should  be  large  enough 
to  intercept  the  yield  of  a  major  deposition¬ 
al  event,  which  would  occur  during  a  50- 
year  storm. 

The  annual  transport  rate  into  Lake  San 
Marcos  reservoir  during  1980,  for  instance, 
was  200  cu  yd  per  hour.  A  similar  solids 
fraction  pumping  capacity  would  be  re¬ 
quired.  The  estimated  top-set  flux  of  the  ex¬ 
treme  flood  event  of  Feb.  16, 1980,  is  27,000 
cu  yd.  This  capture  volume  could  be  gener¬ 
ated  by  excavating  five  400  ft  long  open  flu- 
idizers  trenches  to  a  depth  of  15  ft  In  the 
design  of  our  proposed  fluidizer,  we  use  a 
flexible,  noncollapsing  hose  that  workers 
can  retrieve  after  excavation.  This  would  al¬ 
low  the  excavation  of  each  fluidization 
trench  to  be  accomplished  with  the  same 
mechanical  components.  Past  experience 
has  shown  that  the  systems  are  sensitive  to 
debris-fouling  in  the  fluidization  trenches. 
So  it  is  important  that  the  mechanical  com¬ 
ponents  are  retrievable  for  defouling. 

Value  engineering _ 

The  precise  economics  of  bypassing  sedi¬ 
ment  depend  on  site  factors  unique  to  each 


FIG.  1 .  DAM  BYPASS  SYSTEM 


SINK  FLUIDIZATION  SYSTEM  POSITIONED  OVER  A 
DEPOSITION  DELTA. 


FIG.  2.  OCEANSIDE  LITTORAL  CELL 


reservoir.  However,  we  can  make  a  reason¬ 
able  estimate  of  the  costs  using  convention¬ 
al  means,  based  on  numbers  fi*om  a  major 
beach  nourishment  project  undertaken  by 
the  Corps  of  Engineers  from  November 
1981  through  May  1982,  at  the  San  Luis 
Key  river.  TTie  Corps  extracted  923, 090  cu 
yd  of  sediment  from  borrow  pits,  then  dry 
hauled  it  by  truck  7  mi  to  nourishment 
sites  south  of  the  Oceanside  Harbor.  The 
process  cost  $2.71  per  cubic  yard.  Assum¬ 
ing  a  5%  annual  inflation  rate,  a  similar 
nourishment  project  today  would  cost  $4.41 
per  cubic  yard. 

The  average  distance  from  a  reservoir 
in  the  Oceanside  littoral  cell  to  the  beach  is 
30  mi.  Extra  transportation  costs  would  be 
15-30C  per  mile.  This  would  increase  unit 
costs  anywhere  from  $7.86  to  $11.31  per 
cubic  yard.  Fully  restoring  the  beaches  in 
the  Oceanside  littoral  cell  would  require 
14.5  million  cu  yd  of  sediment  at  a  cost  of 
$114  million-$164  million. 

Crater-sink  systems  must  operate  on  at 
least  three  major  reservoirs  in  order  to  by¬ 
pass  the  required  14.5  cu  yd  of  sediment. 
The  three  most  promising  sites  are  Lake 
Henshaw,  Vail  Lake  and  the  San  Marcos 
reservoirs.  Assuming  a  10-year  system  life 
span,  and  that  each  system  will  bypass 
300,000  cu  yd  of  sediment  each  year,  unit 
capital  costs  would  be  about  85<:  per  yard 
of  sediment.  The  annual  maintenance 
costs  are  estimated  at  5%  of  initial  capital 
costs,  adding  5C  per  cubic  yard  of  sedi¬ 
ment.  Energy  costs  for  equipment  opera¬ 
tion  will  add  roughly  20C  per  cubic  yard  to 
the  cost  of  bypassing  the  sediment  around 
the  dams.  The  unit  cost  totals  $1.10  per 
cubic  yard. 

The  three  systems  would  completely  re¬ 
plenish  the  14.5  million  cu  yd  sand  defidt  in 
the  Oceanside  littoral  cell  in  24  years,  assum¬ 
ing  that  sand  is  continually  lost  in  sediment 
sinks  such  as  La  Jolla  Canyon  at  the  historic 
littoral  drift  rate  of  300,000  cu  yd  per  year.  If 
the  sediments  are  transported  hydraulically 
in  the  streambeds  with  natural  stream  flow, 
the  total  cost  of  this  nourishment  would 
come  to  $15.9  million.  This  compares  very  fa¬ 
vorably  with  the  $114  niillion-“$164  million 
estimated  for  bypassing  with  conventional 
earth-moving  equipment  t 
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BACKGROUND  OF  THE  INVENTION 
5 

Although  the  available  power  of  ocean  wind  wave  has  been  calculated 
at  2.5  X  10’^  W  globally,  present  wave  energy  recovering  technology,  is  only 
marginally  economical  even  along  a  high  energy  coastline  such  as  the  United 
Kingdom.  This  is  in  large  part  due  to  limited  available  energy  density,  which 
10  requires  large  numbers  of  energy  collectors  to  be  strung  out  along  a  given 
length  of  coastline  in  order  to  recover  economically  feasible  quantities  of 
wave  energy. 

A  number  of  attempts  have  been  made  to  focus  tidal  and  wave  energy 
into  a  limited  area  to  allow  for  efficient  conversion  for  power  production, 
15  including  those  disclosed  in  the  patents  of  Arpin  (U.S.  Patent  No. 
2,069,715),  Danel  (U.S.  Patent  No.  2,435,576),  Gay  (U.S.  Patent  No. 
2,441,759),  Southwick  (U.S.  Patent  No.  2,820,148),  Conn,  et  al.  (U.S. 
Patent  No.  4,034,231),  Borgren  (U.S.  Patent  No.  4,363,564)  and  Mehlum 
(Re  32054).  The  Danel,  Gay,  Southwick  and  Borgren  patents  each  describe 
20  the  use  of  sea  walls  of  various  geometries  which  are  designed  to  channel 
wave  energy  to  one  or  more  focal  points.  The  Mehlum  patent  discloses  a 
variation  on  the  solid  converging  sea  walls  or  channels  which  involves  the 
creation  of  a  diffraction  grating  formed  from  discreet  structures  arranged  in 
a  grid  pattern.  The  grid  pattern  becomes  a  Fresnel  lens,  focusing  the  energy 
25  by  diffracting  the  incoming  waves.  Finally,  the  Conn,  et  al.  patent  relies  only 
partially  on  redirection  of  wave  energy.  Instead,  this  device  is  a  V-shaped 
pier-like  structure  with  a  number  of  water  turbines  distributed  at  different 
points  of  the  structure.  Rather  than  bringing  the  waves  to  the  turbine,  or  the 
conversion  device,  the  turbines  are  distributed  to  utilize  a  section  of  the 
3  0  wave  front. 

A  common  concern  of  each  of  these  devices  is  the  sea  wall  or 
channel.  The  construction  of  permanent  structures  can  be  expensive  and 
potentially  environmentally  destructive  by  altering  the  local  ocean  current 
patterns.  Also,  changes  in  sea  topography  due  to  storms  and  erosion  can 
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cause  the  ideal  location  for  energy  collector  placement  to  shift,  rendering 
an  expensive  permanent  structure  less  effective. 

The  need  for  a  portable  reef  or  sea  wall  has  been  addressed  in  U.S. 
Patents  No.  2,069,715  of  Arpin  and  No.  3,886,751  of  Porraz  Jimenez 
5  Labora. 

The  artificial  reef  of  Arpin  is  limited  in  that  its  structure  is  relatively 
shallow,  either  in  a  triangular  or  semicircular  cross  section.  While  some 
modifications  of  wave  patterns  in  shallow  water  is  possible,  the  entire 
purpose  of  this  artificial  reef  is  to  provide  a  nucleation  point  to  accumulate 
10  sand  to  form  a  sandbar.  This  type  of  structure  would  not  be  effective  in 
providing  sufficient  wave  channeling  for  energy  conversion  purposes. 

The  aquatic  construction  modules  of  Poraz  Jimenez  Labora  come 
closer  to  meeting  the  needs  of  a  portable,  sturdy  structure  for  modifying 
water  flow  patterns.  However,  the  interlocking  edges  of  the  modules  must 
15  be  individually  aligned  and  linked  together  during  assembly  of  the  reef  or 
other  protective  barrier.  This  procedure,  while  suitable  for  construction  of 
temporary  levees  and  channels  well  in  advance  of  flooding,  where  rapid 
assembly  of  the  barrier  is  not  critical,  is  extremely  impractical  for 
construction  in  existing  bodies  of  water,  especially  in  the  presence  of  ocean 
20  swells,  using  conventional  floating  cranes  and  riggings. 

SUMMARY  OF  THE  INVENTION 

It  is  an  advantage  of  the  present  invention  to  provide  a  structural  unit 
25  for  formation  of  artificial  reefs  for  diverting  and  channeling  wave  energy  and 
water  in  general. 

It  is  another  advantage  of  the  present  invention  to  provide  an  artificial 
reef  structure  for  channeling  wave  energy  for  energy  collection. 

It  is  still  another  advantage  of  the  present  invention  to  provide  an 
3  0  optimized  artificial  reef  structure  for  enhancing  the  energy  of  a  collected 
wave. 

In  an  exemplary  embodiment  the  artificial  reef  of  the  present  invention 
is  constructed  from  a  plurality  of  structural  units  each  of  which  is  a 
tetrahedron-shaped  sandbag,  which  are  arranged  in  a  talus  pile.  With  the 
35  exception  of  the  tetrahedron-shape,  the  sandbags  themselves  are 


constructed  according  to  conventional  sandbag  technology,  using  the  same 
materials  and  assembly  techniques.  Typically,  sandbags  are  made  from 
polypropylene  or  polyester  fabric  bags  which  are  filled  with  sand.  A 
minimum  of  five  layers  thick  of  the  structural  units  is  required  in  the  talus 
pile  to  reduce  the  permeability  of  the  reef  to  less  than  five  percent.  This 
level  of  permeability  is  important  to  obtain  a  solid  body  response. 

The  reefs  are  formed  by  dropping  the  structural  units  in  the  desired 
locations.  The  tetrahedron-shape  of  the  structural  units  provides  a  self¬ 
interlocking  structure  which  can  be  placed  en  mass  on  the  sea  floor  without 
individual  alignment  of  each  bag. 

The  geometry  of  the  reef  structure  is  a  sloped  bar,  a  parallel  bar,  or 
a  V-shaped  bar,  or  combinations  thereof.  The  seaward  slope  is  important  in 
determining  the  amount  of  energy  which  can  be  focused,  while  the  parallel 
bar  or  V-shaped  bar  focuses  energy  to  discreet  locations.  The  parallel  bars, 
which  run  parallel  to  the  wave  front,  trap  standing  waves  between  the  bars 
to  reinforce  the  incoming  waves  with  a  maximum  energy  focused  where  the 
rays  of  trapped  waves  cross  the  wave  rays  of  the  incoming  swell. 

A  multiple  V-shaped  bar  configuration,  with  the  apex  of  the  V  pointing 
diagonal  to  the  incoming  waves  creates  Bragg  reflection  patterns  which 
cause- crossing  wave  rays  to  form  local  concentrations  of  wave  energy. 

BRIEF  DESCRIPTION  OF  THE  DRAWINGS 

Understanding  of  the  present  invention  will  be  facilitated  by 
consideration  of  the  following  detailed description  of  the  preferred 
embodiment  of  the  present  invention,  taken  in  conjunction  with  the 
accompanying  drawings  in  which  like  reference  numerals  refer  to  like  parts 
and  in  which: 

Figure  1  is  a  perspective  view  of  a  structural  unit  of  the  present 
invention,  partially  cut  away; 

Figure  2a  is  a  front  elevation  and  Figure  2b  is  a  side  elevation  of  an 
artificial  reef  according  to  the  present  invention; 

Figure  3  is  a  top  view  of  a  parallel  bar  arrangement; 

Figure  4  is  a  top  view  of  a  multiple  V-shaped  bar  arrangement; 
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Figures  5a-5c  are  plots  of  wave  enhancement  by  a  parallel  reef  at 
frequencies  of  0.85  Hz  (Figure  5a)  and  0.95  Hz  (Figure  5b)  with  a  reef 
spacing  of  125  cm  and  at  .70  Hz  with  a  reef  spacing  of  105  cm  (Figure  5c); 

Figure  6a  is  a  plot  of  the  power  spectrum  of  the  control  wave  from 
5  the  structure  of  Figure  5b; 

Figure  6b  is  a  plot  of  the  power  spectrum  of  the  enhanced  wave  from 
~  the  structure  of  Figure  5b; 

Figure  7  is  a  comparison  of  predicted  resonance  separation  for  a 
parallel  bar  reef  from  linear  theory  versus  measured  laboratory  results; 

10  Figure  8  is  a  plot  illustrating  the  blue  shift  hypothesis; 

Figures  9a-9c  are  plots  of  alterations  in  wave  characteristics  caused 
by  a  V-bar  reef  for  wave  frequencies  of  0.5  Hz  (Figure  9a),  0.7  Hz  (Figure 
9b)  and  for  a  wave  period  of  1.6  seconds  (Figure  9c). 

15  DESCRIPTION  OF  THE  PREFERRED  EMBODIMENT(S) 

As  illustrated  in  Figure  1 ,  the  structural  unit  2  is  a  sandbag  formed  by 
sewing  together  triangular  sheets  of  polypropylene  or  polyester  fabric, 
typically  woven,  to  form  a  tetrahedron-shaped  bag  within  which  sand  3, 
20  gravel,  or  a  sand/gravel  mixture  is  retained.  Generally,  any  type  of  filler  with 
a  similar  weight  and  conformability  may  be  used,  including  gels  or  liquids. 
For  simplicity,  further  descriptions  will  identify  the  filler  as  "sand",  which  can 
include  any  of  the  above-identified  fillers  and  their  equivalents.  The 
dimensions  of  each  structural  unit  2  may  vary,  however,  in  the  preferred 
25  embodiment,  structural  units  with  heights  of  4  feet  and  6  feet  are  intermixed 
to  form  the  desired  artificial  reef  structures.  In  Figures  2a,  2b  and  3,  the 
different  sizes  of  structural  unit  are  indicated  by  the  hashed  triangles  for  6 
feet  and  open  triangles  for  4  feet. 

In  Figure  2,  a  single  bar  structure  artificial  reef  formed  from  plurality 
30  of  structural  units  2  is  illustrated.  As  can  be  seen,  the  structural  units  have 
intermixed  sizes  which  are  randomly  distributed  by  dropping  the  structural 
units  over  the  location  at  which  the  reef  is  to  be  formed.  As  illustrated  in 
Figure  3,  the  bar  arrangement  does  not  consist  of  a  single  thickness  wall  of 
structural  units,  but  is  essentially  a  talus  pile  formed  by  dropping  the 
35  structural  units  around  the  target  point  to  form  a  bar  with  a  width 


corrBsponding  to  the  combined  widths  of  structural  units.  In  the  prefered 
embodiment,  a  reef  width  of  at  least  five  structural  units  is  needed  to  reduce 
the  permeability  of  the  reef  to  less  than  5%.  The  height  of  the  reef  Is 
determined  by  the  ocean  depth  at  the  location  of  formation  as  well  as  the 
range  of  wave  heights  and  wave  periods.  The  seaward  slopes  of  the 
artificial  reefs  are  important  in  determining  the  amount  of  energy  which  ban 
be  focused  on  a  given  energy  collector  by  the  given  reef  shape. 

The  optimal  wave  form  for  both  energy  recovery  and  shore  protection 
was  found  to  be  a  plunging  wave.  The  plunging  wave  was  obtained  when 
the  slope  of  the  reef  was  sufficiently  steep  to  ensure  a  Carrier-Greenspan 
parameter  less  than  0.4  according  to  the  equation: 


tan  3  > 


0.4g) 


(1) 


Here,  tan  p  is  the  slope  of  the  bar,  Hb  is  the  height  of  the  breaking  wave,  a 
is  the  radian  frequency  equal  to  2/r/T,  where  T  is  the  wave  period,  and  g  is 
the  acceleration  of  gravity.  Consequently,  the  bar  must  be  steeper  as  the 
wave  heights  increase  or  the  wave  periods  decrease.  If  the  reef  slope  is 
constrained  for  structural  regions  to  some  particular  angle  of  repose,  the 
tendency  for  plunging  waves  will  increase  as  the  wave  period  expands  or  the 
wave  height  decreases. 

To  accommodate  both  average  and  extreme  conditions,  multiple  reefs 
such  as  the  parallel  bar  of  Figure  3  and  V-shaped  reefs  of  Figure  4  may  be 
required.  Multiple  reef  configurations  have  the  added  advantage  of  providing 
unique  wave  structure  interactions  that  can  dramatically  enhance  the  local 
wave  height,  as  much  as  80%. 

In  the  case  of  the  parallel  bars  4  and  4'  of  Figure  3,  standing  waves, 
indicated  by  wave  rays  6,  may  become  trapped  between  the  bars  and 
reinforce  the  incoming  waves,  indicated  by  wave  rays  8.  For  non-normal 
incidence,  the  strongest  reinforcement  occurs  where  the  rays  of  trapped 


waves  6,  shown  by  the  dotted  lines,  cross  the  wave  rays  8  of  the  incoming 
swell,  indicated  by  the  solid  lines.  The  wave  energy  collection  devices 
should  be  placed  in  the  region  of  the  loci  30  for  maximum  energy  recovery. 
Other  wave  effects  resulting  from  the  placement  of  the  parallel  bars  include 
specular  reflected  waves,  shown  as  wave  rays  7,  and  non-specular  leaky 
waves,  shown  as  wave  rays  9. 

In  the  multiple  V-bar  configuration  illustrated  in  Figure  4,  there  is  the 
possibility  for  appropriate  spacings  of  the  V-bars  42  to  create  Bragg 
reflection  patterns.  These  complex  interference  patterns  result  in  crossing 
wave  rays  which  provide  local  concentrations  in  wave  energy,  indicated  by 
"bright  spots"  10,  after  multiple  reflections  between  the  various  elements  of 
the  V-bar  system.  The  bright  spots  1 0  in  this  complex  interference  pattern 
are  the  focal  points  of  the  incident  wave  beams  1 2  and  become  the  most 

efficient  wave  energy  collection  points. 

In  Figure  3,  wave  rays  are  diagrammed  over  a  parallel  bar 
configuration  consisting  of  4  and  4'.  It  has  been  found  that  for  a  given 
wave  direction  relative  to  the  reefs,  there  is  a  particular  frequency  at  which 
wave  enhancement  is  especially  vigorous.  Figures  5a-5c  provide 
comparative  data  of  wave  enhancement  due  to  the  parallel  reefs  at  three 
different  wave  frequencies.  In  these  plots,  the  solid  lines  (labeled  REEF  ) 
are  measured  by  wave  staffs  in  the  trapping  zone  1 2  between  the  parallel 
reefs  4  and  4',  while  the  dotted  lines  (labeled  "CONTROL")  correspond  to 
waves  at  a  comparable  depth  along  the  undisturbed  shore.  In  all  cases,  the 
heights  of  the  waves  are  increased  due  to  the  presence  of  the  parallel  reefs 
4  and  4',  while  at  the  same  time,  the  peaks  are  made  steeper  and  sharper. 

For  the  higher  frequency  waves  in  Figure  5a,  the  heights  have  been 
increased  by  nearly  87%  over  those  of  the  control  wave  heights  as  a  result 
of  the  energy  trapping  effects  of  the  parallel  reefs.  A  similar  enhancement 
is  seen  at  the  slightly  higher  frequency  illustrated  in  Figure  5b. 

For  more  complex  sea  states  at  lower  frequencies,  such  as  the 
conditions  of  Figure  5c,  the  parallel  reefs  are  shown  to  greatly  increase  the 
steepness  of  the  waves.  Again,  the  increase  in  wave  heights  relative  to  the 
undisturbed  shore  (control)  ranges  from  78-91%.  By  examining  the  power 
spectra  of  the  data  in  Figure  5c,  it  is  apparent  that  increases  in  wave  heights 


are  the  result  of  harmonics  which  are  enhanced  or  trapped  by  the  parallel 
reefs. 

Figure  6a  is  a  plot  of  the  power  spectra  in  the  control  area  of  the 
undisturbed  beach,  while  Figure  6b  shows  the  power  spectra  obtained 
between  the  two  parallel  reefs  4  and  4'  in  the  energy  trapping  zone  14. 
Comparing  the  two  plots,  the  harmonics  in  Figure  6b  are  signicantly 
enhanced  from  those  of  Figure  6a.  Because  these  harmonics  are  phase- 
locked  to  the  fundamental  frequency  of  the  Incoming  waves,  they  act  to 
reinforce  the  predominant  peak,  resulting  in  the  very  steep  faces  and 
increased  wave  heights  shown  in  Figure  5C. 

Initial  laboratory  testing  involved  direct  measurement  and  photography 
of  the  response  of  1-2.5  second  laboratory  waves  to  model  talus  pile  bars 
and  reefs.  These  model  reefs  were  constructed  of  concrete  tetrahedrons  of 
two  sizes,  2"  and  3".  In  certain  modeled  cases  that  involved  local  wave 
height  doubling,  the  interactions  induced  in  the  laboratory  from  incident 
waves  of  1  second  period  and  10  cm  amplitude  formed  jets  of  water  which 
rose  as  high  as  three  feet  in  the  air.  This  condition  was  observed  for  1  m 
spacings  between  the  parallel  bars  and  a  12°  wave  angle.  Video  data  of 
this  condition  was  combined  with  wave  staff  measurements  to  provide  the 
optimal  reef  separation  as  a  function  of  the  frequency  of  the  incoming 
waves.  This  relationship  is  plotted  in  Figure  7.  The  solid  line  16 
corresponds  to  the  prediction  based  on  linear  theory  for  the  formation  of  a 
trapped  wave  at  some  multiple  of  the  incident  wave  lengths,  as  given  by. 


s  = 


r 

2 


y^cosa 


(2) 

Here,  S  is  the  separation  between  the  bars  (the  width  of  trapping  zone  14) 
in  the  on/offshore  direction  and  a  is  the  approach  angle  of  the  incoming 
swells. 

Two  physical  phenomena  combine  to  create  the  vigorous  wave 
enhancement.  The  first  is  a  trapping  of  the  incident  wave  energy  between 
the  two  bars  4  and  4'  in  the  form  of  a  standing  wave.  If  this  were  the  only 
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mechanism,  the  data  in  Figure  7  should  allow  the  solve  line  prediction  of 
linear  theory.  Therefore,  another  mechanism  is  clearly  at  work.  That 
mechanism  is  the  creation  of  a  double  caustic  as  a  result  of  the  change  in 
the  wave  lengths. 

5  The  wave  length  is  reduced  by  dissipative  effects  over  the  porous 

breakwater  surfaces  as  the  incoming  and  reflected  waves  refract  over  the 
-  tops  of  the  bars.  This  is  illustrated  schematically  in  Figure  8,  in  which  an 
undisturbed  incident  wave  20  is  compared  against  a  wave  22  propagating 
over  the  reef  4.  The  phase  speed  of  the  incident  wave  20  is  Cq.  Phase 
)  speed  of  the  waves  after  passing  over  the  reef  4  is  C',  which  is  slowed 
down  due  to  the  effects  of  friction  over  the  reef. 

Since  the  phase  speed  is  simply  a  ratio  of  the  radian  frequency  to  the 
wave  number,  the  phase  speed  can  only  be  reduced  by  either  reducing  the 
frequency,  or  by  increasing  the  wave  number  through  wave  length 
5  reduction.  (The  wave  number  is  2tt  over  wavelengths.)  However,  according 
to  the  principle  of  stationary  phases 


20 


9k  +  Va  =  0 

ar 


in  steady  state:  9k  =  0  so  Va  =  0, 


(3) 

(4) 


so  If  the  incoming  wave  train  is  steady,  the  frequency  cannot  change  in 
space,  i.e.. 
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C'  (  Co  so  k'  >  ko. 


(5) 


Therefore,  the  wave  number  must  be  increased  due  to  friction,  and  the 
wave  length  is  correspondingly  reduced.  This  reduction  in  wave  length  can 
be  estimated  by  applying  the  Lorentz  Principle: 


(Cd„  =  C,  [ul  )  u  ~  afu 


(6) 


-9- 


hence, 


k'^kA 


"  (7) 

The  Lorentz  Principle  linearizes  the  drag  effects  of  the  reef  in  terms  of 
5  units  of  virtual  mass  force  or  first  order  wave  pressure.  The  resulting 
formulation  gives  rise  to  a  blue-shift  in  wave  number  k'  proportional  to  a 
dissipation  coefficient  f. 

Because  of  this  change  in  wave  number  or  wave  length,  the  data  in 
Figure  7  did  not  obey  linear  theory.  Instead,  the  optimal  separation  distance 
10  in  space  14  of  Figure  2  is  proportional  to  the  dissipation  coefficient.  The 
result  is  a  double  caustic  at  the  point  indicated  by  reference  numeral  30, 
formed  by  the  constructive  interference  between  the  incoming  swell  waves, 
the  trapped  standing  waves  32  between  the  reefs  4  and  4'  and  the  refracted 
waves  with  altered  wave  lengths.  This  constructive  interference  actually 
15  occurs  between  the  two  bars  4  and  4',  although  generally  closer  to  the 
Inshore  bar  4.  This  positioning  of  the  double  caustic  30  insures  that  the 
optimal  location  34  of  the  wave  energy  collectors  occurs  in  the  relatively 
deep  area  between  parallel  bars  4  and  4',  rather  than  over  them,  thereby 
limiting  the  possibility  of  potential  damage  to  the  moorings  of  the  collectors. 
20  Measurements  of  the  onshore  and  inshore  wave  heights  showed  wave 
height  nearly  doubling  due  to  the  energy  trapping  and  an  80%  reduction  of 
wave  height  Inshore  of  the  parallel  bar  system.  Therefore,  the  double 
caustic  30  offers  both  enhanced  wave  energy  recovery  as  well  as  an 
effective  mechanism  for  providing  shore  protection. 

25  The  other  wave  height  enhancing  reef  geometry  tested  was  the  V-bar 

system  illustrated  in  Figure  4.  The  vertex  40  of  the  V-bar  42  is  arranged 
seaward  with  the  axis  before  of  the  V  orthogonal  to  the  local  depth 
contours.  Incoming  waves  are  focused  near  the  vertex  40  of  the  V-bar  42 
through  local  difraction,  acting  to  bend  wave  rays  perpendicular  to  the  local 
30  surfaces  of  the  bar  and  thereby  creating  a  convergence  of  energy  at  the 
vertex  40.  Time  series  measured  for  the  wave  enhancing  effects  of  the  V- 
bar  are  shown  in  Figures  9a-9c.  In  these  plots,  the  curves  labeled  CONTROL 


-10- 


correspond  to  the  control  area  on  the  undisturbed  beach.  The  time  series 
denoted  with  REEF  were  measured  in  the  breaker  zone  of  the  V-bar.  Waves 
having  a  low  frequency  on  the  order  of  0.5  Hz  (Figure  9a)  experience  only 
modest  increase  in  the  breaker  heights,  on  the  order  of  25-30%,  as  a  result 
5  of  the  placement  of  the  V-bar  arrangement.  At  higher  frequencies,  in  Figures 
10  and  11,  the  breaker  heights  are  nearly  doubled  by  the  convergerice 
effects  of  the  V-bar.  Again,  the  wave  form  is  steepened  with  the  higher 
frequency  lumps  (surface  chop)  measured  in  the  control  area  not  being 
eliminated  from  those  breakers  in  the  V-bar. 

3^0  In  some  cases  at  0.7  Hz,  shown  in  Figure  9b,  the  V-bar  more  than 

doubles  the  incident  wave  height.  Vertex  angles  ranging  from  30-160°  were 
evaluated.  The  optimal  range  of  vertex  angles  for  the  V-bar  was  found  to 
be  120°.  The  actual  selection  of  vertex  angles  at  a  given  site,  however,  will 
be  controlled  by  the  beach  slope,  with  steep  beaches  requiring  larger  vertex 
15  angles.  The  large  vertex  angles  on  the  steep  beaches  will  cause  the  wave 
to  shoal  farther  offshore  with  characteristics  desirable  for  efficient  energy 
recovery. 

In  order  to  obtain  the  highly  dissipative  plunging  breakers  that  are 
optimal  for  energy  recovery  and  shore  protection  alike,  the  reef  slope  and  the 
20  water  depth  over  the  reef  must  be  tailored  to  suit  4  particular  sight 
parameters:  (1)  the  wave  height  and  (2)  wave  period  resulting  from  the  local 
wave  climate,  H  and  T;  (3)  the  local  tidal  elevation  77;  and  (4)  the  wave 
direction  a.  The  reef  slope  must  be  adjusted  according  to  some  prevalent 
wave  height/  period  and  tidal  range.  To  take  advantage  of  wave  height 
25  enhancing  techniques  such  as  the  parallel  bar  and  the  V-bar,  the  design  must 
also  be  tuned  to  the  most  prevalent  wave  directions  a.  Once  these  site- 
specific  design  wave  conditions  have  been  established,  together  with  the 
tidal  conditions,  the  reef  geometry  can  be  posed  by  Equation  (1)  and 


s  = 


2 


)** 


cosa 


30 


(8) 
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The  wave  forces  on  the  structural  units  2,  which  act  to  roll  them,  vary 
in  response  to  local  wave  conditions  according  to  the  equation: 


(9) 

5 

where  c,  is  the  drag  coefficient,  c^,  is  the  virtual  mass  coefficient,  p  is  the 
density  of  water,  and  I  is  the  dimension  of  the  structural  unit  2. 
Experiments  were  conducted  in  the  wind/wave  channel  to  determine  the 
drag  and  virtual  mass  coefficients  in  the  above  equation.  During  these 
10  experiments,  frequencies  varied  0.73  Hz-1 .93  Hz.  Wave  heights  were  varied 
from  5-20  cm.  This  produced  variations  in  wave  lengths  corresponding  to 
deep  water  waves  for  the  0.73  Hz  waves  to  the  shallow  water  waves  for 
the  0.195  Hz  waves.  The  resulting  virtual  mass  (cj  and  drag  (p  ) 
coefficients  obtained  during  these  tests  are  compiled  in  Tables  la-c  and  2a-c 
15  for  deep  water  waves,  intermediate  water  waves  and  shallow  water  waves. 
Tables  1  a  and  2a  provide  the  conditions  under  which  the  tests  were  run. 
Tables  1b  and  2b  are  the  virtual  mass  coefficients  c^,  corresponding  to 
Tables  la  and  2a,  respectively.  Tables  1c  and  2c  are  the  drag  coefficients 
Cd  corresponding  to  Tables  la  and  2a. 


20 
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TABLE  1A 


Freq  (Hz) 

H{cm) 

k  (cm’’) 

^2 

u„,D/u 

U.T/D 

DeeD  Water 

0.732 

5 

0.0216 

0.216 

0.409  X  10^ 

0.167 

bar  V#  1^  •  " 

Waves 

0.732 

10 

0.0216 

0.216 

0.785  X  IO4 

0.322 

0.732 

15 

0.0216 

0.216 

1.128  X  10^ 

0.462 

2n  >1/2 

0.732 

20 

0.0216 

0.216 

1.408  X  10^ 

0.577 

Intermediate 

0.585 

5 

0.0142 

0.142  ' 

0.617  X  lO'^ 

0.316 

Water  Waves 

0.585 

10 

0.0142 

0.0142 

0.173  X  lO'^ 

0.600 

i_ 

0.585 

15 

0.0142 

0.142 

1.749  X  lO'^ 

0.895 

200<2n<1/2 

0.585 

20 

0.0142 

0.142 

2.203  X  10^ 

1.128 

Intermediate  0.391 

5 

0.0076 

0.0757 

0.835  X  lO'^ 

0.640 

Water  Waves 

0.391 

10 

0.0076 

0.0757 

1.756  X  10^ 

1.347 

1_  kh 

0.391 

15 

0.0076 

0.0757 

2.608  X  lO'^ 

2.000 

200<2n<1/2 

0.391 

20 

0.0076 

0.0757 

3.373  X  10^ 

2.580 

Shallow  Water 

0.195 

5 

0.0033 

0.0333 

1.027  X  10^ 

1.580 

Waves 

0.195 

10 

0.0033 

0.0333 

2.091  X  10^ 

3.216 

kh 

0.195 

15 

0.0033 

0.0333 

2.965  X  10^ 

4.560 

20  <  1/2000.195 

20 

0.0033 

0.0333 

3.609  X  lO'^ 

5.553 
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TABLE  1 B 


C„, 


run  1 

run  2 

run  3 

run  4 

Deep  Water  0.639 

0.643 

0.642 

0.645 

Waves 

0.699 

0.689 

0.699 

0.695 

kh 

0.704 

0.700 

0.704 

0.703 

2U  >  1/2 

0.763 

0.758 

0.752 

0.754 

Intermediate  0.636 

0.635 

0.634 

0.636 

Water  Waves 

0.635 

0.641 

0.634 

0.636 

1  kh 

0.694 

0.694 

0.690 

0.696 

200<2n<1/2 

0.706 

0.710 

0.719 

0.718 

Intermediate  0.859 

0.865 

0.873 

0.878 

Water  Waves 

0.811 

0.807 

0.810 

0.808 

1  kh 

0.745 

0.753 

0.751 

0.751 

200<2n<1/2 

0.803 

0.772 

0.781 

0.720 

Shallow  Water 

0.417 

0.419 

0.421 

0.421 

Waves 

0.408 

0.404 

0.409 

0.409 

kh 

0.448 

0.447 

0.453 

0.446 

in  <  1/2000.491 

0.481 

0.495 

0.491 
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TABLE  1C 


Cd  =  2C, 


run  1 

run  2 

run  3 

run  4 

Deep  Water  45.31 

45.58 

45.67 

46.17 

Waves 

26.29 

25.57 

25.16 

25.80 

kh 

18.41 

18.16 

18.41 

18.26 

2U  >  1/2 

15.99 

15.76 

15.48 

15.64 

Intermediate  24. 1 9 

24.45 

24.35 

24.11 

Water  Waves 

12.44 

12.66 

12.61 

12.69 

1  kh 

9.15 

9.16 

9.06 

9.22 

200<2n<1/2 

7.26 

7.37 

7.56 

7.58 

Intermediate  1 5.61 

15.72 

15.99 

16.24 

WaterWaves 

7.16 

7.12 

7.16 

7.11 

1  kh 

4.39 

4.46 

4.46 

4.45 

200<2n<1/2 

4.07 

3.79 

3.83 

3.29 

Shallow  Water 

3.16 

3.17 

3.19 

3.19 

Waves 

1.50 

1.48 

1.49 

1.51 

kh 

1.16 

1.19 

1.17 

1.18 

20  <  1/200 

1.06 

1.08 

1.04 

1.06 
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TABLE  2A 


Freq  (Hz) 

H  (cm) 

k  (cm  ’) 

Deep  Water 
Waves 

mi 

5 

10 

mi 

kh 

0.732 

15 

0.0216 

2n  >  1/2 

0.732 

20 

0.0216 

Intermediate  0.586 

5 

0.0142 

Water  Waves 

0.586 

10 

0.0142 

1_  kh 

0.586 

15 

0.0142 

200<2n<1/2 

0.586 

20 

0.0142 

Intermediate  0.391 

5 

0.0076 

Water  Waves 

0.391 

0.391 

10 

15 

0.0076 

0.0076 

200<2n<1/2 

0.391 

20 

0.0076 

Shallow  Water 

0.195 

5 

0.0033 

Waves 

kh 

0.195 

0.195 

10 

15 

0.0033 

0.0033 

2n  <  1/2000.195 

20 

0.0033 

kj  U„D/u  U„T/D 


0.108 

0.204  X  10^ 

0.334 

0.108 

0.393  X  104 

0.644 

0.108 

0.564  X  104 

0.924 

0.108 

0.704  X  104 

1.154 

0.071 

0.308  X  104 

0.632 

0.071 

0.586  X  10^ 

1.200 

0.071  ' 

0.874  X  10^ 

1.790 

0.071 

1.101  X  10^ 

2.256 

0.037 

90.417  X  104 

1.280 

0.037 

90.878  X  10 

2.694 

0.037 

91.304  X  104 

4.001 

0.037 

91.686  X  lO'^ 

5.160 

0.0167 

0.513  xIO^ 

3.161 

0.0167 

1.045  X  10^ 

6.432 

0.0167 

1.482x  104 

9.122 

0.0167 

1.804  X  10^ 

11.106 
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TABLE  2B 


C. 


run  1 

run  2 

run  3 

run  4 

Deep  Water 

0.676 

0.674 

0.678 

0.680 

Waves 

0.726 

0.726 

0.727 

0.732 

kh- 

0.754 

0.748 

0.761 

0.761 

2n  >  1/2 

0.749 

0.746 

0.754 

0.746 

Intermediate 

0.675 

0.688 

0.689 

0.691 

Water  Waves 

0.708 

0.685 

0.700 

0.709 

1  kh 

0.662 

0.656 

0.659 

0.663 

200<2n<1/2 

0.689 

0.689 

0.691 

0.697 

Intermediate 

0.814 

0.816 

0.818 

0.819 

Water  Waves 

0.792 

0.790 

0.792 

0.790 

1  kh 

0.716 

0.713 

0.714 

0.714 

200<2n<1/2 

0.671 

0.665 

0.672 

0.666 

Shallow  Water 

0.405 

0.401 

0.403 

0.400 

Waves 

0.367 

0.381 

0.383 

0.379 

kh 

0.406 

0.404 

0.407 

2fl  <  1/2000.416 

0.411 

0.426 

0.414 

TABLE  2C 


run  1 

run  2 

run  3 

run  4 

Deep  Water 

23.78 

23.77 

23.99 

24.21 

Waves 

13.91 

13.78 

13.75 

13.94 

kh- 

10.13 

10.01 

10.32 

10.34 

2n  >  1/2 

7.30 

7.22 

7.41 

7.52 

Intermediate 

14.57 

14.98 

15.13 

15.20 

Water  Waves 

8.43 

7.93 

8.25 

8.48 

1  kh 

4.41 

4.34 

4.37 

4.42 

200<2n<1/2 

3.57 

3.57 

3.60 

3.65 

Intermediate 

7.72 

7.64 

7.67 

7.70 

Water  Waves 

3.52 

3.50 

3.51 

3.49 

1  kh 

2.13 

2.12 

2.11 

2.12 

200<2n<1/2 

1.55 

1.53 

1.55 

1.55 

Shallow  Water 

1.55 

1.53 

1.53 

1.50 

Waves 

0.711 

0.697 

0.706 

0.692 

1^ 

0.542 

0.540 

0.544 

2n  <  1/2000.468 

0.466 

0.483 

0.470 
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The  virtual  mass  coefficients  were  found  to  be  surprisingly  small,  ranging 
from  0.4  to  0.9.  However,  the  drag  coefficients  were  found  to  be 
remarkably  large,  reaching  as  high  as  Cd  =  45.  In  order  to  provide  a  design 
safety  margin  to  ensure  that  the  reef  maintains  its  integrity  and  remains  in 
5  place,  the  most  extreme  values  in  drag  coefficient  and  virtual  mass 
coefficient  were  used  for  estimating  the  strength  requirements  of  a 
prototype  reef. 

To  satisfy  the  requirements  of  the  extreme  design  wave,  the 
minimum  structural  unit  count  for  the  talus  pile  according  to  Equation  (9) 
10  with  values  from  the  tables  is  estimated  to  be  860  bags  to  construct  a 
single  V-bar  reef,  as  shown  in  Figure  4,  or  1,825  structural  units  to 
construct  a  single  parallel  bar  field. 

The  structural  units  and  artificial  reefs  constructed  therefrom  of  the 
present  invention  effectively  concentrate  or  focus  the  incident  wave 
1 5  energy  at  discrete  locations  so  that  wave  energy  can  be  efficiently 

channeled  into  collectors  for  conversion  into  usable  power.  The  structural 
units  are  inexpensive  and  durable,  can  be  easily  installed  at  the  desired 
site,  and  can  be  resposition  without  major  demolition  and  reconstruction 
efforts,  as  compared  to  a  fixed  form  artificial  reef.  The  parallel  bar  and  V- 
20  bar  configurations  are  designed  to  concentrate  and  amplify  energy  for 
improved  efficiency  of  energy  conversion.  The  structural  units  and 
artifical  reefs  constructed  therefrom  may  also  be  used  for  protection  of 
the  shorleine  from  wave  damage  and  as  an  alternative  for  levee  and  flood 
barrier  construction. 

25  It  will  be  evident  that  there  are  addtional  embodiments  which  are 

not  illustrated  above  but  which  are  clearly  within  the  scope  and  spirit  of 
the  present  invention.  The  above  description  and  drawings  are  therefore 
intended  to  be  exemplary  only  and  the  scope  of  the  invention  is  to  be 
limited  solely  by  the  appended  claims. 
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CLAIMS 

1 .  An  artificial  reef  structure  for  focusing  ocean  wave  energy  on 
2  collection  devices  comprising: 

a  plurality  of  structural  units  disposed  in  at  least  one  talus  pile,  , 

4  each  structural  unit  comprising  a  sandbag  having  a  tetrahedron-shape 
-  with  four  corners  wherein  said  corners  substantially  interlock  one  said 
6  structural  unit  with  other  structural  units  of  said  plurality. 

2.  An  artificial  reef  structure  as  in  Claim  1 ,  wherein  said  sandbag 

2  comprises  a  durable  polymer  fabric  formed  into  a  tetrahedron-shaped  bag 
filled  with  sand. 

3.  An  artificial  reef  structure  as  in  Claim  1 ,  wherein  said  plurality  of 
2  structural  units  comprises  at  least  two  groups  of  structural  units,  each 

group  having  different  dimensions  from  other  groups. 

4.  An  artificial  reef  structure  as  in  Claim  1 ,  wherein  said  at  least  one 

2  talus  pile  comprises  two  parallel  bar  aligned  generally  parallel  to  a 

shoreline  and  separated  by  spacing  to  enhance  wave  height  of  incoming 
4  waves  within  said  spacing  at  a  location  of  said  collection  devices. 

5.  An  artificial  reef  structure  as  in  Claim  1,  wherein  said  at  least  one 
2  talus  pile  comprises  a  plurality  of  V-shaped  bars,  each  V-shaped  bar 

having  a  vertex  and  an  axis,  wherein  said  axis  is  disposed  orthogonal  to 
4  local  depth  contours. 

6.  An  artificial  reef  structure  as  in  Claim  5,  wherein  said  vertex  has  an 

2  angle  within  the  range  of  300-1600  and  said  angle  is  selected  according 
to  a  slope  of  an  adjacent  beach,  with  a  steep  beach  requiring  a  larger 
4  vertex  angle. 

7.  An  artificial  reef  structure  as  in  Claim  1,  wherein  said  at  least  talus 
pile  has  a  height  and  a  slope  selected  according  to  wave  height,  wave 
period,  local  tidal  elevation  and  wave  direction. 
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8.  An  artificial  reef  structure  as  in  Claim  1,  wherein  said  at  least  one 
2  talus  pile  has  a  width  comprising  at  least  five  structural  units. 

9.  A  structural  unit  for  an  artificial  reef  structure  comprising: 

2  a  sandbag  having  a  tetrahedron-shape  with  four  corners  wherein 

said  artificial  reef  structure  is  formed  by  creating  a  talus  pile  from  a  " 

4  .:^  plurality  of  structural  units  so  that  corners  of  said  structural  units 

substantially  interlock  one  said  structural  unit  with  other  structural  units 
6  of  said  plurality.  ‘ 

10.  A  method  for  focusing  ocean  wave  energy  on  a  collection  device 
2  which  comprises: 

forming  a  plurality  of  structural  units  each  comprising  a 
4  tetrahedron-shaped  bag  filled  with  sand; 

identifying  a  location  adjacent  a  shoreline  at  which  the  ocean  wave 

6  energy  is  to  be  collected;  and 

disposing  said  plurality  of  structural  units  at  said  location  to  form  at 

8  least  one  talus  pile. 

11.  A  method  for  focusing  wave  energy  as  in  Claim  1 0,  wherein  the 
2  step  of  identifying  a  location  includes  measuring  wave  height,  wave 

period,  local  tidal  elevation  and  wave  direction. 

12.  A  method  for  focusing  ocean  wave  energy  as  in  Claim  10,  wherein 
2  the  step  of  disposing  said  plurality  of  structural  units  includes  forming  two 

parallel  bars  separated  by  a  spacing  to  enhance  wave  height  of  incoming 
4  waves  within  said  spacing  at  a  location  of  the  collection  device. 

13.  A  method  for  focusing  ocean  wave  energy  as  in  Claim  10,  wherein 
2  the  step  of  disposing  said  plurality  of  structural  units  includes  forming  a 

plurality  of  V-shaped  bars  to  create  a  diffraction  pattern  for  enhancing 
4  wave  height. 
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e  ABSTRACT 

8  Ths  srtificisi  reef  of  the  present  invention  is  constructed  from  e 

plurality  of  structural  units  each  of  which  is  a  tetrahedron-shaped 

10  sandbag,  which  are  arranged  in  a  talus  pile.  The  reefs  are  formed  by 
dropping  the  structural  units  in  the  desired  locations.  The  tetrahedron- 
12  ""  shape  of  the  structural  units  provides  a  self-interlocking  structure  which 
can  be  placed  en  mass  on  the  sea  floor  without  individual  alignment  of 
14  each  bag.  The  geometry  of  the  reef  structure  is  a  sloped  bar,  a  parallel 
bar,  or  a  V-shaped  bar,  or  combinations  thereof.  The  parallel  bars,  which 
1 6  run  parallel  to  the  shoreline,  trap  standing  waves  between  the  bars  to 

reinforce  the  incoming  waves  with  a  maximum  energy  focused  where  the 
18  rays  of  trapped  waves  cross  the  wave  rays  of  the  incoming  swell.  A 
multiple  V-shaped  bar  configuration,  with  the  apex  of  the  V  pointing 
20  diagonal  to  the  incoming  waves  creates  Bragg  reflection  patterns  which 
cause  crossing  wave  rays  to  form  local  concentrations  of  wave  energy. 
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April  11,  1990 


Captain  Richard  W.  Mies 

COMSUBDEVRON-12 

Box  70 

Naval  Submarine  Base  New  London 
Groton,  Connecticut  06340 


Dear  Commodore- Mies ,  ■  ,  i  „ 

Enclosed  please  find  a  copy  of  the  test  and  evaluation  plan  for 
the  gliding  submarine  tactic  which  we  discussed  with  Admiral 
Pittenger  at  the  Tactical  Oceanography  Symposium  in  Monterey.  The 
plan  includes  a  brief  discussion  of  the  concept  and  its  possible 
application,  as  well  as  my  thoughts  on  how  to  best  go  about 
attempting  to  glide  a  submarine  in  the  initial  sea  trials.  These 
thoughts  have  been  based  on  some  discussions  I  have  had  with 
Commander  Harold  Skelly,  USN  retired,  who  is  now  living  in  San 
Diego.  The  test  and  evaluation  plan  also  inciludes  some  discussion 
of  the  optimal  glide  speeds  and  crab  angles  when  the  tactic  is 
attempted  in  ocean  currents.  These  optimal  calculations  are  based 
on  a  technical  paper  of  mine  that  will  be  published  shortly  in  the 
AIAA  Journal  of  Aircraft.  I  have  included  a  copy  of  this  paper  to 
assist  in  further  evaluation  of  the  proposed  tactic  by  you  and  your 

staff.  ,  .  ,  ^ 

I  realize  that  SSN's  and  SSBN's  have  not  been  designed  for 
operation  in  a  gliding  mode.  No  doubt  the  tacptic  would  be  more 
effective  if  submarines  had  somewhat  larger  dive  planes. 

Nonetheless  I  was  impressed  by  your  willingness  to  consider  an  "off- 
the-wall”  idea  such  as  this.  Although  the  first  attempts  at  gliding 
a  submarine  may  not  be  altogether  successful,  I  hope  enough  of  a 
positive  result  is  achieved  to  encourage  further  pursuit  of  this 
idea.  I  am  willing  and  eager  to  contribute  further  to  its 
evolution,  (619)  534-6480. 


Sincerely 


Rear  Admiral  Pittenger,  OP-96 
Dr.  Steven  Ramberg,  ONR,  Code  1121 
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THE  GLIDING  SUBMARINE  TACTIC 
by  Scott  A.  Jenkins,  Ph.D. 

Scripps  Institution  of  Oceanography 

I  INTRODUCTION 

The  purpose  of  this  test  and  evaluation  plan  is  to 
explore  the  feasibility  of  low  speed  translation  between  any 
two  points  without  turning  the  screw.  It  is  believed  that  this 
might  be  possible  by  ballasting  in  a  slightly  positive  or 
negative  state  and  then  gliding  on  sail-plane  lift  in  the 
manner  of  an  airplane  or  bird.  The  notion  here  is  that 
passive  acoustic  detectability  can  be  reduced  by  trading  the 
continual  emlttance  of  hydrodynamic  noise  from  propel lor 
swirl  for  an  occasional  "burp”  of  ballast  pump  noise  as  the 
buoyancy  is  alternately  adjusted  between  positive  and 
negative  states  during  a  dolphin  style  glide  path,  see 
Figure  1.  Tactical  advantages  would  likely  arise  following 
high  speed  transits  to  some  patrol  area,  whence  low  speed 
dolphin  glides  could  be  invoked  to  quietly  slip  away  from 
suspected  contact. 

II  CONCEPT, 

Two  distinct  glide  modes  appear  to  be  possible,  a 
sinking  glide  as  shown  in  Figure  2  or  a  rising  glide  as 
shown  in  Figure  3.  A  steady  state  sinking  glide  can 
probably  be  achieved  by  ballasting  negatively  in  a  nose  down 
attitude.  Similarly  a  steady  state  rising  glide  mode  can 
probably  be  achieved  by  ballasting  positively  in  a  nose  up 
attitude.  In  either  case  the  center  of  buoyancy  must  be 


FIGURE  1:  A)  DOLPHIN  STYLE  GLIDE  PATH 
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FIGURE  2:  SINKING  GLIDE  MODE 
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moved  forward  to  some  point  in  the  neighborhood  of  the 
center  of  pressure  on  the  sail-plane  as  shown  in  Figures  2 
and  3.  The  exact  vertical  position  of  the  center  of 
buoyancy  relative  to  the  longitudinal  axis  of  the  hull  is 
not  particularly  critical  from  the  standpoint  of  gliding/ 
dynamics. 

Naturally,  a  submarine  ballasted  nose  heavy  as 
diagrammed  in  Figure  2  will  begin  to  sink  along  a  trajectory 
that  will  continually  steepen  to  a  nearly  vertical  nose  down 
descent.  If,  however,  the  sail-planes  are  adjusted  to  a 
positive  angle  of  incidence,  >  0,  then  they  will  generate 
a  lift  force,  L,  acting  at  right  angles  to  the  descent 
trajectory  as  shown  in  Figure  2.  The  sail-plane  lift  will 
prevent  the  trajectory  from  becoming  vertical.  Instead,  the 
trajectory  will  assume  some  glide  slope,  9,  along  which  the 
vertical  component  of  the  lift  force,  Lcos  9,  balances  the 
gravitational  force,  Fg,  due  to  the  negative  buoyancy.  The 
submarine  will  glide  at  some  speed,  u,  along  this  glide 
slope  for  which  the  drag  on  the  submarine,  D,  exactly 
balances  th^/^ravitational  thrust,  T,  due  to  the  component 
of  negative  buoyancy  acting  along  the  glide  slope,  Fgsin  9. 

A  nose-down  pitching  moment  will  result  from  the  dynamic 
pressure  distribution  associated  with  the  lift  of  the  sail¬ 
planes.  This  nose-down  pitching  moment  must  be  corrected  by 
a  slight  downward  trim  force,  S,  exerted  by  a  negative  angle 
of  incidence,  7,  on  the  stern-planes.  The  condition  of  such 
a  force  balance  is  known  as  a  steady  state  glide.  The 


gravitational  force  due  to  the  negative  buoyancy  is  pulling 
the  submarine  down  along  its  glide  slope  and  causing  it  to 
translate  forward  in  the  process  at  a  speed,  ucos  $. 

A  similar  steady  state  glide  is  possible  when  the 
submarine  is  ballasted  positively  in  a  nose  up  attitude  as 
shown  in  Figure  3.  In  this  case  the  sail-planes  are  set  at 
a  negative  angle  of  incidence,  ^  <0,  giving  rise  to  a  lift 
force  acting  at  right  angles  to  the  ascent  trajectory  in  the 
downward  direction.  This  downward  oriented  lift  force 
prevents  the  ascent  trajectory  from  becoming  vertical  and 
establishes  a  steady  state  ascending  glide  slope  along  which 
the  vertically  downward  component  of  lift,  Lcos  tf,  balances 
the  positive  buoyancy,  Fb*  An  upward  pitching  moment  due  to 
the  dynamic  pressures  on  the  sail-planes  must  be  corrected 
by  a  slight  upward  trim  force,  S,  resulting  from  a  positive 
angle  of  incidence,  7,  on  the  stern-planes.  When  trimmed  in 
this  manner,  the  submarine  ascends  along  a  constant  glide 
slope  at  a  glide  speed,  u,  such  that  the  drag  on  the 
submarine  is  balanced  by  the  component  of  buoyancy  acting 
upward  alorigy|-he  glide  slope,  Fjjsin  6.  While  the  submarine 
ascends  along  this  rising  glide  slope,  it  will  translate 
forward  at  a  speed,  ucos  e. 

While  translating  forward  in  either  the  sinking  or 
rising  glide  mode,  the  submarine  must  execute  a  certain 
depth  excursion.  Because  depth  excursions  must  remain 
limited  to  take  advantage  of  sound  ducts,  a  succession  of 
alternating  rising  and  sinking  glides  Will  be  necessary  in 


order  to  translate  over  significant  distances.  This  will 
result  in  a  dolphin  style  glide  path  as  diagrammed  in  Figure 
1,  during  which  it  will  be  necessary  to  periodically  operate 
the  ballast  pumps  to  reconfigure  the  buoyancy  and  trim  for 
alternating  rising  and  sinking  glide  modes.  If  a 
controllable  pitch  screw  is  available,  it  should  be 
feathered  to  reduce  drag  and  hydrodynamic  noise  during  each 
glide  mode. 

Ill  GLIDE  POLARS 

Consider  the  sinking  glide  mode  in  Figure  2.  For  any 
giv4n  degree  of  negative  buoyancy,  Fgi,  Fg2/  Fg3,...,  there 
will  be  a  unique  relation  between  the  sail-plane  angle  or 
sinking  speed  and  the  maximum  obtainable  steady  state  glide 
speed.  That  relation  is  known  as  a  glide  polar.  Each 
degree  of  negative  buoyancy  will  have  its  own  specific  glide 
polar.  For  any  given  ballast  configuration,  smaller  sail¬ 
plane  angles  will  result  in  higher  gliding  speeds,  but  at  a 
cost  of  greater  sinking  speeds  or  greater  depth  excursions. 

Conversely  the  smallest  sinking  speeds  or  depth  excursions 

1 

result  froni.the  largest  sail-plane  angles  but  at  a  cost  of  a 
slower  glide  speed.  The  smallest  depth  excursions  and 
sinking  speeds  result  from  the  largest  sail-plane  angle  that 
can  be  achieved  without  the  sail-planes  stalling.  Sail¬ 
plane  stall  therefore  limits  the  glide  polar  at  some  minimum 
gliding  speed  and  maximum  sail-plane  angle  setting.  Greater 
degrees  of  negative  buoyancy  result  in  higher  glide  speeds 
for  any  given  sail-plane  angle.  However,  the  higher  glide 


speeds  resulting  from  greater  negative  buoyancy  are  achieved 
at  a  cost  of  higher  sinking  speeds 'or  greater  depth 
excursions. 

The  rising  glide  mode,  Figure  3,  also  has  a  family  of 
glide  polars  associated  with  it  depending  upon  the  degree,  of 
positive  buoyancy.  In  this  case  however  larger  (less 
negative)  sail-plane  angles  result  in  higher  gliding  speeds. 
The  smallest  rising  speeds  or  depth  excursions  result  from 
the  smallest  (most  negative)  sail-plane  angles  that  can  be 
achieved  before  sail-plane  stall  ensues.  For  any  given 
sail-plane  angle,  the  glide  speed  will  be  increased  by 
increasing  the  degree  of  positive  buoyancy,  but  at  a  cost  of 
greater  rising  speeds  or  depth  excursions. 

Despite  the  wealth  of  design  data  on  SSN's  and  SSBN's, 
the  glide  polars  can  probably  hot  be  calculated  exactly.  In 
all  likelihood  they  will  have  to  be  measured  directly  in  sea 
trials.  This  is  due  to  the  fact  that  the  glide  polar 
results  from  the  complex  variations  of  lift  and  drag  with 
glide  speed.  Nonetheless  the  glide  polar  can  probably  be 
represented.^^  a  parabola  of  some  sort  once  the  data  is 
obtained  because  a  large  portion  of  the  drag  is  due  to 
profile  drag  on  the  hull  and  induced  drag  on  the  dive  planes 
whose  general  behavior  with  glide  speed  is  known.  Therefore 
it  will  not  be  necessary  to  measure  a  glide  polar  for  each 
particular  degree  of  positive  or  negative  buoyancy.  The 
glide  polar  measured  for  a  particular  ballast  configuration 
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can  be  used  to  calculate  the  glide  polars  for  all  other 

possible  ballast  configurations. 

IV  GLIDE  SLOPE  MINIMIZATION 

In  an  operational  scenario  one  seeks  to  ininiinize  the 
numbers  of  times  that  the  ballast  pumps  must  be  operated  / 
during  a  dolphin  style  glide  path  between  any  two  points. 
Therefore,  it  is  desirable  to  operate  at  a  combination  of 
dive  plane  angles  and  glide  speeds  on  the  glide  polar  which 
minimizes  the  glide  slope  in  either  the  rising  or  sinking 
mode.  These  optimal  points  are  shown  in  Figures  2  and  3  as 
the  optimal  glide  tangent  lines  which  intersect  the  glide 
polars  for  particular  degrees  of  positive  and  negative 
buoyancy.  These  optimal  tangent  lines  are  shown  for  the 
case  of  an  ocean  with  no  motion.  However  in  a  real  ocean 
with  currents,  the  glide  polars  become  distorted  depending 
on  whether  the  current  is  acting  in  the  direction  of  the 
glide  path  or  at  some  angle,  a,  to  it.  For  the  general  case 
of  rising  or  sinking  glides  at  some  arbitrary  angle,  a,  to 
the  current,  there  will  be  an  optimal  combination  of  glide 
speed  and  grab  angle  which  minimizes  the  glide  slope,  rising 
speed,  and  depth  excursion  while  maintaining  the  desired 
course  heading.  Some  hypothetical  calculations  of  the 
optimal  glide  speeds  and  crab  angles  are  shown  in  Figures  2 
and  3.  The  curves  are  calculated  for  an  assumed  parabolic 
glide  polar.  The  algorithms  are  rather  complicated  and  are 
based  on  reference  1.  However,  the  qualitative  features  of 
these  calculations  are  perfectly  general.  For  instance,  the 
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optimal  glide  speed  directly  into  a  current  (a  =  0)  is 
faster  than  for  the  case  of  a  glide  with  the  current.  In 
order  to  minimize  the  glide  slope  and  vertical  depth 
excursions,  an  11  knot  glide  is  required  directly  against  a 
5  knot  current  whereas  only  a  2.5  knot  glide  is  required  for 
gliding  along  with  ^5  knot  current.  If  the  glides  are 
performed  at  a  =90°  to  the  currents,  then  faster  glide 
speeds  are  required  for  higher  currents  in  order  to  minimize 
the  depth  excursions.  For  glides  performed  at  o  =  45° 
either  against  or  with  the  current,  a  sufficiently  fast 
glide  speed  is  required  such  that  the  crab  angle  does  not 
have  to  exceed  45°  in  order  to  maintain  the  required  course 
heading.  Otherwise  the  submarine  is  not  making  adequate 
progress  towards  its  goal.  There  are  additional  corrections 
to  the  optimal  glide  speeds  and  crab  angles  when  vertical 
water  motions  are  encountered  such  as  along  fronts,  eddies 
and  internal  waves.  These  can  also  be  calculated  as 
described  in  reference  1. 

The  optimal  glide  speeds  (speed-to-glide)  and  the 
optimal  crdb%hgies  for  varying  current  regimes  are  all 
based  on  the  rising  and  sinking  glide  polars  in  still  water. 
In  other  words,  once  the  glide  polar  is  known,  it  is 
possible  to  make  all  of  the  necessary  subsequent 
calculations  for  glides  in  a  moving  ocean.  These 
calculations  can  be  easily  summarized  in  either  nomographic 
tables  or  on  graphs  such  as  those  shown  in  Figures  2  and  3. 
Adherence  to  these  optimal  combinations  of  glide  speeds  and 


crab  angles  is  necessary  to  minimize  the  numbers  of  times 
that  the  ballast  pumps  must  be  operated  when  translating 
between  any  two  points  along  a  dolphin  style  glide  path. 

V  ACHIEVING  STEADY  STATE  GLIDE 

The  first  attempt  should  probably  be  made  by  coasting 
into  glide  after  setting  the  buoyancy  and  ballast  trim  while 
under  power  and  then  stopping  the  screw.  This  will  ensure 
sufficient  control  authority  while  making  the  transition 
from  power,  to  glide.  If  on  the  other  hand  the  attempt  is 
made  fromr  a  static  condition,  then  the  submarine  must  either 
fall  or  rise  from  a  state  of  rest  and  transition  through  a 
low  speed  acceleration  during  which  control  authority  could 
be  marginal.  Furthermore  the  first  attempts 'to  transition 
from  power  to  glide  should  be  made  for  the  rising  glide  mode 
in  the  event  of  unforeseen  difficulties. 

The  following  sequence  of  events  for  the  first  glide 
attempts  are  suggested: 

a)  Establish  a  constant  speed  under  power  of  5  knots; 

b)  Adjust  the  ballast  while  under  power  for  positive 
buoyartcj^i  and  nose  up  trim.  (The  center  of  buoyancy 
must  be  moved  at  least  as  far  forward  as  the  mean 
hydrodynamic  center  of  pressure  on  the  sail-planes. 

The  point  at  which  the  sail-planes  pivot  is 
sufficiently  far  forward  for  this  condition.); 

c)  Adjust  the  sail-planes  for  negative  angle  of 
incidence  while  the  ballast  is  being  adjusted  for 
positive  buoyancy.  Establish  a  negative  sail-plane 


angls  and  positive  stern  plane  angle  of  incidence 
sufficient  to  maintain  constant  depth  once  ballasting 
for  positive  buoyancy  and  a  forward  center  of  buoyancy 
has  been  achieved. 

d)  Stop  and  feather  the  screw.  Allow  the  submarine, to 
decelerate  to  a  constant  steady  state  glide  speed. 
Adjust,  the  stern-plane  trim  to  some  positive  angle  of 
Incidence  sufficient  to  maintain  a  constant  nose  up 
attitude  along  the  glide  slope.  The  idea  here  is  to 
keep  the  nose  pointed  directly  into  the  flow,  and  to 
avoid  flying  the  hull  at  some  angle  to  the  flow. 

Record  the  sail-plane  angle,  the  rate  of  change  of 
depth,  and  the  glide  speed. 

e)  Increase  the  negative  incidence  of  the  sail-planes 
and  the  positive  incidence  of  the  stern-planes  by  a 
small  amount.  Allow  the  submarine  to  decelerate  to  a 
new  steady  state  glide  speed.  Record  the  sail-plane 
angle,  the  rate  of  change  of  depth,  and  the  glide 
speed. 

vf  ’  ■  ,  • 

f)  Report  steps  d  and  e  in  small  increments  of 
increasingly  negative  sail-plane  angle  and  positive 
stern  plane  angle  until  sail-plane  stall  ensues. 

g)  Return  to  the  original  5  knot  powered 
configuration.  Adjust  ballast  and  trim  to  the  same 
settings  as  above.  Once  the  submarine  has  decelerated 
to  a  steady  state  glide,  decrease  the  degree  of 
negative  angle  of  incidence  on  the  sail-planes.  Retrim 


the  stern  planes  to  fly  the  hull  directly  into  the  flow 
along  the  glide  slope.  Repeat  steps  e-f  for  successive 
reductions  in  sail-plane  angle  until  unreasonably  large 
ascending  glide  slopes  (greater  than  30  or  45°)  result, 
h)  The  above  procedure  will  essentially  detail  tho  y 
glide  polar.  If  the  experiment  is  performed  with 
little  or  no  current,  then  the  data  will  be  sufficient 
for  following  optimization  calculations. 

Once  steady  state  rising  glide  modes  have  been 
demonstrated,  then  the  above  procedures  can  be  followed  for 
ballasting  for  negative  buoyancy  and  nose  down  trim.  It 
will  be  necessary  to  repeat  the  above  steps  for  sinking 
glide  modes  to  determine  the  sinking  glide  polar.  The 
sinking  glide  polar  will  not  necessarily  be  the  exact 
inverse  of  the  rising  glide  polar,  since  the  sail-planes 
atop  the  sail  do  not  produce  the  same  interference  drag 
effects  with  the  hull  for  both  positive  and  negative  angles 
of  incidence. 

VI  DOLPHIN  GLIDES 

Transitions  from  power  to  glide  while  underway  should 
be  repeated  a  sufficient  number  of. times  to  develop  crew 
proficiency.  Once  such  proficiency  is  achieved,  then 
transitions  from  rising  glide  modes  to  sinking  glide  modes 
can  be  attempted  in  order  to  demonstrate  dolphin  style 
gliding  over  large  distances.  Initially  these  transitions 
should  be  attempted  at  fairly  high  glide  speeds  and  steep 
glide  slopes,  in  order  to  insure  substantial  control 


authority.  Once  dolphin  gliding  has  been  proven  feasible  at 
relatively  inefficient  high  glide  speeds,  then  the  maneuver 
can  be  practiced  at  lower  glide  speeds  and  flatter  glide 
slopes  for  which  the  tactic  will  be  operationally  most  _ 
effective.  Once  such  dolphin  gliding  proficiency  has  been 
demonstrated,  then  detailed  studies  of  the  acoustic 
signature  during  long  distance  dolphin  glides  can  be 
undertaken.  Such  tests  should  be  undertaken  only  for  a 
pitch  controllable  screw.  The  drag  and  acoustic  noise  while 
gliding  with  an  unfeathered  screw  would  likely  prove  to  be 
prohibitive. 
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Theory  and  experiment  are  presented  for  optimal  combinations  of  glide  speeds  and  crab  angles  that  minimize 
the  glide  slope  at  any  arbitrary  angle  to  the  wind  along  a  constant  course  heading.  The  optimization  scheme  is 
formulated  for  constant  wind  speeds,  wind  direction,  and  air  mass  sink  rates.  An  analytic  solution  is  found  in 
the  asymptotic  limit  of  a  small  crosswind  component.  A  general  solution  is  obtained  by  seeded  iterations  with 
a  Taylor  series  expansion  about  that  limit.  The  solution  was  verified  in  test  flights  with  quartering  tailwinds  and 
direct  crosswinds  at  the  700,  500,  and  300  mb  levels.  It  was  concluded  that  relatively  small  pilot-induced  speed 
errors  result  in  significant  glide  slope  degradations  with  unnecessary  and  potentially  dangerous  altitude  losses. 
Applications  of  the  results  are  considered  for  wind  fields  that  vary  slowly  in  space,  such  as  lee  waves  and  Rossby 
waves. 


Nomenclature 


Ao 

=  wing  area,  ft^ 

A,  B, 

A',  B',  G 

=  flight  path  way  points 

a,  b,  c 

-  glide  polar  coefficients  at  altitude 

a 3  6,  c 

=  glide  polar  coefficients  at  sea  level 

Co 

=  quadratic  drag  coefficient 

Co.i 

=  quadratic  induced  drag  coefficient 

Cd,p 

=  quadratic  profile  drag  coefficient 

Cl 

=  quadratic  lift  coefficient 

Jny  fn 

=  resolvent  cubic  parameters  for  iterative 

^o»/o 

general  solution 

=  resolvent  cubic  parameters  for  weak  crosswinds 

K^ 

=  profile  drag  factor 

Ki 

=  induced  drag  factor 

ATa 

=  glide  polar  factor  due  to  profile  drag 

K, 

=  glide  polar  factor  due  to  induced  drag 

L/D 

=  inverse  glide  slope 

L/Djt^^ 

=  inverse  glide  slope  due  to  speed  to  fly 

L/D(u  = 
U-^bU) 

=  inverse  glide  slope  due  to  speed-to-fly  errors 

1 

=  wing  cord,  n.mi. 

m 

=  gross  weight,  lb 

n 

=  integral  index  number 

Pn 

=  second-order  coefficient  for  general  iterative  case 

Po 

=  second-order  coefficient  for  weak  cross  winds 

Qn 

=  first-order  coefficient  for  general  iterative  case 

Qo 

=  first-order  coefficient  for  weak  cross  winds 

R 

=  Rossby  radius  of  deformation,  n.mi. 

rn 

=  zero-order  coefficient  for  general  iterative  case 

ro 

=  zero-order  coefficient  for  weak  crosswinds 

S 

=  air  mass  sink  rate,  kt 

u 

=  speed-to-fly  at  altitude,  kt 

0 

=  speed-to-fly  at  sea  level  in  still  air,  kt 

Ui 

=  indicated  speed-to-fly,  kt 

Un 

=  general  speed-to-fly  for  nth  wind  increment,  kt 

Uo 

=  speed  to  fly  in  weak  crosswinds,  kt 
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u  =  horizontal  component  of  glide  velocity  relative 

to  air  at  altitude,  kt 

u  =  horizontal  component  of  glide  velocity 

relative  to  air  at  sea  level,  kt 
Us  =  stall  speed,  kt 

V  =  wind  speed,  kt 

Vo  =  weak  wind  speed,  kt 

w  =  vertical  component  of  glide  velocity  relative  to 

air  at  altitude,  relative  sinking  speed  in  knots 
w  =  vertical  component  of  glide  velocity  relative  to 

air  at  sea  level,  relative  sinking  speed  in  knots 
X  =  position  along  flight  path  relative  to  any  given 

way  point 

a.  =  wind  angle  relative  to  course  line 
=  crab  angle  relative  to  course  line 
7  =  profile  drag  exponent 

bU  =  speed-to- fly  error,  kt 

bv  =  wind  speed  increment,  kt 

7}  =  relative  squared  crosswind  velocity 

difference,  kt^ 

6  =  minimum  glide  slope  in  still  air  at  sea  level 

X  =  wave  length,  n.mi. 

p  =  density  at  altitude,  slug/ft^ 

p  =  density  at  sea  level,  slug/  ft^ 

Introduction 

HE  first  treatment  of  the  problem  of  the  optimal  glide 
speed  that  minimizes  the  glide  slope  in  a  moving  atmo¬ 
sphere  was  due  to  MacCready.**^  The  solutions  have  since 
been  referred  to  as  “speed-to- fly.”  MacCready’s  original 
work  was  valid  only  for  a  convective  atmosphere  in  which 
regions  of  rising  or  sinking  air  translate  horizontally  with  the 
wind.  His  speed-to-fly  solutions  were  thus  independent  of  the 
wind  speed  and  determined  uniquely  by  the  rate  of  vertical  air 
mass  movement  and  the  glide  polar  in  still  air  at  sea  level. 

Subsequently,  Kuettner,^  following  arguments  introduced 
by  Reichmann,^  considered  the  effects  on  glide  slope  mini¬ 
mization  when  regions  of  rising  or  sinking  air  do  not  translate 
with  the  wind  over  the  ground.  Such  a  condition  is  typically 
encountered  in  both  hydrostatic  and  nonhydrostatic  lee  waves 
(see  Ref.  5).  Kuettner  considered  only  the  cases  of  glides 
directly  with  the  wind.  The  solutions  for  speed-to-fly  by  both 
Kuettner  and  Reichmann  were  graphical  in  nature,  involving 
tangent  plotting  for  a  few  selected  cases  with  no  matherriatics 
as  such  presented.  Even  so,  some  striking  departures  from 
MacCready’s  original  findings  were  uncovered  by  both  au¬ 
thors.  Speed-to-fly  was  found  to  be  dependent  on  the  wind 
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speed  with  the  flattest  glide  slopes  achieved  for  tailwinds  at 
speeds  less  than  the  MacCready  speed-to-fly.  Kuettner  also 
considered  the  degradation  in  glide  slope  at  high  altitudes  with 
a  direct  tailwind.  He  concluded  that  the  practice  of  ballasting 
to  high  wing  loadings  is  disadvantageous  since  the  reduced  air 
density  at  high  altitudes  has  already  effectively  raised  the  wing 
loading  considerably.  The  increase  in  relative  sinking  speeds  at 
these  higher  wing  loadings  cannot  be  compensated  in  down¬ 
wind  glides  because  slow  speeds-to-fly  are  required. 

In  the  present  paper,  the  speed-to-fly  problem  is  treated  for 
the  general  case  of  glides  along  constant  course  headings  at  an 
arbitrary  angle  to  the  wind  for  any  given  altitude.  Optimal 
flight  paths  for  crossing  lee  waves  are  shown  to  result  from  an 
orthogonal  series  of  constant  course  glides.  The  formulation 
of  the  problem  is  based  on  a  quadratic  glide  polar  in  a  con¬ 
stant  wind  field.  The  constant  wind-field  assumption  provides 
acceptable  accuracy  for  treating  glides  in  lee  waves  and  Rossby 
waves,  which  vary  over  geophysical  length  scales.  The  result¬ 
ing  speed-to-fly  equation  is,  nonetheless,  found  to  be  tran¬ 
scendental.  An  asymptotic  analytic  solution  is  found  in  the 
limit  of  small  crosswind  components.  A  general  solution  for 
arbitrary  crosswind  components  is  then  developed  based  on  a 
seeded  iteration  with  a  Taylor  series  expansion  about  the  weak 
crosswind  solution.  The  results  are  compared  with  data  from 
skewed  test  glides  performed  between  the  700-  and  300-mb 
levels.  These  comparisons  show  that  speed-to-fly  errors  can 
degrade  the  glide  slope  to  such  an  extent  that  the  performance 
gains  achieved  by  modern  airfoils  and  composite  construction 
are  lost. 

Formulation 

Consider  glides  performed  at  l-g  load  factor  across  a  field 
of  lee  waves  in  such  a  manner  that  a  constant  course  heading 
is  maintained  between  any  two  arbitrary  points  (A,B)  in  an 
Eulerian  frame.  Let  the  course  heading  connecting  (A,B)  be  at 
an  angle  a  to  the  wind  whose  speed  is  v  as  shown  in  Fig.  1 .  If 
the  horizontal  component  of  the  glide  velocity  relative  to  the 
air  is  w,  then  it  must  be  directed  at  some  crab  angle  jS  in  order 
to  maintain  a  constant  course  heading  between  (A,B).  The 
crab  angle  is  thus  given  by: 

jS  =  sin“^ [  —  (v/w)  sina]  (1) 

While  traversing  the  lee  wave  field,  regions  of  rising  flow 
(lift)  will  be  encountered  along  the  upwind  faces  of  the  lee 
wave  crests,  and  descending  flow  (sink)  along  the  downwind 
faces.  Let  the  vertical  velocity  component  of  the  air  mass  be  S, 
which  shall  be  taken  as  positive  in  the  downward  direction  in 
an  Eulerian  frame.  If  the  vertical  component  of  the  glide 
velocity  relative  to  the  air  is  w  (positive  downward),  then  the 
inverse  glide  slope  (L/D)  with  respect  to  an  Eulerian  frame 
such  as  the  ground  will  be 


j-LEE  WAVE  FIELD  H 

Fig.  1  Schematic  of  constant  course  strategies  with  lee  waves. 


Any  optimal  strategy  for  crossing  a  field  of  lee  waves  is 
given  by  one  or  more  glides  along  constant  course  headings. 
This  is  because  the  lee  wave  field  structures  the  air  mass  sink 
rates  in  parallel  bands  having  high  cross-stream  coherence. 
Within  the  lee  wave  field,  an  orthogonal  series  of  constant 
course  glides  yields  the  optimum  flight  path  between  any  two 
points  (A,B)  that  are  separated  by  some  streamwise  excursion 
(see  Fig.  1).  It  is  clear  from  inspection  of  Eq.  (2)  that  straight 
line  glides  with  a  direct  tailwind,  a  =  0,  give  the  maximum 
possible  L/D  when  entering  the  lee  waves  along  (A,B')  or 
when  traversing  the  lee  wave  field  along  (A',B).  These  two 
glides  can  be  joined  by  a  direct  crosswind  glide  at  constant 
course  heading,  (B%A'),  utilizing  lee  wave  lift  to  achieve  an 

infinite  L/D  when  5 - w.  However,  sufficiently  long  cross- 

wind  glides  in  uninterrupted  lift  such  as  (B',A')  are  not  always 
possible.  Variations  in  the  alignment  and  width  of  wave  gener¬ 
ating  mountains  often  destroy  the  necessary  cross-stream  co¬ 
herence  of  the  lee  wave  field.  If  the  ultimate  goal  is  signifi¬ 
cantly  far  downwind,  say  point  G  in  Fig.  1,  then  skewed, 
constant  course  glides  such  as  (B,  G)  will  become  an  inevitable 
part  of  cross-country  strategy. 

Before  proceeding  to  optimize  any  of  the  constant  course 
glides  set  forth  in  Fig.  1 ,  it  is  desirable  to  eliminate  w  from  Eq. 
(2)  in  order  to  save  labor.  The  vertical  and  horizontal  compo¬ 
nents  of  the  glide  velocity  relative  to  the  air  mass  are  related  by 
the  glide  polar  at  sea  level  in  still  air.  This  relation  can  be 
posed  from  the  usual  assumption  that  the  total  drag  is  a  linear 
combination  of  profile  drag  and  induced  drag: 

Cd  =  Cd.p  +  Cdj  (3) 

Because  of  Reynolds  number  dependence,  the  profile  drag 
coefficient  is  not  constant  but  will  vary  with  the  glide  velocity 
in  still  air  at  sea  level  U  according  to 

(4) 

where  Ki  is  a  factor  that  varies  with  the  kinematic  viscosity 
and  characteristic  length  scale  and  y  will  depend  on  the  rela¬ 
tive  proportion  of  the  total  wetted  surface  area  that  is  sub¬ 
jected  to  laminar  vs  turbulent  boundary  layers.  The  induced 
drag  coefficient  on  the  other  hand  may  be  inferred  from  lifting 
line  theory^’*^  as 

Cdj=K2CI  (5) 

where  the  factor  Ki  is  a  function  of  wing  aspect  ratio,  twist, 
taper,  and  aeroelasticity. 

Based  upon  a  drag  formulation  by  Eqs.  (3~5),  the  glide 
polar  at  sea  level  may  be  written  as 


w  =  (Cfl/Ci)  a  =  +  (K^/u) 

(6a) 

Kj  =  (fiAoKi)/(2m) 

(6b) 

Kt  =  {2mK2)/(pAo) 

(6c) 

If  we  expand  Eqs.  (6)  in  a  Taylor  series  about  some  optimal 
speed-to-fly  in  still  air,  u-U,  then  the  glide  polar  is  quadratic 
to  second  order  according  to 

+  {Ki/u)-~ciU^  +  5a  +  & 

(7) 

where 

a  =  (3-i)(2-y)K30'-y  +  aKA/C^) 

6  =  C3-y)Ki0^-y-(K,/0)-^20di 

+  (K4/0)  +  W^-60 
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Fig.  2  Quadratic  approximation  to  glide  polar  data  at  sea  level  and 
at  altitude. 


Here  a,  6,  and  c  are  based  on  the  air  density  at  sea  level.  We 
may  then  correct  the  quadratic  glide  polar  in  Eq.  (7)  to  any 
altitude  where  the  air  density  is  p  by  taking 


w  =  au^  +  bu  +  c 

(8a) 

a  =  a  (p/p)-''^ 

(8b) 

11 

(8c) 

c  =  c 

OO 

0(/?/),  where  n  --  5  according  to  Ref.  9.  Therefore,  the  relative 
sizes  of  the  terms  in  Eq.  (10)  are: 

=  (lla) 

du  du  dx  \nl/ 

=  =  (lib) 

dv  dv  dx  \R  / 

d(L/D)  ^  d(L/D)doi^  Jdx\ 

3(L/D)  „  _  aiL/D)  as 


In  the  atmosphere,  R  -  0(100  n.mi.),  whereas  X  ~  0(10 
n.mi.)  for  hydrostatic  lee  waves  or  X  0(1  n.mi.)  for  non¬ 
hydrostatic  lee  waves.  These  are  immense  compared  to  the 
wing  cord.  Therefore,  changes  in  L/D  with  respect  to  wind 
speed,  direction,  and  air  mass  sink  rate  are  negligible  com¬ 
pared  to  those  with  respect  to  changes  in  glide  speed  over  any 
incremental  distance  dx.  Hence,  the  speed-to-fly  U,  which 
maximizes  the  L/D  over  the  ground,  is  given  to  0(2xrt//X) 
accuracy  by 

d  L 

—  ^  =  0atw  =  l/  (12) 

du  D 


Equation  (12)  yields  the  general  speed-to-fly  equation  for 
glides  skewed  relative  to  the  wind,  which  may  be  written  as 

-  aU^  -  lavU  cosWL/^  -  sin^a 


A  comparison  is  shown  in  Fig.  2  between  the  quadratic 
approximation  according  to  Eq.  (7)  and  measured  glide  polar 
data  for  the  test  aircraft,  a  Schempp-Hirth  Nimbus  Ilb,  taken 
from  Johnson®  and  corrected  to  sea  level.  The  solid  curve  is  a 
least-squares  best  fit  of  a  quadratic  to  the  data.  The  values  for 
the  sea  level  polar  coefficients  that  result  from  this  fit  are 


a  =  0.0012155  kt-' 

(9a) 

6=  -0.1106912 

(9b) 

a  =  3.564157  kt 

(9c) 

-f*  (c  +  S  +  2av'^  sin^a)  U  —  vb 

+  v'^b  sin^o:  =  0  (13) 

The  transcendental  nature  of  the  general  speed-to-fly  equation 
[Eq.  (13)]  follows  from  the  fact  that  speed-to-fly  is  a  function 
of  the  crab  angle  by  way  of  Eqs.  (2)  and  (12),  and  crab  angle 
is,  in  turn,  a  function  of  speed-to-fly  as  required  by  Eq.  (1). 
Analytic  solutions  to  Eq.  (13)  arise  in  the  asymptotic  limit  of 
an  indefinitely  small  crosswind  component,  sin^a,  for 

which  the  speed-to-fly  w  =  i/o  is  given  by 


Corrected  glide  polar s  for  the  700- ,  500- ,  and  300-mb  altitudes 
are  also  shown  in  Fig.  2  as  dashed  and  dotted  curves.  Note 
that  the  maximum  inverse  glide  slope  of  47.7  remains  un¬ 
changed  with  altitude  for  still  air.  The  deformation  of  the 
polar  at  altitude  is  equivalent  to  increasing  the  wing  loading  by 
factors  of  1.20,  1.42,  and  1.84  for  the  700-,  500-,  and  300-mb 
levels,  respectively. 

With  Eq.  (2)  expressed  in  terms  of  Eqs.  (8),  the  inverse  glide 
slope  is  reduced  to  dependence  on  four  variables,  u,  v,  a,  5. 
Consequently,  incremental  changes  in  L/D  are  expressed  as 


T/o'+Pot/o'  +  ^ot/o  +  /-o  =  0  (14) 

where 

Po  =  2vo  cosa 

^0  =  [  -  (c  +  S  +  lavl  sin^Qf  -  Vq^  cosa)/fl] 

To  =  (Vq^  sinW^z) 

Vo  =  V  V (/o/sin^a 


d(L/D)  = 


d{L/D) 

du 


du  -1 - :: - 


dv  -I- 


djL/D) 

da 


da 


The  speed-to-fly  equation  thus  collapses  to  an  ordinary  cubic 
equation  for  weak  crosswinds,  yielding  the  following  solution: 


d{L/D) 
^  dS 


dS 


(10) 


Now,  each  of  the  variables  in  Eq.  (10)  changes  with  incremen¬ 
tal  changes  in  position  along  the  flight  path  dx  according  to 
characteristic  length  scales.  The  air  mass  sink  rates  vary  over 
distances  comparable  to  the  wave  length,  0(X/27r),  of  the  lee 
waves.  The  wind  speed  and  direction  vary  over  the  lowest 
mode  length  scales  of  the  upper  level  storms,  typically  the 
Rossby  radius  of  deformation,  0(R),  as  discussed  in  Ref.  5. 
On  the  other  hand,  the  horizontal  component  of  the  glide 
velocity  can  vary  over  distances  of  only  a  few  cord  lengths. 


t/o  = 


y  A  2l\  L  2  A/  4  27_ 


1/3 


2vo  cosa 


(15) 


where 

=  {\/21){2pl  ^  9poQo  +  27ro) 
/o  =  (l/3)(3go-Po') 
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The  remaining  two  roots  of  Eq.  (14)  are  less  than  zero  (back¬ 
wards  flight)  and  are  therefore  not  sensible.  As  a  particular 
case  of  Eq.  (15),  the  speed-to-fly  solution  in  a  direct  headwind 
or  tailwind  (a  =  0),  which  corresponds  to  the  graphical  result 
by  tangent  plotting  due  to  Kuettner,^  is  found  to  be 

Uo(a  =  0)  =  (l/2a)[4a^v^  +  4a(c  +  S  -  bv)Y^^  -  v  (16) 

For  a  general  solution  to  Eq.  (13),  a  seeded  iteration  is 
begun  from  the  weak  crosswind  solution  Uq  for  some  small 
wind  speed  v  =  Vq.  With  each  iterative  step  n  thereafter,  the 
wind  speed  is  stepped  incrementally  by  ndv.  The  numerical 
computations  for  any  given  step  n  are  based  on  a  Taylor  series 
expansion  of  (t/^  “  sin^a)*^^  about  the  speed-to-fly  solution 

from  the  previous  iteration,  =  By  this  numerical 

scheme,  the  general  speed-to-fly  equation  for  the  nth  iteration 
is 

(17a) 

where 


metric  function  about  zero  wind  speed  and  is  slower  than  all 
headwind  cases  and  faster  than  all  tailwind  cases.  All  remain¬ 
ing  possible  crosswind  solutions  are  intermediate  between 
these  two  extremes,  with  small  differences  among  skewed 
glides  into  the  wind,  but  significantly  faster  speeds-to-fly 
for  increasing  crosswind  components  with  the  wind.  When 
0<q!<90  deg,  the  maximum  inverse  glide  slopes  apparently 
result  from  keeping  the  crab  angle  small,  /3<0(45  deg). 
Hence,  the  speeds-to-fly  during  skewed  glides  with  the  wind 
increase  with  increasing  crosswind  component  in  order  to 
maintain  the  course  line  without  resorting  to  exceedingly  large 
crab  angles. 

It  is  also  interesting  to  note  from  Figs.  3  and  4  that  the 
indicated  speed  to  fly  (based  on  pitot-static  airspeed  systems) 
varies  with  altitude  for  any  given  wind  speed  and  direction. 
This  effect  was  not  considered  by  MacCready.*’^  Generally, 
slower  indicated  airspeeds  are  required  at  higher  altitudes 
for  any  nonzero  wind  speed  regardless  of  direction.  Conse¬ 
quently,  the  indicated  speeds-to-fly  and  crab  angles  for  60-kt 
winds  at  700-mb  are  comparable  to  those  required  for  90-kt 
winds  at  the  300-mb  level.  The  true  speeds-to-fly  increase  with 


1  -  vb  cosa[(l/2)7?-^^^  -  -  2av  cosajU^^iV  ^^^)1 

{  -  a  -  av  cosairij-^^^  - 

{c  S  +  2av^  sin^o!  —  vb  cosaU^^  ”  2gv  cosah^^^  —  (1/2)C/^^  ~  (l/2)Un- I 

[  -ff-JV  COSq:(i7“^^^  - 


1  v^b  sin^g  -  vb  cosaW'^-  -  (1/2)U„^  -  (l/2)f/„^  ) 

~  I  -  a  -  av  cosa(r}~^^^  - 


(17d) 


where  n  =  1,2,3...;  rj  =  —  v^sin^a;  v  =  Vq  +  ndv.  Hence, 

the  general  speed-to-fly  solution  for  any  given  wind  speed  and 
arbitrary  a  is  given  by 


increasing  altitudes  but  by  less  than  a  constant  factor  (p/p)^^^. 
This  results  from  the  fact  that  the  higher  effective  wing  load¬ 
ing  at  altitude  enables  more  efficient  penetration  into  a  head¬ 
wind  component  but  requires  slower  glides  at  a  higher  Ci  with 
a  tailwind  component  in  order  to  minimize  an  already  high 
relative  sinking  speed. 

In  Fig.  5,  a  comparison  is  made  of  the  maximum  possible 
inverse  glide  slopes  between  the  700-  and  300-mb  levels  if  exact 


where  d„  =  (l/27)(2pi  -  9p„(}„  +  27r„):/„  =  (l/3)(39„  -p„^); 
n  =  1,2,3,... .  Again,  the  two  remaining  roots  to  Eq.  (17a)  are 
less  than  zero  and  consequently  have  no  physical  significance. 

Numerical  Results 

The  general  speed-to-fly  solutions  for  arbitrary  wind  speed 
and  direction  were  computed  between  the  700-  and  300-mb 
levels  by  seeded  iterations  using  Eqs.  (14-18).  The  computa¬ 
tions  were  based  on  an  unballasted  Nimbus  Ilb  glide  polar 
according  to  Eqs.  (8)  and  (9).  The  computational  sweeps  were 
stepped  in  both  the  positive  and  negative  wind  directions  using 
increments  of  6v  —  1.0  kt.  This  procedure  was  selected  in 
preference  to  some  general  iterative  search  routine  because  the 
speed-to-fly  equation  [Eq.  (13)]  has  multiple  roots  and  it  was 
desirable  to  find  cause  and  effect  relationships  between  speed- 
to-fly,  optimal  crab  angle,  air  mass  movements,  and  altitude 
effects  on  the  glide  polar. 

The  speed-to-fly  and  optimal  crab  angle  dependence  on  the 
wind  speed  and  direction  is  computed  in  Fig.  3  for  the  700-mb 
level  and  in  Fig.  4  for  the  300-mb  level.  These  computations 
assume  no  net  vertical  motion  in  the  atmosphere,  5  =  0,  as 
would  generally  be  the  case  during  long  glides  in  between  wave 
geiiesrating  mountains.  Negative  wind  speeds  correspond  to 
headwinds,  whereas  positive  values  denote  tailwinds.  Positive 
crab  angles  denote  deviations  from  the  course  line  in  the  same 
sense  as  the  wind  vector  (headwinds),  whereas  negative  crab 
angles  deviate  Trom  the  course  line  against  the  wind  angle 
(tailwinds).  The  fastest  and  slowest  speeds-to-fly  arise  from 
glides  directly  against  and  with  the  wind  (a  =  0),  respectively. 
The  speed-to-fly  in  a  direct  crosswind  (a  =  90  deg)  is  a  sym- 
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Fig.  4  Speed-to-fly  and  crab  angle  vs  winds  at  300  mb. 


Fig.  5  Maximum  inverse  glide  slopes  vs  winds  at  700  and  300  mb. 
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Fig.  6  Speed-to-fly  and  crab  angle  vs  subsidence  in  quartering  winds. 


adherence  to  speed-to-fly  and  optimal  crab  angle  is  made 
according  to  Figs.  3  and  4.  A  higher  L/D  is  possible  into  the 
wind  at  higher  altitudes  for  any  given  headwind  and  direction. 
Nonetheless,  strong  headwind  conditions  clearly  should  be 
avoided,  as  the  L/D  falls  catastrophically  from  47.7  in  still  air 
to  less  than  15.  Conversely,  flatter  glide  slopes  are  achieved 
for  any  given  tailwind  and  direction  at  lower  altitudes.  In  a 
direct  crosswind,  a  -  90  deg,  it  is  clearly  advantageous  to  fly 
at  a  higher  altitude  for  any  given  wind  speed.  Again,  the 
inverse  glide  slopes  for  60-kt  winds  at  700-mb  are  found  to  be 
comparable  to  those  in  90  kt  of  wind  at  the  300-mb  level. 
These  features  are  consistent  with  the  fact  that  a  higher  effec¬ 
tive  wing  loading  at  higher  altitudes  is  an  advantage  when 
penetrating  into  the  wind  or  in  strong  crosswinds  but  a  disad¬ 
vantage -when  gliding  slowly  at  minimum  relative  sinking 


speeds  with  the  wind.  Even  so,  a  90-kt  direct  tailwind  at  300 
mb  will  improve  the  L/V  from  47.7  in  still  air  to  92.7.  This 
value  degrades  to  65.0  if  the  300-mb  winds  become  quartering 
tailwinds,  a  =  45  deg. 

Figure  6  shows  computations  of  speed-to-fly  and  optimal 
crab  angles  for  the  case  of  quartering  headwinds  and  tailwinds 
at  the  500-mb  level  when  vertical  air  mass  movements  are 
encountered.  Contour  plots  are  shown  for  particular  air  mass 
sink  rates  of  -  1 .0  <  S  <  8.0  kt,  as  would  be  typically  encoun¬ 
tered  while  traversing  mountains  in  fields  of  lee  waves. 
Computations  for  an  air  mass  rising  faster  than  the  aircraft’s 
minimum  relative  sinking  speed  are  not  possible  because  a 
state  of  perpetual  motion  is  achieved  at  the  start  of  the  seeded 
iteration,  v  =  Vq,  for  which  the  problem  has  not  been  correctly 
posed  by  Eq.  (12).  The  problem  in  this  case  becomes  one  of 
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maximizing  climb  rate  for  which  the  speed  to  fly  is  U— Wj. 
Otherwise,  the  speed-to-fly  solutions  in  Fig.  6  require  one  to 
speed  up  in  sinking  air  and  to  slow  down  in  rising  air,  similar 
to  MacCready.*’^  However,  unlike  MacCready’s  original  re¬ 
sult,  the  speeds-to-fly  for  any  given  air  mass  sink  rate  are 
dependent  upon  the  wind  speed,  requiring  faster  glides  for  a 
headwind  component  and  slower  glides  with  a  tailwind  com¬ 
ponent.  Regardless  of  the  air  mass  sink  rate,  speeds-to-fly  in 
strong  quartering  tailwinds  tend  to  increase  slightly  with  in¬ 
creasing  wind  speed  in  order  to  prevent  glide  slope  degrada¬ 
tion  resulting  from  large  crab  angles,  )3>  0(45  deg).  This  result 
is  a  modification  to  Kuettner’s  original  downwind  strategy 
(see  Ref.  3).  In  strong  quartering  headwinds  with  a  rapidly 
sinking  air  mass,  one  finds  in  Figs.  6  that  rather  large  speeds 
to  fly,  0(200  kt),  are  required  that  exceed  the  structural  load 
limits  of  the  test  aircraft.  In  1-kt  rising  air  with  a  quartering 
tailwind,  the  speed-to-fly  is  less  than  the  best  L/D  speed  for 
still  air. 

Figure  7  shows  the  effect  of  vertical  air  mass  movements 
on  speeds-to-fly  and  optimal  crab  angles  in  direct  crosswinds 
at  the  500-mb  level.  In  strong  wind  with  rapidly  sinking  air, 
the  speeds-to-fly  are  significantly  less  than  the  corresponding 
quartering  headwind  cases  in  Fig.  6  but  also  significantly 
greater  than  the  corresponding  tailwind  cases.  For  weak  direct 
crosswinds  (v  <  30  kt)  there  is  little  change  in  speed  to  fly  with 
increasing  or  decreasing  wind  speed  since  most  of  the  compen¬ 
sation  is  done  with  the  crab  angle.  Unlike  all  other  cases 
considered,  the  optimal  crab  angle  for  strong  direct  cross- 
winds  (v  >  60  kt)  does  indeed  exceed  45  deg  in  nonsinking  air. 

Figure  8  shows  the  maximum  obtainable  inverse  glide  slopes 
with  vertical  air  mass  motion  in  quartering  headwinds  and 
tailwinds  and  in  direct  crosswinds  at  the  500-mb  level.  The 
computations  are  based  on  exact  adherence  to  speed  to  fly  and 
optimal  crab  angle  according  to  Figs.  6  and  7.  One  finds  that 
slowly  rising  air,  S  =  — 1,0  kt,  yields  impressive  L/D  for  both 
strong  quartering  tailwinds  {L/D  «  190)  as  well  as  weak  cross- 
winds  (L/D  « 140).  However,  strong  quartering  headwinds 
degrade  the  L/D  to  <10  regardless  of  the  vertical  air  mass 
motions.  In  fact,  the  inverse  glide  slope  in  this  case  is  almost 


as  bad  for  8  kt  of  air  mass  sink  as  it  would  be  for  ndnsinking 
air.  The  degradation  in  L/D  for  strong  direct  crosswinds  is 
only  slightly  less.  For  weak  quartering  headwinds/tailwinds 
and  direct  crosswinds,  any  rate  of  sinking  of  the  air  mass 
exerts  a  pronounced  reduction  in  L/D  compared  to  the  case 
S  =  0,  Most  of  this  performance  loss  occurs  during  the  first 
few  knots  of  air  mass  sink.  Very  little  additional  loss  occurs 
for  S  >  4.0  kt. 

Experiment 

To  ascertain  the  accuracy  of  the  speed-to-fly  and  optimal 
crab  angle  computations,  a  series  of  test  flights  were  con¬ 
ducted  with  an  unballasted  Nimbus  Ilb.  The  test  flights  were 
performed  using  the  active  lee  wave  generating  topography  of 
the  Laguna  Mountains  in  Southern  California.  Wave  entry 
was  accomplished  by  aerotow  to  a  point  in  the  neighborhood 
of  the  Julian  VORTAC.  Upon  release,  climb  was  initiated  to 
the  500-mb  level  according  to  visual  flight  rules,  with  addi¬ 
tional  climbs  to  the  300-mb  level  upon  receipt  of  clearance 
from  the  Los  Angeles  Center.  After  reaching  maximum  per¬ 
missible  altitude,  the  test  aircraft  departed  the  wave  generating 
region  for  a  series  of  way  points  whose  exact  position  was 
known.  The  altitude  lost  in  gliding  between  the  way  points  was 
recorded  by  means  of  a  precision  barograph  and  used  together 
with  the  known  distances  between  way  points  to  calculate  the 
inverse  glide  slope  with  respect  to  the  ground.  Passage  over 
each  way  point  was  determined  by  visual  recognition  and 
annotated  on  the  barograph  by  executing  a  rapid  dive  and 
recovery  maneuver  to  produce  a  mark  on  the  trace  known  as 
a  “notch.” 

Wind  data  for  the  tests  were  derived  from  National  Weather 
Service  soundings  at  San  Diego,  California,  and  Blythe,  Ari¬ 
zona.  Test  flights  were  conducted  about  the  700-  ,  500- ,  and 
300-mb  levels  for  which  sounding  data  and  forecast  charts 
were  available.  Way  points  were  selected  with  separations 
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Fig.  9  Effects  of  speed-to-fly  errors  in  quartering  tailwinds  and  di¬ 
rect  crosswinds. 


ranging  from  11.4  to  37.8  statute  miles  so  that  the  maximum 
altitude  variance  about  any  one  of  these  levels  did  not  exceed 
±  2000  ft.  Since  the  test  aircraft  lacked  accurate  positioning 
equipment,  no  attempt  was  made  to  test  speed-to-fly  theories 
within  the  variable  vertical  velocities  of  the  lee  wave  field. 
Instead,  the  way  points  were  selected  upwind  of  the  wave 
generating  topography  to  test  the  direct  crosswind  calculations 
and  a  considerable  distance  downwind  of  the  wave  generating 
topography  for  quartering  tailwind  tests.  This  arrangement 
was  aided  by  the  fact  that  the  lee  waves  were  generated  over  a 
broad  tilted  plateau  and  were,  therefore,  essentially  hydro¬ 
static  in  nature  with  vertically  upward  propagation  and  no 
secondary  or  tertiary  phases  downwind  (see  Ref.  5). 

Figure  9  shows  the  results  of  these  tests  during  which  known 
speed-to-fly  errors  bU  were  flown  to  yield  measured  inverse 
glide  slopes  L/£)(t/  =  6/ +  617).  The  speed-to-fly  errors  are 
corrected  for  the  airspeed  system  errors  measured  by  John¬ 
son.®  In  every  case  for  which  a  nonzero  speed-to-fly  error  was 
flown,  the  resulting  inverse  glide  slopes  are  found  to  be  less 
than  the  theoretical  maximum  L /Dmax  when  exact  adherence 
to  speed-to-fly  and  optimal  crab  angle  is  maintained.  The 
curves  in  Fig.  9  are  generated  by  a  numerical  search  about 
either  side  of  the  speed-to-fly  solutions  computed  from  Eqs. 
(14-18)  for  the  respective  altitudes  and  winds.  The  particular 
wind  speeds  for  those  computations  are  selected  in  between 


the  San  Diego  and  Blythe  soundings  based  on  a  minimization 
of  the  variance  between  the  data  and  resulting  theory  curves. 
A  check  of  the  L/Z>max  values  is  provided  by  the  data  points 
gathered  for  zero  speed-to-fly  errors.  The  loss  in  L/D  due  to 
any  given  speed-to-fly  error  is  greater  at  lower  altitudes  for 
both  quartering  tailwinds  and  direct  crosswinds.  Furthermore, 
flying  too  slowly  produces  a  greater  degradation  in  L/D  than 
flying  too  fast.  Speed-to-fly  errors  have  a  greater  adverse 
impact  in  quartering  tailwinds  than  in  direct  crosswinds.  How¬ 
ever,  in  either  case,  the  speed-to-fly  errors  resulting  from 
flying  too  slowly  can  diminish  the  L/D  by  15  points  or  more. 
This  is  equivalent  to  turning  the  clock  on  aerodynamic  design 
back  by  as  much  as  50  years.  In  other  words,  pilotage  must  be 
considered  a  leading-order  process. 

Conclusions 

1)  Speed-to-fly  during  glides  skewed  relative  to  the  wind  by 
less  than  90  deg  requires  flying  faster  through  sinking  air  or 
into  a  headwind  component  and  flying  slower  through  rising 
air  or  with  a  tailwind  component  such  that  the  crab  angle  does 
not  become  excessively  large,  j3<G(45  deg). 

2)  Speed-to-fly  in  direct  crosswinds  requires  flying  faster 
through  sinking  air  or  with  increasing  wind  speed  and  flying 
slower  through  rising  air  or  with  decreasing  wind  speed,  using 
crab  angles  that  may  exceed  45  deg  for  the  case  of  strong 
winds  in  a  nonsinking  air  mass. 

3)  Maximum  L/D  with  any  given  tailwind  component  is 
achieved  at  lower  altitudes. 

4)  Maximum  L/D  with  any  given  headwind  component  or 
any  given  direct  crosswind  is  achieved  at  higher  altitudes. 

5)  Speed-to-fly  errors  resulting  from  flying  too  slowly  yield 
the  greatest  losses  in  L  /D. 

6)  Losses  in  L/D  for  any  given  speed-to-fly  error  are 
greater  at  lower  altitudes. 
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